23char chb_cossinc_leg_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_cossinc_leg.C,v 1.5 2014/10/13 08:53:10 j_novak Exp $" ;
124void chb_cossinc_leg(
const int* deg ,
const double* cfi,
double* cfo) {
128int ip, k2, l, jmin, j, i, m ;
138 double* som =
new double[nr] ;
141 double* aa = mat_cossinc_leg(np, nt) ;
167 const double* cc = cfi ;
183 for (ip=0; ip < np/4 + 1 ; ip++) {
190 for (k2=0; k2 < 2; k2++) {
192 if ( (k == 1) || (k == np+1) ) {
194 for (l=0; l<nt; l++) {
195 for (i=0; i<nr; i++) {
205 for (l=0; l<m; l++) {
206 for (i=0; i<nr; i++) {
212 for (l=m; l<nt; l++) {
213 for (i=0; i<nr; i++) {
217 jmin = ( m == 0 ) ? l : 0 ;
218 for (j=jmin; j<nt; j++) {
219 double amlj = aa[nt*l + j] ;
220 for (i=0; i<nr; i++) {
221 som[i] += amlj * cc[nr*j + i] ;
225 for (i=0; i<nr; i++) {
248 for (k2=0; k2 < 2; k2++) {
252 for (l=0; l<nt; l++) {
253 for (i=0; i<nr; i++) {
263 for (l=0; l<m; l++) {
264 for (i=0; i<nr; i++) {
271 for (l=m; l<nt-1; l++) {
272 for (i=0; i<nr; i++) {
276 jmin = ( m == 1 ) ? l : 0 ;
278 for (j=jmin; j<nt-1; j++) {
279 double amlj = aa[nt*l + j] ;
280 for (i=0; i<nr; i++) {
281 som[i] += amlj * cc[nr*j + i] ;
286 for (i=0; i<nr; i++) {
294 for (i=0; i<nr; i++) {