23char chb_cossincp_legp_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_cossincp_legp.C,v 1.7 2014/10/13 08:53:10 j_novak Exp $" ;
135void chb_cossincp_legp(
const int* deg ,
const double* cfi,
double* cfo) {
139int ip, k2, l, jmin, j, i, m ;
149 double* som =
new double[nr] ;
152 double* aa = mat_cossincp_legp(np, nt) ;
159 const double* cc = cfi ;
178 for (l=m/2; l<nt; l++) {
179 for (i=0; i<nr; i++) {
184 for (j=jmin; j<nt; j++) {
185 double amlj = aa[nt*l + j] ;
186 for (i=0; i<nr; i++) {
187 som[i] += amlj * cc[nr*j + i] ;
191 for (i=0; i<nr; i++) {
201 for (i=0; i<2*ntnr; i++) {
226 for (ip=0; ip < np/4 + 1 ; ip++) {
233 for (k2=0; k2 < 2; k2++) {
235 if ( (k == 1) || (k == np+1) ) {
237 for (l=0; l<nt; l++) {
238 for (i=0; i<nr; i++) {
248 for (l=0; l<m/2; l++) {
249 for (i=0; i<nr; i++) {
255 for (l=m/2; l<nt; l++) {
256 for (i=0; i<nr; i++) {
260 jmin = ( m == 0 ) ? l : 0 ;
261 for (j=jmin; j<nt; j++) {
262 double amlj = aa[nt*l + j] ;
263 for (i=0; i<nr; i++) {
264 som[i] += amlj * cc[nr*j + i] ;
268 for (i=0; i<nr; i++) {
291 for (k2=0; k2 < 2; k2++) {
295 for (l=0; l<nt; l++) {
296 for (i=0; i<nr; i++) {
306 for (l=0; l<(m-1)/2; l++) {
307 for (i=0; i<nr; i++) {
314 for (l=(m-1)/2; l<nt-1; l++) {
315 for (i=0; i<nr; i++) {
319 jmin = ( m == 1 ) ? l : 0 ;
321 for (j=jmin; j<nt-1; j++) {
322 double amlj = aa[nt*l + j] ;
323 for (i=0; i<nr; i++) {
324 som[i] += amlj * cc[nr*j + i] ;
328 for (i=0; i<nr; i++) {
336 for (i=0; i<nr; i++) {