gsl库中的插值用法

#include <stdlib.h>
#include <stdio.h>
#include <iostream>
#include <sstream>
#include <algorithm>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_interp.h>
#include <gsl/gsl_spline.h>
#include <gsl/gsl_splinalg.h>
#include <gsl/gsl_spline2d.h>
#include <mutex>
#include <vector>
#include <fstream>

gsl_interp_accel* g_acc = gsl_interp_accel_alloc();
std::vector<double> g_time;
std::vector<double> g_lon;

void main()
{
	g_time.push_back(0);
	g_lon.push_back(0);

	g_time.push_back(1);
	g_lon.push_back(10);

	g_time.push_back(2);
	g_lon.push_back(20);

	g_time.push_back(3);
	g_lon.push_back(30);


	gsl_spline* spline
		= gsl_spline_alloc(gsl_interp_cspline, g_time.size());



	gsl_spline_init(spline, &g_time[0], &g_lon[0], g_time.size()); // init the splin
	 
	double r = gsl_spline_eval(spline, 2.5, g_acc);     //get the value from splin (at time=2.5)

	int xx = 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值