25char op_sxpun_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_sxpun.C,v 1.6 2014/10/13 08:53:26 j_novak Exp $" ;
99void _sxpun_pas_prevu(Tbl * tb,
int& base) {
100 cout <<
"sxpun pas prevu..." << endl ;
101 cout <<
"Tbl: " << tb <<
" base: " << base << endl ;
110void _sxpun_identite(Tbl* ,
int& ) {
118void _sxpun_r_jaco02(Tbl* tb,
int&)
121 if (tb->get_etat() == ETATZERO) {
127 int nr = (tb->dim).dim[0] ;
128 int nt = (tb->dim).dim[1] ;
129 int np = (tb->dim).dim[2] ;
133 double* xo =
new double [tb->get_taille()];
136 for (
int i=0; i<tb->get_taille(); i++) {
145 int borne_phi = np + 1 ;
150 for (
int k=0 ; k< borne_phi ; k++)
156 for (
int j=0 ; j<nt ; j++) {
160 for (
int i = 0 ; i < nr-1 ; i++) {
162 for (
int m = i+1 ; m < nr ; m++) {
163 int signe = ((m-1-i)%2 == 0 ? 1 : -1) ;
164 somme += signe*((m+1)*(m+2)/double((i+1)*(i+2))-(i+1)*(i+2)/double((m+1)*(m+2)))*xci[m] ;
166 xco[i] = (2*i+3)/
double(4)*somme ;