23char chb_legpi_sini_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legpi_sini.C,v 1.6 2014/10/13 08:53:11 j_novak Exp $" ;
133void chb_legpi_sini(
const int* deg ,
const double* cfi,
double* cfo) {
143 assert( cfi != cfo ) ;
146 double* som =
new double[nr] ;
149 double* bb = mat_legpi_sini(np, nt) ;
156 const double* cc = cfi ;
169 for (j=0; j<nt-1; j++) {
172 for (i=0; i<nr; i++) {
176 for (l=0; l<nt-1; l++) {
177 double bmjl = bb[nt*j + l] ;
178 for (i=0; i<nr; i++) {
179 som[i] += bmjl * cc[nr*l + i] ;
183 for (i=0; i<nr; i++) {
191 for (i=0; i<nr; i++) {
199 for (i=0; i<2*ntnr; i++) {
214 for (l=0; l<nt; l++) {
215 for (i=0; i<nr; i++) {
230 for (j=0; j<nt-1; j++) {
233 for (i=0; i<nr; i++) {
237 for (l=0; l<nt-1; l++) {
238 double bmjl = bb[nt*j + l] ;
239 for (i=0; i<nr; i++) {
240 som[i] += bmjl * cc[nr*l + i] ;
244 for (i=0; i<nr; i++) {
252 for (i=0; i<nr; i++) {
267 for (m=3; m < np ; m+=2) {
269 for (k2=0; k2 < 2; k2++) {
273 for (j=0; j<nt-1; j++) {
276 for (i=0; i<nr; i++) {
280 for (l=(m-1)/2; l<nt-1; l++) {
281 double bmjl = bb[nt*j + l] ;
282 for (i=0; i<nr; i++) {
283 som[i] += bmjl * cc[nr*l + i] ;
287 for (i=0; i<nr; i++) {
295 for (i=0; i<nr; i++) {
315 for (l=0; l<nt; l++) {
316 for (i=0; i<nr; i++) {
324 assert(resu == cfo + (np+2)*ntnr) ;