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)();