23char chb_cosi_legip_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_cosi_legip.C,v 1.7 2014/10/13 08:53:10 j_novak Exp $" ;
135void chb_cosi_legip(
const int* deg ,
const double* cfi,
double* cfo) {
137int k2, l, jmin, j, i, m ;
145 assert( cfi != cfo ) ;
148 double* som =
new double[nr] ;
151 double* aa = mat_cosi_legip(np, nt) ;
158 const double* cc = cfi ;
177 for (l=m/2; l<nt-1; l++) {
178 for (i=0; i<nr; i++) {
184 for (j=jmin; j<nt-1; 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++) {
199 for (i=0; i<nr; i++) {
207 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 (l=0; l<nt; l++) {
234 for (i=0; i<nr; i++) {
245 for (l=0; l<m/2; l++) {
246 for (i=0; i<nr; i++) {
252 for (l=m/2; l<nt-1; l++) {
253 for (i=0; i<nr; i++) {
259 for (j=jmin; j<nt-1; j++) {
260 double amlj = aa[nt*l + j] ;
261 for (i=0; i<nr; i++) {
262 som[i] += amlj * cc[nr*j + i] ;
266 for (i=0; i<nr; i++) {
274 for (i=0; i<nr; i++) {
294 assert(resu == cfo + (np+2)*ntnr) ;