23char chb_legi_cossinci_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legi_cossinci.C,v 1.5 2014/10/13 08:53:10 j_novak Exp $" ;
130void chb_legi_cossinci(
const int* deg ,
const double* cfi,
double* cfo) {
132int ip, k2, l, j, i, m ;
142 double* som =
new double[nr] ;
145 double* bb = mat_legi_cossinci(np, nt) ;
152 const double* cc = cfi ;
168 for (ip=0; ip < np/4 + 1 ; ip++) {
175 for (k2=0; k2 < 2; k2++) {
177 if ( (k == 1) || (k == np+1) ) {
179 for (j=0; j<nt; j++) {
180 for (i=0; i<nr; i++) {
191 for (j=0; j<nt-1; j++) {
193 for (i=0; i<nr; i++) {
197 for (l=m/2; l<nt-1; l++) {
199 double bmjl = bb[nt*j + l] ;
200 for (i=0; i<nr; i++) {
201 som[i] += bmjl * cc[nr*l + i] ;
205 for (i=0; i<nr; i++) {
213 for (i=0; i<nr; i++) {
234 for (k2=0; k2 < 2; k2++) {
238 for (j=0; j<nt; j++) {
239 for (i=0; i<nr; i++) {
251 for (i=0; i<nr; i++) {
258 for (j=1; j<nt-1; j++) {
260 for (i=0; i<nr; i++) {
264 for (l=(m+1)/2; l<nt-1; l++) {
265 double bmjl = bb[nt*j + l] ;
266 for (i=0; i<nr; i++) {
267 som[i] += bmjl * cc[nr*l + i] ;
271 for (i=0; i<nr; i++) {
279 for (i=0; i<nr; i++) {
302 assert(resu == cfo + (np+2)*ntnr) ;