23char chb_leg_cossinc_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_leg_cossinc.C,v 1.4 2014/10/13 08:53:10 j_novak Exp $" ;
128void chb_leg_cossinc(
const int* deg ,
const double* cfi,
double* cfo) {
130int ip, k2, l, j, i, m ;
140 double* som =
new double[nr] ;
143 double* bb = mat_leg_cossinc(np, nt) ;
170 const double* cc = cfi ;
186 for (ip=0; ip < np/4 + 1 ; ip++) {
193 for (k2=0; k2 < 2; k2++) {
195 if ( (k == 1) || (k == np+1) ) {
197 for (j=0; j<nt; j++) {
198 for (i=0; i<nr; i++) {
208 for (j=0; j<nt; j++) {
211 for (i=0; i<nr; i++) {
215 for (l=m; l<nt; l++) {
217 double bmjl = bb[nt*j + l] ;
218 for (i=0; i<nr; i++) {
219 som[i] += bmjl * cc[nr*l + i] ;
223 for (i=0; i<nr; i++) {
246 for (k2=0; k2 < 2; k2++) {
250 for (j=0; j<nt; j++) {
251 for (i=0; i<nr; i++) {
262 for (j=0; j<nt-1; j++) {
265 for (i=0; i<nr; i++) {
269 for (l=m; l<nt-1; l++) {
270 double bmjl = bb[nt*j + l] ;
271 for (i=0; i<nr; i++) {
272 som[i] += bmjl * cc[nr*l + i] ;
276 for (i=0; i<nr; i++) {
284 for (i=0; i<nr; i++) {