XLPack 7.0
XLPack Numerical Library (C API) Reference Manual
Loading...
Searching...
No Matches

◆ cfftmf()

void cfftmf ( int  lot,
int  jump,
int  n,
int  inc,
doublecomplex  c[],
int  lc,
double  wsave[],
int  lwsave,
double  work[],
int  lwork,
int *  info 
)

One-dimensional complex Fourier transform (multiple sequences)

Purpose
cfftmf computes the one-dimensional Fourier transform of multiple periodic sequences within a real array. This is referred to as the forward transform or Fourier analysis, transforming the sequence from physical to spectral space.
c[l*jump+k] = (1/n)Σc[l*jump+j]exp(-2πijk/n) (Σ for j = 0 to n-1) (l = 0 to lot-1, k = 0 to n-1) (i is imaginary unit)
This transform is normalized since a call to cfftmf followed by a call to cfftmb (or vice-versa) reproduces the original array subject to algorithmic constraints, roundoff error, etc.
Parameters
[in]lotThe number of sequences to be transformed. (lot >= 1)
[in]jumpIncrement between the locations, in array c[], of the first elements of two consecutive sequences to be transformed. (jump >= 1)
[in]nThe length of the sequence to be transformed. (n >= 1)
The transform is most efficient when n is a product of small primes.
[in]incInteger increment between the locations, in array c[], of two consecutive elements within the sequence. (inc >= 1)
[in,out]c[]Array c[lenc]
[in] The sequence to be transformed.
[out] The Fourier forward transformed sequences of data.
[in]lcThe length of the array c[]. (lc >= inc*(n - 1) + 1)
[in]wsave[]Array wsave[lwsave]
Work data. Its contents must be initialized with a call to cfftmi before the first call to cfftmf or cfftmb for a given transform length n.
[in]lwsaveThe length of the array wsave[]. (lwsave >= 2*n + ln(n)/ln(2) + 4)
[out]work[]Array work[lwork]
Work array.
[in]lworkThe length of the array work[]. (lwork >= 2*lot*n)
[out]info= 0: Successful exit
= -1: The argument lot had an illegal value (lot < 1)
= -2: The argument jump had an illegal value (jump < 1)
= -3: The argument n had an illegal value (n < 1)
= -4: The argument inc had an illegal value (inc < 1)
= -6: The argument lc had an illegal value (lc not big enough)
= -8: The argument lwsave had an illegal value (lwsave not big enough)
= -10: The argument lwork had an illegal value (lwork not big enough)
Reference
FFTPACK 5.1