118 fftw_plan prepare_fft(
int, Tbl*&) ;
121 void cfpcossin(
const int* deg,
const int* dim,
double* cf)
133 cout <<
"cfpcossin: np+2 > n1 : np+2 = " << np+2 <<
" , n1 = " 141 fftw_plan p = prepare_fft(np, pg) ;
146 double fac = 2./double(np) ;
149 for (
int j=0; j<n2; j++) {
150 for (
int k=0; k<n3; k++) {
153 double* debut = cf + index ;
155 for (
int i=0; i<np; i++) {
165 double* psin = g.t + np - 1 ;
166 (*debut) = (*pcos)/double(np) ;
167 debut += n2n3 ; pcos++ ;
170 for (
int i=1; i<np/2; i++){
171 *debut = (*pcos)*fac ;
172 debut += n2n3 ; pcos++ ;
173 *debut = -(*psin)*fac ;
174 debut += n2n3 ; psin-- ;
176 (*debut) = (*pcos)/double(np) ;