26 #define NFFT_PRECISION_DOUBLE
30 static void simple_test_nfst_1d(
void)
34 const
char *error_str;
40 NFST(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 NFST(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 = NFST(check)(&p);
57 printf(
"Error in nfst module: %s\n", error_str);
62 NFST(trafo_direct)(&p);
63 NFFT(vpr_double)(p.f,p.M_total,
"ndst, vector f");
67 NFFT(vpr_double)(p.f,p.M_total,
"nfst, vector f");
70 NFST(adjoint_direct)(&p);
71 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint ndst, vector f_hat");
75 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint nfst, vector f_hat");
83 printf(
"Computing one dimensional ndct, nfct, adjoint ndct, and adjoint nfct...\n\n");
84 simple_test_nfst_1d();