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);
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_);
Integrator< T > & I(int k)