NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ 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().

165  {
166  int k1;
167  if( this->k_==I.k_ ) return *this;
168  delete [] poly_interpolation_;
169  delete [] poly_differentiation_;
170  delete [] poly_integration_;
171  delete [] bd_weights_;
172  delete [] gregory_weights_;
173  delete [] rcorr_;
174  k_=I.k_;
175  k1=k_+1;
176  poly_interpolation_= new T [k1*k1];
177  poly_differentiation_= new T [k1*k1];
178  poly_integration_= new T [k1*k1*k1];
179  bd_weights_= new T [k1+1];
180  gregory_weights_= new T [ 4*k1*k1 ];
181  gregory_omega_=gregory_weights_ + ((2*k_+1)*2*k1); /*pointer to last line*/
182  if(k_>1) rcorr_ = new T [ (k_-1)*k1*k1] ; else rcorr_=0;
183  integration::read_poly_interpolation(k_,poly_interpolation_);
184  integration::read_poly_differentiation(k_,poly_differentiation_);
185  integration::read_poly_integration(k_,poly_integration_);
186  integration::read_bd_weights(k_+1,bd_weights_);
187  integration::read_gregory_weights(k_,gregory_weights_);
188  integration::read_rcorr(k_,rcorr_);
189  return *this;
190 }
Integrator< T > & I(int k)
+ Here is the call graph for this function: