24 #define NFFT_PRECISION_DOUBLE
28 static void simple_test_nfct_1d(
void)
32 const
char *error_str;
38 NFCT(init_1d)(&p,N,M);
41 NFFT(vrand_real)(p.x, p.M_total, NFFT_K(0.0), NFFT_K(0.5));
45 NFCT(precompute_one_psi)(&p);
48 NFFT(vrand_real)(p.f_hat, p.N_total, NFFT_K(0.0), NFFT_K(1.0));
49 NFFT(vpr_double)(p.f_hat,p.N_total,"given Fourier coefficients, vector f_hat");
52 error_str = NFCT(check)(&p);
55 printf(
"Error in nfct module: %s\n", error_str);
60 NFCT(trafo_direct)(&p);
61 NFFT(vpr_double)(p.f,p.M_total,
"ndct, vector f");
65 NFFT(vpr_double)(p.f,p.M_total,
"nfct, vector f");
68 NFCT(adjoint_direct)(&p);
69 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint ndct, vector f_hat");
73 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint nfct, vector f_hat");
81 printf(
"Computing one dimensional ndct, nfct, adjoint ndct, and adjoint nfct...\n\n");
82 simple_test_nfct_1d();