Who konws how to use FFTW guru interface?

From: jamesyx (jamesxy7_at_gmail.com)
Date: 09/27/04 

Date: 27 Sep 2004 06:37:37 -0700

I work with it for a few days, but it does not work even for a simple
test, such as transfer use fftw_plan_guru_r2r for cos(x).

There always causes cores when I excute the plan, and it also warnning
me with "required for this conversiono of constant".

fftw_plan planf;
int rank = 1;
int howmany_rank = 1;
double *f;
fftw_iodim *dims, *howmany_dims;

dims = (fftw_iodim *)malloc(rank*sizeof(fftw_iodim));
howmany_dims = (fftw_iodim *)malloc(howmany_rank*sizeof(fftw_iodim));

for (i=1; i<=rank; i++){
  dims[i-1].n = M;
  dims[i-1].is = 0;
  dims[i-1].os = 0;}
for (i=1; i<=howmany_rank; i++){
  howmany_dims[i-1].n = M;
  howmany_dims[i-1].is = 1;
  howmany_dims[i-1].os = 1;}

planf = fftw_plan_guru_r2r(rank, dims, howmany_rank, howmany_dims, f,
f,
                           FFTW_REDFT01, FFTW_EXHAUSTIVE);

......
fftw_execute_r2r(planf, f, f);
......