23char op_mult_xp1_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_mult_xp1.C,v 1.3 2014/10/13 08:53:26 j_novak Exp $" ;
79void _mult_xp1_pas_prevu(Tbl * tb,
int& base) {
80 cout <<
"mult_xp1 pas prevu..." << endl ;
81 cout <<
"Tbl: " << tb <<
" base: " << base << endl ;
90void _mult_xp1_identite(Tbl* ,
int& ) {
98void _mult_xp1_r_jaco02(Tbl* tb,
int& )
101 if (tb->get_etat() == ETATZERO) {
106 int nr = (tb->dim).dim[0] ;
107 int nt = (tb->dim).dim[1] ;
108 int np = (tb->dim).dim[2] ;
112 double* xo =
new double [tb->get_taille()];
115 for (
int i=0; i<tb->get_taille(); i++) {
124 int borne_phi = np + 1 ;
129 for (
int k=0 ; k< borne_phi ; k++)
135 for (
int j=0 ; j<nt ; j++) {
137 xco[0] = 1.5*xci[0] + 0.3*xci[1] ;
138 for (
int i = 1 ; i < nr-1 ; i++) {
139 xco[i] = i*(i+2)/
double((i+1)*(2*i+1))*xci[i-1] + (i*i+3*i+3)/
double((i+1)*(i+2))*xci[i] + (i+1)*(i+3)/double((i+2)*(2*i+5))*xci[i+1] ;
141 xco[nr-1] = (nr*nr-1)/
double((nr)*(2*nr-1))*xci[nr-2] + (1+1/double((nr)*(nr+1)))*xci[nr-1] ;