Changeset 467 in Sophya for trunk/SophyaLib/NTools/fftserver.cc

Timestamp:

Oct 15, 1999, 2:05:15 PM (26 years ago)

Author:

ansari

Message:

Fixing of bugs, output formats, inverse transforms.

File:

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

trunk/SophyaLib/NTools/fftserver.cc

@@ -23 +23 @@
 sz_cfft = 0;
 ws_cfft = NULL;
+sz_cdfft = 0;
+ws_cdfft = NULL;
 }
 
@@ -29 +31 @@
 if (ws_rfft) delete[] ws_rfft;
 if (ws_cfft) delete[] ws_cfft;
+if (ws_cdfft) delete[] ws_cdfft;
 }
 
@@ -75 +78 @@
   checkint_rfft(l);
   rfftf_(&l, inout, ws_rfft);
+  for (int k= 2;k<=(l+1)/2;k++) inout[2*k-2]=-inout[2*k-2];
 }
 
@@ -81 +85 @@
   checkint_dfft(l);
   dfftf_(&l, inout, ws_dfft);
+  for (int k= 2;k<=(l+1)/2;k++) inout[2*k-2]=-inout[2*k-2];
 }
 
@@ -92 +97 @@
   }
   cfftf_(&l, foo, ws_cfft);
-  for (int i=0;i<l;i++) inout[i]= complex<float> (foo[2*i], foo[2*i+1]);
+  inout[0]=complex<float> (foo[0],foo[1]);
+  for (int i=1;i<l;i++) inout[l-i]= complex<float> (foo[2*i], foo[2*i+1]);
   delete[] foo;
 }
@@ -99 +105 @@
 {
   checkint_cdfft(l);
-  double* foo = new double[2*l];
+  double* foo=new double[2*l];
   for (int i=0;i<l;i++){
     foo[2*i]=inout[i].real();
@@ -105 +111 @@
   }
   cdfftf_(&l, foo, ws_cdfft);
-  for (int i=0;i<l;i++) inout[i]= complex<double> (foo[2*i],foo[2*i+1]);
+  inout[0]=complex<double> (foo[0],foo[1]);
+  for (int i=1;i<l;i++) {
+    inout[l-i]= complex<double> (foo[2*i],foo[2*i+1]);
+  }
   delete[] foo;
 }
@@ -112 +121 @@
 {
   checkint_rfft(l);
-  rfftb_(&l, inout, ws_rfft);
+  rfftf_(&l, inout, ws_rfft);
 }
 
@@ -118 +127 @@
 {
   checkint_dfft(l);
-  dfftb_(&l, inout, ws_dfft);
+  dfftf_(&l, inout, ws_dfft);
 }
 
@@ -129 +138 @@
     foo[2*i+1]=inout[i].imag();
   }
-  cfftb_(&l, foo, ws_cfft);
+  cfftf_(&l, foo, ws_cfft);
   for (int i=0;i<l;i++) inout[i]=complex<float> (foo[2*i],foo[2*i+1]);
   delete[] foo;
@@ -142 +151 @@
     foo[2*i+1]=inout[i].imag();
   }
-  cdfftb_(&l, foo, ws_cdfft);
+  cdfftf_(&l, foo, ws_cdfft);
   for (int i=0;i<l;i++) inout[i]=complex<double> (foo[2*i],foo[2*i+1]);
   delete[] foo;