23char chb_legp_cossincp_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legp_cossincp.C,v 1.7 2014/10/13 08:53:11 j_novak Exp $" ;
141void chb_legp_cossincp(
const int* deg ,
const double* cfi,
double* cfo) {
143int ip, k2, l, j, i, m ;
153 double* som =
new double[nr] ;
156 double* bb = mat_legp_cossincp(np, nt) ;
163 const double* cc = cfi ;
181 for (j=0; j<nt; j++) {
184 for (i=0; i<nr; i++) {
188 for (l=m/2; l<nt; l++) {
190 double bmjl = bb[nt*j + l] ;
191 for (i=0; i<nr; i++) {
192 som[i] += bmjl * cc[nr*l + i] ;
196 for (i=0; i<nr; i++) {
206 for (i=0; i<2*ntnr; i++) {
230 for (ip=0; ip < np/4 + 1 ; ip++) {
237 for (k2=0; k2 < 2; k2++) {
239 if ( (k == 1) || (k == np+1) ) {
241 for (j=0; j<nt; j++) {
242 for (i=0; i<nr; i++) {
252 for (j=0; j<nt; j++) {
255 for (i=0; i<nr; i++) {
259 for (l=m/2; l<nt; l++) {
261 double bmjl = bb[nt*j + l] ;
262 for (i=0; i<nr; i++) {
263 som[i] += bmjl * cc[nr*l + i] ;
267 for (i=0; i<nr; i++) {
290 for (k2=0; k2 < 2; k2++) {
294 for (j=0; j<nt; j++) {
295 for (i=0; i<nr; i++) {
306 for (j=0; j<nt-1; j++) {
309 for (i=0; i<nr; i++) {
313 for (l=(m-1)/2; l<nt-1; l++) {
314 double bmjl = bb[nt*j + l] ;
315 for (i=0; i<nr; i++) {
316 som[i] += bmjl * cc[nr*l + i] ;
320 for (i=0; i<nr; i++) {
328 for (i=0; i<nr; i++) {