23char op_mult2_xm1_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_mult2_xm1.C,v 1.2 2014/10/13 08:53:25 j_novak Exp $" ;
66#include "type_parite.h"
77void _mult2_xm1_identite(Tbl* ,
int& ) {
85void _mult2_xm1_cheb(Tbl *tb,
int& )
89 if (tb->get_etat() == ETATZERO) {
94 int nr = (tb->dim).dim[0] ;
95 int nt = (tb->dim).dim[1] ;
96 int np = (tb->dim).dim[2] ;
101 double* trav =
new double[nr] ;
104 for (k=0 ; k<np+1 ; k++) {
106 for (j=0 ; j<nt ; j++) {
108 double* cf = tb->t + k*ntnr + j*nr ;
110 mult2_xm1_1d_cheb(nr, cf, trav) ;
112 for (i=0; i<nr; i++) {