1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
# AFFTXX -- Forward fourier transform (complex transform, complex data).
# The fourier transform of the complex array A of length NPIX pixels is
# computed and the NPIX complex transform coefficients are returned in the
# complex array B. The transformation may be performed in place if desired.
# NPIX must be a power of 2.
procedure afftxx (a, b, npix)
complex a[ARB] # data (input)
complex b[ARB] # transform (output)
int npix
int ier
pointer sp, xr, xi
begin
call smark (sp)
call salloc (xr, npix, TY_REAL)
call salloc (xi, npix, TY_REAL)
# Rearrange the elements of the A array as required by FFT842.
# Convert the array A of complex values into an array of reals
# and an array of imaginaries.
call aupxr (a, Memr[xr], Memr[xi], npix)
# Compute the forward transform.
call fft842 (0, npix, Memr[xr], Memr[xi], ier)
if (ier == 1)
call fatal (1, "afftxx: npix not a power of 2")
# Repack the real and imaginary arrays to form the complex output
# array.
call apkxr (Memr[xr], Memr[xi], b, npix)
call sfree (sp)
end
|
