23char cipcossin_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Coef/FFTW3/cipcossin.C,v 1.3 2014/10/13 08:53:19 j_novak Exp $" ;
129fftw_plan back_fft(
int, Tbl*&) ;
132void cipcossin(
const int* deg,
const int* dimc,
const int* dimf,
133 double* cf,
double* ff)
149 cout <<
"cipcossin: np+2 > n1c : np = " << np <<
" , n1c = "
155 cout <<
"cipcossin: np > n1f : np = " << np <<
" , n1f = "
161 cout <<
"cipcossin: n3f > n3c : n3f = " << n3f <<
" , n3c = "
167 cout <<
"cipcossin: n2f > n2c : n2f = " << n2f <<
" , n2c = "
175 fftw_plan p = back_fft(np, pg) ;
178 int n2n3c = n2c*n3c ;
179 int n2n3f = n2f*n3f ;
183 for (
int j=0; j<n2c; j++) {
184 for (
int k=0; k<n3c; k++) {
185 index = n3c * j + k ;
186 double* debut = cf + index ;
189 for (
int i=1; i<np/2; i++) {
191 g.set(i) = 0.5 * (*debut) ; debut += n2n3c ;
192 g.set(isym) = -0.5 * (*debut) ; debut += n2n3c ;
194 g.set(np/2) = *debut ;
201 if ((j<n2f) && (k<n3f)) {
202 debut = ff + n3f * j + k ;
203 for (
int i=0; i<np; i++) {