/* interpolation.h */

#ifndef INTERPOLATION_H_
#define INTERPOLATION_H_

#include <iostream>
#include <gsl/gsl_spline.h>
#include <memory>
#include <vector>

class interpolation
{
	public:
		interpolation();
		interpolation(std::vector<double>, std::vector<double>);

		auto set(std::vector<double>, std::vector<double>)	-> void;
		auto operator()(double x)														-> double;

		auto min()-> double { return _x_min; }
		auto max()-> double { return _x_max; }


	private:
		std::vector<double> _x, _y;
		double _x_min, _x_max;
		std::unique_ptr<gsl_interp_accel,void(*)(gsl_interp_accel*)> _acc;
		std::unique_ptr<gsl_spline,void(*)(gsl_spline*)> _spline;

};

#endif

/* End of File */