24char mult2_xm1_1d_cheb_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/mult2_xm1_1d_cheb.C,v 1.2 2014/10/13 08:53:24 j_novak Exp $" ;
81void mult2_xm1_1d_cheb(
int nr,
const double* cf,
double* cresu) {
93 cresu[0] = ai*cf[0] + aip1*cf[1] + aip2*cf[2] ;
97 cresu[1] = double(-2)*cf[0] + double(1.75)*cf[1] + aip1*cf[2] + aip2*cf[3] ;
101 cresu[2] = double(0.5)*cf[0] + aim1*cf[1] + ai*cf[2] + aip1*cf[3]
107 for (i=3; i<nr-2; i++) {
108 cresu[i] = aim2*cf[i-2] + aim1*cf[i-1] + ai*cf[i] + aip1*cf[i+1]
114 cresu[nr-2] = aim2*cf[nr-4] + aim1*cf[nr-3] + ai*cf[nr-2] + aip1*cf[nr-1] ;
118 cresu[nr-1] = aim2*cf[nr-3] + aim1*cf[nr-2] + ai*cf[nr-1] ;