24char chb_legmp_cos_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legmp_cos.C,v 1.3 2014/10/13 08:53:11 j_novak Exp $" ;
116void chb_legmp_cos(
const int* deg ,
const double* cfi,
double* cfo) {
128 double* som =
new double[nr] ;
131 double* bb = mat_legmp_cos(np, nt) ;
138 const double* cc = cfi ;
156 for (j=0; j<nt; j++) {
159 for (i=0; i<nr; i++) {
163 for (l=m; l<nt; l++) {
165 double bmjl = bb[nt*j + l] ;
166 for (i=0; i<nr; i++) {
167 som[i] += bmjl * cc[nr*l + i] ;
171 for (i=0; i<nr; i++) {
181 for (i=ntnr; i<3*ntnr; i++) {
199 for (m=0; m < np + 1 ; m+=2) {
201 for (k2=0; k2 < 2; k2++) {
203 if ( (k == 1) || (k == np+1) ) {
205 for (j=0; j<nt; j++) {
206 for (i=0; i<nr; i++) {
216 for (j=0; j<nt; j++) {
219 for (i=0; i<nr; i++) {
223 for (l=m; l<nt; l++) {
225 double bmjl = bb[nt*j + l] ;
226 for (i=0; i<nr; i++) {
227 som[i] += bmjl * cc[nr*l + i] ;
231 for (i=0; i<nr; i++) {
253 assert(resu == cfo + (np+2)*ntnr) ;