void cffti (long n,double wsave[])
{
/* implicit double (a-h,o-z);*/
#undef wsave_1
#define wsave_1(a1) wsave[a1-1]
/*-----implicit-declarations-----*/
long iw2;
long iw1;
/*-----end-of-declarations-----*/
if (n == 1) return;
iw1 = n+n+1;
iw2 = iw1+n+n;
cffti1 (n,&wsave_1(iw1),(long *)&wsave_1(iw2));
return;
}
int cffti1(int *n, double *wa, double *ifac) {
// Casting (double*) to (int*) is dangerous if sizes differ!!!
return cffti1(n, wa, (int *)ifac);
}
