28 char diff_xdsdx2_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Diff/diff_xdsdx2.C,v 1.4 2014/10/13 08:52:51 j_novak Exp $" ;
59 void multxpun_1d(
int,
double**,
int) ;
101 for (indice =0; indice<
max_points; indice++) {
102 if (nr_done[indice] ==
npoints) {
106 if (nr_done[indice] == -1)
111 cerr <<
"Diff_xdsdx2::get_matrice() : no space left!!" <<
'\n'
112 <<
"The value of Diff.max_points must be increased..." << endl ;
120 double* vect =
new double[
npoints] ;
121 double* cres =
new double[
npoints] ;
122 for (
int i=0; i<
npoints; i++) {
128 mult_xm1_1d_cheb(
npoints, vect, cres) ;
130 resu.
set(j,i) = cres[j] ;
135 resu.
set(j,i) = vect[j] ;
140 resu.
set(j,i) = vect[j] ;
151 ost <<
" x d2 / dx2 " << endl ;
virtual ~Diff_xdsdx2()
Destructor.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
static const int max_points
Maximal number of matrices stored per base.
virtual const Matrice & get_matrice() const
Returns the matrix associated with the operator.
int npoints
Number of coefficients.
void operator=(const Diff &)
Assignment to another Diff.
#define R_JACO02
base de Jacobi(0,2) ordinaire (finjac)
#define R_CHEBP
base de Cheb. paire (rare) seulement
#define R_CHEBI
base de Cheb. impaire (rare) seulement
#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.
Class for the elementary differential operator (see the base class Diff ).
void initialize()
Initializes arrays.
Diff_xdsdx2(int base_r, int nr)
Standard constructor.
void operator=(const Diff_xdsdx2 &)
Assignment to another Diff_xdsdx2.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
#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