◆ Cfftmb()

One-dimensional complex Fourier backward transform (multiple sequences)

Purpose: This routine computes the one-dimensional Fourier transform of periodic sequences within a real array. This is referred to as the backward transform or Fourier synthesis, transforming the sequence from spectral to physical 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 Cfftmb followed by a call to Cfftmf (or vice-versa) reproduces the original array subject to algorithmic constraints, roundoff error, etc.

Parameters

[in]	Lot	Number of sequences to be transformed. (Lot >= 1)
[in]	Jump	Increment between the locations, in array C(), of the first elements of two consecutive sequences to be transformed. (Jump >= 1)
[in]	N	Length of the sequence to be transformed. (N >= 1) (The transform is most efficient when N is a product of small primes)
[in,out]	C()	Array C(LC - 1) (LC >= (Lot - 1)Jump + IncC(N - 1) + 1) [in] The sequences to be transformed. [out] The Fourier backward transformed sequences of data.
[in]	Wsave()	Array Wsave(LWsave - 1) (LWsave >= 2*N + ln(N)/ln(2) + 4) 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.
[out]	Info	= 0: Successful exit. = -1: The argument Lot had an illegal value. (Lot < 1, or, Lot, Jump, N and Inc are inconsistent) = -2: The argument Jump had an illegal value. (Jump < 1) = -3: The argument N had an illegal value. (N < 1) = -4: The argument C() had an illegal value. (Array C() is not big enough) = -5: The argument Wsave() had an illegal value. (Array Wsave() is not big enough) = -7: The argument Inc had an illegal value. (Inc < 1)
[in]	Inc	(Optional) Integer increment between the locations, in array C(), of two consecutive elements within the sequence. (Inc >= 1) (default = 1)