Now that I'm getting this filtering stuff implemented, I want to have
options to export frequency-domain data to files so that I can look at
them easily with kappa (for example). The FFT of the data, and the
filters, are complex valued. FFTW has the option to store the real and
imaginary parts in separate real-valued arrays. However, the default
(which I'm currently using) is to define arrays of type
typedef double fftw_complex[2];
It looks like much of starlink knows what to do with complex numbers
because they are a native fortran type. Does the C interface support some
implementation of complex numbers? (a cursory glance at the C interface
documentation for NDF leads me to suspect that it doesn't).
I'm also going to want to change the units of the output WCS time axis to
frequency units (Hz)...
Ed
+----------------------------------------------------------------+
| Dr. Edward Chapin (604) 822-6709 (work) |
| Dept. of Physics and Astronomy |
| University of British Columbia |
| 6224 Agricultural Road, Vancouver, B.C. |
| V6T 1Z1, Canada [log in to unmask] |
+----------------------------------------------------------------+
|