23char circheb_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/FFTW3/circheb.C,v 1.3 2014/10/13 08:53:19 j_novak Exp $" ;
137fftw_plan back_fft(
int, Tbl*&) ;
138double* cheb_ini(
const int) ;
141void circheb(
const int* deg,
const int* dimc,
double* cf,
const int* dimf,
160 cout <<
"circheb: nr > n3c : nr = " << nr <<
" , n3c = "
166 cout <<
"circheb: nr > n3f : nr = " << nr <<
" , n3f = "
172 cout <<
"circheb: n1c > n1f : n1c = " << n1c <<
" , n1f = "
178 cout <<
"circheb: n2c > n2f : n2c = " << n2c <<
" , n2f = "
190 fftw_plan p = back_fft(nm1, pg) ;
194 double* sinp = cheb_ini(nr);
198 int n2n3f = n2f * n3f ;
199 int n2n3c = n2c * n3c ;
207 int borne_phi = ( n1c > 1 ) ? n1c-1 : 1 ;
209 for (j=0; j< borne_phi; j++) {
213 for (k=0; k<n2c; k++) {
215 int i0 = n2n3c * j + n3c * k ;
216 double* cf0 = cf + i0 ;
218 i0 = n2n3f * j + n3f * k ;
219 double* ff0 = ff + i0 ;
237 for ( i = 3; i < nr; i += 2 ) {
238 ff0[i] = cf0[i] - c1 ;
242 double fmoins0 = - nm1s2 * c1 - som ;
247 for ( i = 3; i < nr; i += 2 ) {
248 g.set(nm1-i/2) = -0.25 * ( ff0[i] - ff0[i-2] ) ;
259 for (i=1; i<nm1s2; i ++ ) g.set(i) = 0.5 * cf0[2*i] ;
260 g.set(nm1s2) = cf0[nm1] ;
271 for ( i = 1; i < nm1s2 ; i++ ) {
275 double fp = .5 * ( g(i) + g(isym) ) ;
276 double fm = .5 * ( g(i) - g(isym) ) / sinp[i] ;
279 ff0[isym] = fp - fm ;
283 ff0[0] = g(0) + fmoins0 ;
284 ff0[nm1] = g(0) - fmoins0 ;
285 ff0[nm1s2] = g(nm1s2) ;