#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;
}
04-28