Changeset 1401 in Sophya for trunk/SophyaExt/IFFTW

Timestamp:

Feb 13, 2001, 3:30:04 PM (25 years ago)

Author:

ansari

Message:

Modifs-debug FFTWServer (suite) - Reza 13/2/2001

File:

• trunk/SophyaExt/IFFTW/fftwserver.cc (modified) (3 diffs)

trunk/SophyaExt/IFFTW/fftwserver.cc

@@ -209 +209 @@
     if (_pndb) _pndb->Recreate(in.NbDimensions(), sz);
     else _pndb = new FFTWServerPlan(in.NbDimensions(), sz, FFTW_BACKWARD, false);
+    fftwnd_one(_pndb->pnd, (fftw_complex *)(in.Data()) , (fftw_complex *)(out.Data()) );
   }
 
@@ -218 +219 @@
 {  
   int rank = ckR8.CheckResize(in, out);
-  TArray<r_8> outtemp(in, false); 
-
   if (rank == 1) { // One dimensional transform 
     if (_p1drf) _p1drf->Recreate(in.Size());
     else _p1drf = new FFTWServerPlan(in.Size(), FFTW_REAL_TO_COMPLEX, true);
+    TArray<r_8> outtemp; 
+    outtemp.ReSize(in);
     rfftw_one(_p1drf->rp, (fftw_real *)(in.Data()) , (fftw_real *)(outtemp.Data()));
     ReShapetoCompl(outtemp, out);
@@ -252 +253 @@
   int rank = ckR8.CheckResize(in, out);
   if (rank == 1) { // One dimensional transform 
-    TArray<r_8> intemp(out, false); 
+    TArray<r_8> intemp;
+    intemp.ReSize(out);
     if (_p1drb) _p1drb->Recreate(out.Size());
     else _p1drb = new FFTWServerPlan(out.Size(), FFTW_COMPLEX_TO_REAL, true);
 
     ReShapetoReal(in, intemp);
-    cerr << " DEBUG-FFTWServer::FFTBackward() in = \n" << in << endl;
-    cerr << " DEBUG-FFTWServer::FFTBackward() intemp = \n" << intemp << endl;
+    //    cerr << " DEBUG-FFTWServer::FFTBackward() in = \n" << in << endl;
+    //    cerr << " DEBUG-FFTWServer::FFTBackward() intemp = \n" << intemp << endl;
     rfftw_one(_p1drb->rp, (fftw_real *)(intemp.Data()) , (fftw_real *)(out.Data()));