26 #define NFFT_PRECISION_DOUBLE
30 static void simple_test_nfct_1d(
void)
34 const
char *error_str;
40 NFCT(init_1d)(&p,N,M);
43 NFFT(vrand_real)(p.x, p.M_total, NFFT_K(0.0), NFFT_K(0.5));
46 if( p.flags & PRE_ONE_PSI)
47 NFCT(precompute_one_psi)(&p);
50 NFFT(vrand_real)(p.f_hat, p.N_total, NFFT_K(0.0), NFFT_K(1.0));
51 NFFT(vpr_double)(p.f_hat,p.N_total,"given Fourier coefficients, vector f_hat");
54 error_str = NFCT(check)(&p);
57 printf(
"Error in nfct module: %s\n", error_str);
62 NFCT(trafo_direct)(&p);
63 NFFT(vpr_double)(p.f,p.M_total,
"ndct, vector f");
67 NFFT(vpr_double)(p.f,p.M_total,
"nfct, vector f");
70 NFCT(adjoint_direct)(&p);
71 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint ndct, vector f_hat");
75 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint nfct, vector f_hat");
83 printf(
"Computing one dimensional ndct, nfct, adjoint ndct, and adjoint nfct...\n\n");
84 simple_test_nfct_1d();