26 #define NFFT_PRECISION_DOUBLE
33 const int N = 1000000;
34 const int M = 1000000;
37 printf(
"nthreads = " NFFT__D__
"\n", NFFT(get_num_threads)());
41 NFFT(init_1d)(&p,N,M);
44 NFFT(vrand_shifted_unit_double)(p.x,p.M_total);
47 t0 = NFFT(clock_gettime_seconds)();
49 NFFT(precompute_one_psi)(&p);
50 t1 = NFFT(clock_gettime_seconds)();
51 fprintf(stderr,
"precompute elapsed time: %.3" NFFT__FIS__
" seconds\n",t1-t0);
54 NFFT(vrand_unit_complex)(p.f_hat,p.N_total);
57 t0 = NFFT(clock_gettime_seconds)();
59 t1 = NFFT(clock_gettime_seconds)();
60 fprintf(stderr,
"compute elapsed time: %.3" NFFT__FIS__
" seconds\n",t1-t0);
65 FFTW(cleanup_threads)();