29 static void simple_test_nsfft(
int d,
int J,
int M)
34 nsfft_init(&p, d, J, M, 6, NSDFT);
36 nsfft_init_random_nodes_coeffs(&p);
41 nsfft_trafo_direct(&p);
49 nsfft_adjoint_direct(&p);
60 int main(
int argc,
char **argv)
65 printf(
"1) computing a two dimensional nsdft, nsfft and adjoints\n\n");
68 M=(J+4)*nfft_exp2i(J+1);
69 simple_test_nsfft(d,J,M);
73 printf(
"2) computing a three dimensional nsdft, nsfft and adjoints\n\n");
76 M=6*nfft_exp2i(J)*(nfft_exp2i((J+1)/2+1)-1)+nfft_exp2i(3*(J/2+1));
77 simple_test_nsfft(d,J,M);
void nfft_vpr_complex(fftw_complex *x, const NFFT_INT n, const char *text)
Print complex vector to standard output.
data structure for an NSFFT (nonequispaced sparse fast Fourier transform) plan with double precision ...
fftw_complex * f_hat
Fourier coefficients.