Integrator() [1/2]

template<typename T = double>

integration::Integrator< T >::Integrator ( int  k = 0 )

inline

Initializes the Integrator class for a given order k.

Purpose

Initializes the Integrator class for given order \(k\) between 1 and 8. If \(k>5\), Gregory integration up to order \(k\) is used, while all other operations are limited to maximum order 5.

Parameters

k	The order of the weights.

Definition at line 147 of file integration.hpp.

References integration::Integrator< T >::k().

                     {
  int k1=k+1;
  if(k<0 || k > GREGORY_KMAX ){ std::cout << "Integrator: k out of range " << std::endl; abort();}
  poly_interpolation_= new T [k1*k1];
  poly_differentiation_= new T [k1*k1];
  poly_integration_= new T [k1*k1*k1];
  bd_weights_= new T [k1+1]; 
  gregory_weights_= new T [ 4*k1*k1 ];
  gregory_omega_=gregory_weights_ + ((2*k+1)*2*k1); /*pointer to last line*/
  if(k>1) rcorr_ = new T [ (k-1)*k1*k1] ; else rcorr_=0;
  k_=k;
  integration::read_poly_interpolation(k_,poly_interpolation_); 
  integration::read_poly_differentiation(k_,poly_differentiation_); 
  integration::read_poly_integration(k_,poly_integration_); 
  integration::read_bd_weights(k+1,bd_weights_);
  integration::read_gregory_weights(k,gregory_weights_); 
  integration::read_rcorr(k_,rcorr_);
}    

Here is the call graph for this function: