23char chb_cossinci_legi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_cossinci_legi.C,v 1.5 2014/10/13 08:53:10 j_novak Exp $" ;
131void chb_cossinci_legi(
const int* deg ,
const double* cfi,
double* cfo) {
135int ip, k2, l, jmin, j, i, m ;
145 double* som =
new double[nr] ;
148 double* aa = mat_cossinci_legi(np, nt) ;
155 const double* cc = cfi ;
171 for (ip=0; ip < np/4 + 1 ; ip++) {
178 for (k2=0; k2 < 2; k2++) {
180 if ( (k == 1) || (k == np+1) ) {
182 for (l=0; l<nt; l++) {
183 for (i=0; i<nr; i++) {
194 for (l=0; l<m/2; l++) {
195 for (i=0; i<nr; i++) {
202 for (l=m/2; l<nt-1; l++) {
203 for (i=0; i<nr; i++) {
209 for (j=jmin; j<nt-1; j++) {
210 double amlj = aa[nt*l + j] ;
211 for (i=0; i<nr; i++) {
212 som[i] += amlj * cc[nr*j + i] ;
216 for (i=0; i<nr; i++) {
225 for (i=0; i<nr; i++) {
246 for (k2=0; k2 < 2; k2++) {
250 for (l=0; l<nt; l++) {
251 for (i=0; i<nr; i++) {
264 for (l=0; l<(m+1)/2; l++) {
265 for (i=0; i<nr; i++) {
272 for (l=(m+1)/2; l<nt-1; l++) {
273 for (i=0; i<nr; i++) {
280 for (j=jmin; j<nt-1; j++) {
281 double amlj = aa[nt*l + j] ;
282 for (i=0; i<nr; i++) {
283 som[i] += amlj * cc[nr*j + i] ;
287 for (i=0; i<nr; i++) {
295 for (i=0; i<nr; i++) {
319 assert(resu == cfo + (np+2)*ntnr) ;