// spline test
#include <iostream>
#include <iomanip>
#include <vector>
// home made call to GSL spline
#include "gsl_spline.hpp"
using namespace std;
typedef vector<double> vec;


int main(){
  vec x(10),y(10);
  // known points (x,y)
  double t;
  t = -5.0;
  for(unsigned i=0;i<x.size();++i){
    x[i] = t; y[i] = 10.0-t*t;
    t += 1.0;
  }
  y[4] += 10.0; // one point up to demonstrate Akima spline
  cout<<fixed<<setprecision(8);
  for(unsigned i=0;i<x.size();++i)
    cout<<x[i]<<" "<<y[i]<<endl;
  cout<<"\n\n\n";
  
  // points to interpolate xx
  vec xx(100),yy(100);
  double dx = (x[x.size()-1]-x[0])/(xx.size()-1);
  for(unsigned i=0;i<xx.size();++i){ xx[i] = x[0]+i*dx;}

  // interpolate using natural cubic spline
  spline(x,y,xx,yy,"cspline");
  for(unsigned i=0;i<xx.size();++i)
    cout<<xx[i]<<" "<<yy[i]<<endl;
  cout<<"\n\n\n";
  // interpolate using natural Akima spline
  spline(x,y,xx,yy,"Akima");
  for(unsigned i=0;i<xx.size();++i)
    cout<<xx[i]<<" "<<yy[i]<<endl;
 
  return 0;
}