23char chb_sini_legpi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_sini_legpi.C,v 1.6 2014/10/13 08:53:11 j_novak Exp $" ;
130void chb_sini_legpi(
const int* deg ,
const double* cfi,
double* cfo) {
132int k2, l, jmin, j, i, m ;
140 assert( cfi != cfo ) ;
143 double* som =
new double[nr] ;
146 double* aa = mat_sini_legpi(np, nt) ;
153 const double* cc = cfi ;
165 for (l=0; l<nt-1; l++) {
166 for (i=0; i<nr; i++) {
172 for (j=jmin; j<nt-1; j++) {
173 double amlj = aa[nt*l + j] ;
174 for (i=0; i<nr; i++) {
175 som[i] += amlj * cc[nr*j + i] ;
179 for (i=0; i<nr; i++) {
187 for (i=0; i<nr; i++) {
195 for (i=0; i<2*ntnr; i++) {
211 for (l=0; l<nt; l++) {
212 for (i=0; i<nr; i++) {
226 for (l=0; l<nt-1; l++) {
227 for (i=0; i<nr; i++) {
233 for (j=jmin; j<nt-1; j++) {
234 double amlj = aa[nt*l + j] ;
235 for (i=0; i<nr; i++) {
236 som[i] += amlj * cc[nr*j + i] ;
240 for (i=0; i<nr; i++) {
248 for (i=0; i<nr; i++) {
263 for (m=3; m < np ; m+=2) {
265 for (k2=0; k2 < 2; k2++) {
267 for (l=0; l<(m-1)/2; l++) {
268 for (i=0; i<nr; i++) {
275 for (l=(m-1)/2; l<nt-1; l++) {
276 for (i=0; i<nr; i++) {
282 for (j=jmin; j<nt-1; j++) {
283 double amlj = aa[nt*l + j] ;
284 for (i=0; i<nr; i++) {
285 som[i] += amlj * cc[nr*j + i] ;
289 for (i=0; i<nr; i++) {
297 for (i=0; i<nr; i++) {
317 for (l=0; l<nt; l++) {
318 for (i=0; i<nr; i++) {
326 assert(resu == cfo + (np+2)*ntnr) ;