24char mult_xm1_1d_cheb_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/mult_xm1_1d_cheb.C,v 1.2 2014/10/13 08:53:24 j_novak Exp $" ;
82void mult_xm1_1d_cheb(
int nr,
const double* cf,
double* cresu) {
92 cresu[0] = ai*cf[0] + aip1*cf[1] ;
96 cresu[1] = cf[0] + ai*cf[1] + aip1*cf[2] ;
102 for (i=2; i<nr-1; i++) {
103 cresu[i] = aim1*cf[i-1] + ai*cf[i] + aip1*cf[i+1] ;
108 cresu[nr-1] = aim1*cf[nr-2] + ai*cf[nr-1] ;