zfftzm - initialize the trigonometric weight and factor tables or compute the one-dimensional Fast Fourier Transform (forward or inverse) of a set of data sequences stored in a two-dimensional double complex array. =head1 SYNOPSIS
SUBROUTINE ZFFTZM( IOPT, M, N, SCALE, X, LDX, Y, LDY, TRIGS, IFAC, * WORK, LWORK, IERR) DOUBLE COMPLEX X(LDX,*), Y(LDY,*) INTEGER IOPT, M, N, LDX, LDY, LWORK, IERR INTEGER IFAC(*) DOUBLE PRECISION SCALE DOUBLE PRECISION TRIGS(*), WORK(*)
SUBROUTINE ZFFTZM_64( IOPT, M, N, SCALE, X, LDX, Y, LDY, TRIGS, * IFAC, WORK, LWORK, IERR) DOUBLE COMPLEX X(LDX,*), Y(LDY,*) INTEGER*8 IOPT, M, N, LDX, LDY, LWORK, IERR INTEGER*8 IFAC(*) DOUBLE PRECISION SCALE DOUBLE PRECISION TRIGS(*), WORK(*)
SUBROUTINE FFTM( IOPT, [M], [N], [SCALE], X, [LDX], Y, [LDY], TRIGS, * IFAC, WORK, [LWORK], IERR) COMPLEX(8), DIMENSION(:,:) :: X, Y INTEGER :: IOPT, M, N, LDX, LDY, LWORK, IERR INTEGER, DIMENSION(:) :: IFAC REAL(8) :: SCALE REAL(8), DIMENSION(:) :: TRIGS, WORK
SUBROUTINE FFTM_64( IOPT, [M], [N], [SCALE], X, [LDX], Y, [LDY], * TRIGS, IFAC, WORK, [LWORK], IERR) COMPLEX(8), DIMENSION(:,:) :: X, Y INTEGER(8) :: IOPT, M, N, LDX, LDY, LWORK, IERR INTEGER(8), DIMENSION(:) :: IFAC REAL(8) :: SCALE REAL(8), DIMENSION(:) :: TRIGS, WORK
#include <sunperf.h>
void zfftzm(int iopt, int m, int n, double scale, doublecomplex *x, int ldx, doublecomplex *y, int ldy, double *trigs, int *ifac, double *work, int lwork, int *ierr);
void zfftzm_64(long iopt, long m, long n, double scale, doublecomplex *x, long ldx, doublecomplex *y, long ldy, double *trigs, long *ifac, double *work, long lwork, long *ierr);
zfftzm initializes the trigonometric weight and factor tables or computes the one-dimensional Fast Fourier Transform (forward or inverse) of a set of data sequences stored in a two-dimensional double complex array: .Ve
N1-1
Y(k,l)
= SUM W*X(j,l)
j=0 .Ve
where
k ranges from 0 to N1-1 and l ranges from 0 to N2-1
i = sqrt(-1)
isign = 1 for inverse transform or -1 for forward transform
W = exp(isign*i*j*k*2*pi/N1)
.Ve
IOPT = 0 computes the trigonometric weight table and factor table
IOPT = -1 computes forward FFT
IOPT = +1 computes inverse FFT
0 = normal return
-1 = IOPT is not 0, 1 or -1
-2 = M < 0
-3 = N < 0
-4 = (LDX < M)
-5 = (LDY < M) or (LDY not equal LDX when X and Y are same array)
-6 = (LWORK not equal 0) and (LWORK < 2*M*NCPUS)
-7 = memory allocation failed
fft