: linear interpolation (Steffen) of data [x, y] More...

#include <core_power.h>

Inheritance diagram for Interp_obj:

Public Member Functions
	Interp_obj ()

	~Interp_obj ()

double	operator() (double x) const

template<typename T , size_t N>
void	init (const Data_Vec< T, N > &data)

Public Attributes
double	x_min

double	x_max

Private Attributes
bool	is_init

gsl_spline *	spline

gsl_interp_accel *	acc

Detailed Description

: linear interpolation (Steffen) of data [x, y]

Definition at line 96 of file core_power.h.

Constructor & Destructor Documentation

Interp_obj::Interp_obj ( )

inline

Definition at line 99 of file core_power.h.

References init(), operator()(), x, and ~Interp_obj().

99 : is_init(false) {}

Interp_obj::is_init

bool is_init

Definition: core_power.h:107

Interp_obj::~Interp_obj ( )

Definition at line 511 of file core_power.cpp.

Referenced by Interp_obj().

 {
     if (is_init)
     {
         gsl_spline_free(spline);
         gsl_interp_accel_free (acc);
     }
 }

Member Function Documentation

template<typename T , size_t N>

template void Interp_obj::init ( const Data_Vec< T, N > & data )

EXPLICIT INSTANTIATION *

Definition at line 496 of file core_power.cpp.

References Catch::begin, Catch::end, and Data_Vec< T, N >::size().

Referenced by Extrap_Pk< T, N >::Extrap_Pk(), gen_sigma_func_binned_gsl_qawf_nl(), and Interp_obj().

 {
     is_init = true;
     acc = gsl_interp_accel_alloc ();
     spline = gsl_spline_alloc (gsl_interp_steffen, data.size());
 
     std::vector<double> tmp_x(data[0].begin(), data[0].end());
     std::vector<double> tmp_y(data[1].begin(), data[1].end());
 
     gsl_spline_init (spline, tmp_x.data(), tmp_y.data(), data.size());
 
     x_min = data[0].front();
     x_max = data[0].back();
 }