24 #define NFFT_PRECISION_DOUBLE
28 static void simple_test_nfst_1d(
void)
32 const
char *error_str;
38 NFST(init_1d)(&p,N,M);
41 NFFT(vrand_real)(p.x, p.M_total, NFFT_K(0.0), NFFT_K(0.5));
45 NFST(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 = NFST(check)(&p);
55 printf(
"Error in nfst module: %s\n", error_str);
60 NFST(trafo_direct)(&p);
61 NFFT(vpr_double)(p.f,p.M_total,
"ndst, vector f");
65 NFFT(vpr_double)(p.f,p.M_total,
"nfst, vector f");
68 NFST(adjoint_direct)(&p);
69 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint ndst, vector f_hat");
73 NFFT(vpr_double)(p.f_hat,p.N_total,
"adjoint nfst, vector f_hat");
81 printf(
"Computing one dimensional ndct, nfct, adjoint ndct, and adjoint nfct...\n\n");
82 simple_test_nfst_1d();