From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- sys/vops/aiftxr.x | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100644 sys/vops/aiftxr.x (limited to 'sys/vops/aiftxr.x') diff --git a/sys/vops/aiftxr.x b/sys/vops/aiftxr.x new file mode 100644 index 00000000..a9647e7c --- /dev/null +++ b/sys/vops/aiftxr.x @@ -0,0 +1,27 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# AIFTXR -- Inverse fourier transform (complex transform, real arrays). +# The inverse transform of the real arrays FR and FI containing complex data +# pairs is computed and the complex spatial data coefficients are returned in +# the real arrays SR and SI. The transformation may be performed in place if +# desired. NPIX must be a power of 2. + +procedure aiftxr (fr, fi, sr, si, npix) + +real fr[ARB], fi[ARB] # transform, frequency domain (input) +real sr[ARB], si[ARB] # data, spatial domain (output) +int npix +int ier + +begin + # The following are no-ops if the transform is being performed + # in place. + + call amovr (fr, sr, npix) + call amovr (fi, si, npix) + + # Compute the inverse transform. + call fft842 (1, npix, sr, si, ier) + if (ier == 1) + call fatal (1, "afftxr: npix not a power of 2") +end -- cgit