operator=()

template<typename T = double>

Integrator& integration::Integrator< T >::operator= ( const Integrator< T > &  I )

inline

Definition at line 165 of file integration.hpp.

References integration::I().

                                            {
  int k1;
  if( this->k_==I.k_ ) return *this;
  delete [] poly_interpolation_;
  delete [] poly_differentiation_;
  delete [] poly_integration_;
  delete [] bd_weights_;
  delete [] gregory_weights_;
  delete [] rcorr_;
  k_=I.k_;
  k1=k_+1;
  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;
  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_);
  return *this;
}

Here is the call graph for this function: