33 void nfft_benchomp_createdataset(
unsigned int d, 
unsigned int trafo_adjoint, 
int *N, 
int M, 
double sigma)
 
   41   for (t = 0; t < d; t++)
 
   49     n[t] = sigma*nfft_next_power_of_2(N[t]);
 
   63   printf(
"%d %d ", d, trafo_adjoint);
 
   76       printf(
"%.16e ", x[d*j+t]);
 
   82     for (j=0; j < N_total; j++)
 
   83       printf(
"%.16e %.16e\n", creal(f_hat[j]), cimag(f_hat[j]));
 
   88       printf(
"%.16e %.16e\n", creal(f[j]), cimag(f[j]));
 
   96 int main(
int argc, 
char **argv)
 
  106     fprintf(stderr, 
"usage: d tr_adj N_1 ... N_d M sigma\n");
 
  112   fprintf(stderr, 
"d=%d", d);
 
  114   if (d < 1 || argc < 5+d) {
 
  115     fprintf(stderr, 
"usage: d tr_adj N_1 ... N_d M sigma\n");
 
  119   N = malloc(d*
sizeof(
int));
 
  121   trafo_adjoint = atoi(argv[2]);
 
  122   if (trafo_adjoint < 0 && trafo_adjoint > 1)
 
  125   fprintf(stderr, 
", tr_adj=%d, N=", trafo_adjoint);
 
  128     N[t] = atoi(argv[3+t]);
 
  131     fprintf(stderr, 
"%d ",N[t]);
 
  135   sigma = atof(argv[4+d]);
 
  137   fprintf(stderr, 
", M=%d, sigma=%.16g\n", M, sigma);
 
  139   nfft_benchomp_createdataset(d, trafo_adjoint, N, M, sigma);