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