24char chb_sin_legmi_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_sin_legmi.C,v 1.3 2014/10/13 08:53:11 j_novak Exp $" ;
114void chb_sin_legmi(
const int* deg ,
const double* cfi,
double* cfo) {
116int k2, l, jmin, j, i, m ;
124 assert( cfi != cfo ) ;
127 double* som =
new double[nr] ;
130 double* aa = mat_sin_legmi(np, nt) ;
137 const double* cc = cfi ;
149 for (i=0; i<nr; i++) {
155 for (l=1; l<nt-1; l++) {
156 for (i=0; i<nr; i++) {
162 for (j=jmin; j<nt-1; j++) {
163 double amlj = aa[nt*l + j] ;
164 for (i=0; i<nr; i++) {
165 som[i] += amlj * cc[nr*j + i] ;
169 for (i=0; i<nr; i++) {
177 for (i=0; i<nr; i++) {
185 for (i=0; i<2*ntnr; i++) {
201 for (l=0; l<nt; l++) {
202 for (i=0; i<nr; i++) {
216 for (i=0; i<nr; i++) {
222 for (l=1; l<nt-1; l++) {
223 for (i=0; i<nr; i++) {
229 for (j=jmin; j<nt-1; j++) {
230 double amlj = aa[nt*l + j] ;
231 for (i=0; i<nr; i++) {
232 som[i] += amlj * cc[nr*j + i] ;
236 for (i=0; i<nr; i++) {
244 for (i=0; i<nr; i++) {
259 for (m=3; m < np ; m+=2) {
261 for (k2=0; k2 < 2; k2++) {
262 int lmax = (m<nt-1 ? m : nt -1) ;
263 for (l=0; l<lmax; l++) {
264 for (i=0; i<nr; i++) {
271 for (l=m; l<nt-1; l++) {
272 for (i=0; i<nr; i++) {
278 for (j=jmin; j<nt-1; j++) {
279 double amlj = aa[nt*l + j] ;
280 for (i=0; i<nr; i++) {
281 som[i] += amlj * cc[nr*j + i] ;
285 for (i=0; i<nr; i++) {
293 for (i=0; i<nr; i++) {
313 for (l=0; l<nt; l++) {
314 for (i=0; i<nr; i++) {
322 assert(resu == cfo + (np+2)*ntnr) ;