23char cirjaco02_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/cirjaco02.C,v 1.4 2014/10/13 08:53:11 j_novak Exp $" ;
117void cirjaco02(
const int* deg,
const int* dimc,
double* cf,
const int* dimf,
136 cout <<
"circheb: nr > n3c : nr = " << nr <<
" , n3c = "
142 cout <<
"circheb: nr > n3f : nr = " << nr <<
" , n3f = "
148 cout <<
"circheb: n1c > n1f : n1c = " << n1c <<
" , n1f = "
154 cout <<
"circheb: n2c > n2f : n2c = " << n2c <<
" , n2f = "
165 int n2n3f = n2f * n3f ;
166 int n2n3c = n2c * n3c ;
174 int borne_phi = ( n1c > 1 ) ? n1c-1 : 1 ;
176 for (j=0; j< borne_phi; j++) {
180 for (k=0; k<n2c; k++) {
182 int i0 = n2n3c * j + n3c * k ;
183 double* cf0 = cf + i0 ;
185 i0 = n2n3f * j + n3f * k ;
186 double* ff0 = ff + i0 ;
188 Tbl jj = jacobipointsgl(nm1) ;
190 for (i=0 ; i<nr ; i++) {
192 for (
int n = 0 ; n<nr ; n++){
193 som += cf0[n]*jj(n,i) ;