Changeset 538 in Sophya

Timestamp:

Nov 3, 1999, 3:30:40 PM (26 years ago)

Author:

ansari

Message:

From Alex (compte pixel) A.Kim/Reza 3/11/99

File:

• trunk/SophyaLib/NTools/fftserver.cc (modified) (4 diffs)

trunk/SophyaLib/NTools/fftserver.cc

@@ -3 +3 @@
 
 extern "C" {
-  int rffti_(int *n, float *wsave);
+  int rffti_(int *n, float *wsave);             //fftpack functions
   int rfftf_(int *n, float *r__, float *wsave);
   int rfftb_(int *n, float *r__, float *wsave);
@@ -19 +19 @@
 FFTServer::FFTServer()
 {
-sz_rfft = 0;
-ws_rfft = NULL;
-sz_cfft = 0;
-ws_cfft = NULL;
-sz_cdfft = 0;
-ws_cdfft = NULL;
+  sz_rfft = 0;     //the working array and its size for the different
+  ws_rfft = NULL;  //possible numerical types
+  sz_cfft = 0;
+  ws_cfft = NULL;
+  sz_cdfft = 0;
+  ws_cdfft = NULL;
 }
 
@@ -36 +36 @@
 void FFTServer::checkint_rfft(int l)
 {
-  if (sz_rfft == l) return;
-
-  if (ws_rfft) delete[] ws_rfft;
+  if (sz_rfft == l) return;       //checkint functions check and reallocate
+                                  //memory for the work arrays when performing
+  if (ws_rfft) delete[] ws_rfft;  //a transform
   sz_rfft = l;
   ws_rfft = new float[2*l+15];
@@ -73 +73 @@
   cdffti_(&l, ws_cdfft);
 }
+
+/* In general forward transformations are resorted since fftpack functions
+   return inverse transformations */
 
 void FFTServer::fftf(int l, float* inout)