Changeset 16107 for trunk/Ohana/src/opihi/cmd.data/fft1d.c

Timestamp:

Jan 17, 2008, 8:41:08 AM (19 years ago)

Author:

eugene

Message:

converting fft and fftn to code from Jones, with Magnier version of multi-dim

File:

• trunk/Ohana/src/opihi/cmd.data/fft1d.c (modified) (3 diffs)

trunk/Ohana/src/opihi/cmd.data/fft1d.c

@@ -3 +3 @@
 int fft1d (int argc, char **argv) {
   
-  int Npix, Nbit, ZeroImaginary;
+  int N, Npix, Nbit, ZeroImaginary, forward;
   Vector *Ire, *Iim, *Ore, *Oim;
+
+  forward = TRUE;
+  if ((N = get_argument (argc, argv, "-inverse"))) {
+    remove_argument (N, &argc, argv);
+    forward = FALSE;
+  }
 
   if ((argc != 6) || (strcmp (argv[3], "to"))) {
@@ -22 +28 @@
   Npix = Ire[0].Nelements;
   if (!ZeroImaginary && (Npix != Iim[0].Nelements)) {
-    gprint (GP_ERR, "vector mismatch in size\n");
+    gprint (GP_ERR, "vector size mismatch\n");
     return (FALSE);
   }
@@ -50 +56 @@
   }    
 
-  fft (Ore[0].elements, Oim[0].elements, Npix, Nbit, TRUE); 
+  fft1D (Ore[0].elements, Oim[0].elements, Npix, Nbit, forward); 
   
   return (TRUE);