28 char diff_mx2_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Diff/diff_mx2.C,v 1.4 2014/10/13 08:52:50 j_novak Exp $" ;
59 void mult2_xp1_1d(
int,
double **,
int) ;
100 if (nr_done[indice] ==
npoints) {
104 if (nr_done[indice] == -1)
109 cerr <<
"Diff_mx2::get_matrice() : no space left!!" <<
'\n'
110 <<
"The value of Diff.max_points must be increased..." << endl ;
118 double* vect =
new double[
npoints] ;
119 double* cres =
new double[
npoints] ;
120 for (
int i=0; i<
npoints; i++) {
125 mult2_xm1_1d_cheb(
npoints, vect, cres) ;
127 resu.
set(j,i) = cres[j] ;
132 resu.
set(j,i) = vect[j] ;
137 resu.
set(j,i) = vect[j] ;
149 ost <<
" Multiplication by x^2 " << endl ;
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
static const int max_points
Maximal number of matrices stored per base.
Class for the elementary differential operator multiplication by (see the base class Diff ).
int npoints
Number of coefficients.
void initialize()
Initializes arrays.
void operator=(const Diff &)
Assignment to another Diff.
#define R_JACO02
base de Jacobi(0,2) ordinaire (finjac)
virtual const Matrice & get_matrice() const
Returns the matrix associated with the operator.
Diff_mx2(int base_r, int nr)
Standard constructor.
#define TRA_R
Translation en R, used for a bitwise shift (in hex)
Base (abstract) class for 1D spectral differential operators in one domain.
int base
Base in radial direction.
virtual ~Diff_mx2()
Destructor.
void operator=(const Diff_mx2 &)
Assignment to another Diff_mx2.
#define MAX_BASE
Nombre max. de bases differentes.
double & set(int j, int i)
Read/write of a particuliar element.
#define R_CHEBU
base de Chebychev ordinaire (fin), dev. en 1/r