99 if (nr_done[indice] ==
npoints) {
103 if (nr_done[indice] == -1)
108 cerr <<
"Diff_dsdx2::get_matrice() : no space left!!" <<
'\n' 109 <<
"The value of Diff.max_points must be increased..." << endl ;
117 double* vect =
new double[
npoints] ;
118 for (
int i=0; i<
npoints; i++) {
124 resu.
set(j,i) = vect[j] ;
134 ost <<
" d2 / dx2 " << endl ;
Class for the elementary differential operator (see the base class Diff ).
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Base (abstract) class for 1D spectral differential operators in one domain.
static const int max_points
Maximal number of matrices stored per base.
void operator=(const Diff_dsdx2 &)
Assignment to another Diff_dsdx2.
int npoints
Number of coefficients.
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
int base
Base in radial direction.
virtual ~Diff_dsdx2()
Destructor.
virtual const Matrice & get_matrice() const
Returns the matrix associated with the operator.
void operator=(const Diff &)
Assignment to another Diff.
void initialize()
Initializes arrays.
double & set(int j, int i)
Read/write of a particuliar element.
Diff_dsdx2(int base_r, int nr)
Standard constructor.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
#define MAX_BASE
Nombre max. de bases differentes.