Contents
     zfftf - compute  the  Fourier  coefficients  of  a  periodic
     sequence.  The FFT operations are unnormalized, so a call of
     ZFFTF followed by a call of ZFFTB will  multiply  the  input
     sequence by N.
     SUBROUTINE ZFFTF(N, X, WSAVE)
     DOUBLE COMPLEX X(*)
     INTEGER N
     DOUBLE PRECISION WSAVE(*)
     SUBROUTINE ZFFTF_64(N, X, WSAVE)
     DOUBLE COMPLEX X(*)
     INTEGER*8 N
     DOUBLE PRECISION WSAVE(*)
  F95 INTERFACE
     SUBROUTINE FFTF([N], X, WSAVE)
     COMPLEX(8), DIMENSION(:) :: X
     INTEGER :: N
     REAL(8), DIMENSION(:) :: WSAVE
     SUBROUTINE FFTF_64([N], X, WSAVE)
     COMPLEX(8), DIMENSION(:) :: X
     INTEGER(8) :: N
     REAL(8), DIMENSION(:) :: WSAVE
  C INTERFACE
     #include <sunperf.h>
     void zfftf(int n, doublecomplex *x, double *wsave);
     void zfftf_64(long n, doublecomplex *x, double *wsave);
     N (input) Length of the sequence to be  transformed.   These
               subroutines are most efficient when N is a product
               of small primes.  N >= 0.
     X (input) On entry, an array  of  length  N  containing  the
               sequence to be transformed.
     WSAVE (input)
               On entry, WSAVE must be an array of dimension (4 *
               N  + 15) or greater and must have been initialized
               by ZFFTI.