23char op_mult_xm1_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/op_mult_xm1.C,v 1.2 2014/10/13 08:53:26 j_novak Exp $" ;
69#include "type_parite.h"
80void _mult_xm1_identite(Tbl* ,
int& ) {
88void _mult_xm1_cheb(Tbl *tb,
int& )
92 if (tb->get_etat() == ETATZERO) {
97 int nr = (tb->dim).dim[0] ;
98 int nt = (tb->dim).dim[1] ;
99 int np = (tb->dim).dim[2] ;
104 double* trav =
new double[nr] ;
107 for (k=0 ; k<np+1 ; k++) {
109 for (j=0 ; j<nt ; j++) {
111 double* cf = tb->t + k*ntnr + j*nr ;
113 mult_xm1_1d_cheb(nr, cf, trav) ;
115 for (i=0; i<nr; i++) {