NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ Integrator() [2/2]

template<typename T = double>
integration::Integrator< T >::Integrator ( const Integrator< T > &  I)
inline

Definition at line 199 of file integration.hpp.

References integration::I().

199  {
200  int k1;
201  k_=I.k_;
202  k1=k_+1;
203  poly_interpolation_= new T [k1*k1];
204  poly_differentiation_= new T [k1*k1];
205  poly_integration_= new T [k1*k1*k1];
206  bd_weights_= new T [k1+1];
207  gregory_weights_= new T [ 4*k1*k1 ];
208  gregory_omega_=gregory_weights_ + ((2*k_+1)*2*k1); /*pointer to last line*/
209  if(k_>1) rcorr_ = new T [ (k_-1)*k1*k1] ; else rcorr_=0;
210  integration::read_poly_interpolation(k_,poly_interpolation_);
211  integration::read_poly_differentiation(k_,poly_differentiation_);
212  integration::read_poly_integration(k_,poly_integration_);
213  integration::read_bd_weights(k_+1,bd_weights_);
214  integration::read_gregory_weights(k_,gregory_weights_);
215  integration::read_rcorr(k_,rcorr_);
216 }
Integrator< T > & I(int k)
+ Here is the call graph for this function: