23char chb_cosp_legpp_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/chb_cosp_legpp.C,v 1.7 2014/10/13 08:53:10 j_novak Exp $" ;
134void chb_cosp_legpp(
const int* deg ,
const double* cfi,
double* cfo) {
136int k2, l, jmin, j, i, m ;
146 double* som =
new double[nr] ;
149 double* aa = mat_cosp_legpp(np, nt) ;
156 const double* cc = cfi ;
175 for (l=m/2; l<nt; l++) {
176 for (i=0; i<nr; i++) {
181 for (j=jmin; j<nt; j++) {
182 double amlj = aa[nt*l + j] ;
183 for (i=0; i<nr; i++) {
184 som[i] += amlj * cc[nr*j + i] ;
188 for (i=0; i<nr; i++) {
198 for (i=ntnr; i<3*ntnr; i++) {
218 for (m=0; m < np + 1 ; m+=2) {
220 for (k2=0; k2 < 2; k2++) {
222 if ( (k == 1) || (k == np+1) ) {
224 for (l=0; l<nt; l++) {
225 for (i=0; i<nr; i++) {
235 for (l=0; l<m/2; l++) {
236 for (i=0; i<nr; i++) {
242 for (l=m/2; l<nt; l++) {
243 for (i=0; i<nr; i++) {
247 jmin = ( m == 0 ) ? l : 0 ;
248 for (j=jmin; j<nt; j++) {
249 double amlj = aa[nt*l + j] ;
250 for (i=0; i<nr; i++) {
251 som[i] += amlj * cc[nr*j + i] ;
255 for (i=0; i<nr; i++) {
277 assert(resu == cfo + (np+2)*ntnr) ;