23char chb_legii_sinp_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_legii_sinp.C,v 1.4 2014/10/13 08:53:10 j_novak Exp $" ;
119void chb_legii_sinp(
const int* deg ,
const double* cfi,
double* cfo) {
129 assert( cfi != cfo ) ;
132 double* som =
new double[nr] ;
134 double* bb = mat_legii_sinp(np, nt) ;
141 const double* cc = cfi ;
153 for (i=0; i<nr; i++) {
160 for (j=1; j<nt-1; j++) {
163 for (i=0; i<nr; i++) {
167 for (l=1; l<nt-1; l++) {
168 double bmjl = bb[nt*j + l] ;
169 for (i=0; i<nr; i++) {
170 som[i] += bmjl * cc[nr*l + i] ;
174 for (i=0; i<nr; i++) {
182 for (i=0; i<nr; i++) {
190 for (i=0; i<2*ntnr; i++) {
205 for (l=0; l<nt; l++) {
206 for (i=0; i<nr; i++) {
220 for (i=0; i<nr; i++) {
227 for (j=1; j<nt-1; j++) {
230 for (i=0; i<nr; i++) {
234 for (l=1; l<nt-1; l++) {
235 double bmjl = bb[nt*j + l] ;
236 for (i=0; i<nr; i++) {
237 som[i] += bmjl * cc[nr*l + i] ;
241 for (i=0; i<nr; i++) {
249 for (i=0; i<nr; i++) {
264 for (m=3; m < np ; m+=2) {
266 for (k2=0; k2 < 2; k2++) {
271 for (i=0; i<nr; i++) {
276 for (j=1; j<nt-1; j++) {
278 for (i=0; i<nr; i++) {
282 for (l=(m+1)/2; l<nt-1; l++) {
283 double bmjl = bb[nt*j + l] ;
284 for (i=0; i<nr; i++) {
285 som[i] += bmjl * cc[nr*l + 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) ;