23char chb_legip_cosi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legip_cosi.C,v 1.7 2014/10/13 08:53:10 j_novak Exp $" ;
137void chb_legip_cosi(
const int* deg ,
const double* cfi,
double* cfo) {
147 assert( cfi != cfo ) ;
150 double* som =
new double[nr] ;
153 double* bb = mat_legip_cosi(np, nt) ;
160 const double* cc = cfi ;
178 for (j=0; j<nt-1; j++) {
181 for (i=0; i<nr; i++) {
185 for (l=m/2; l<nt-1; l++) {
187 double bmjl = bb[nt*j + l] ;
188 for (i=0; i<nr; i++) {
189 som[i] += bmjl * cc[nr*l + i] ;
193 for (i=0; i<nr; i++) {
201 for (i=0; i<nr; i++) {
209 for (i=ntnr; i<3*ntnr; i++) {
227 for (m=0; m < np + 1 ; m+=2) {
229 for (k2=0; k2 < 2; k2++) {
231 if ( (k == 1) || (k == np+1) ) {
233 for (j=0; j<nt; j++) {
234 for (i=0; i<nr; i++) {
244 for (j=0; j<nt-1; j++) {
247 for (i=0; i<nr; i++) {
251 for (l=m/2; l<nt-1; l++) {
253 double bmjl = bb[nt*j + l] ;
254 for (i=0; i<nr; i++) {
255 som[i] += bmjl * cc[nr*l + i] ;
259 for (i=0; i<nr; i++) {
267 for (i=0; i<nr; i++) {
287 assert(resu == cfo + (np+2)*ntnr) ;