33 void fastsum_benchomp_createdataset(
unsigned int d, 
int L, 
int M)
 
   52       x[k*d+j] = 2.0 * r_max * (
double)rand()/(double)RAND_MAX - r_max;
 
   55       r2 += x[k*d+j] * x[k*d+j];
 
   57     if (r2 >= r_max * r_max)
 
   73       y[k*d+j] = 2.0 * r_max * (
double)rand()/(double)RAND_MAX - r_max;
 
   76       r2 += y[k*d+j] * y[k*d+j];
 
   78     if (r2 >= r_max * r_max)
 
   84   printf(
"%d %d %d\n", d, L, M);
 
   89       printf(
"%.16e ", x[d*j+t]);
 
   94     printf(
"%.16e %.16e\n", creal(alpha[j]), cimag(alpha[j]));
 
   99       printf(
"%.16e ", y[d*j+t]);
 
  108 int main(
int argc, 
char **argv)
 
  115     fprintf(stderr, 
"usage: d L M\n");
 
  123   fprintf(stderr, 
"d=%d, L=%d, M=%d\n", d, L, M);
 
  125   fastsum_benchomp_createdataset(d, L, M);