Changeset 1405 in Sophya for trunk/SophyaExt/IFFTW/fftwserver.cc

Timestamp:

Feb 15, 2001, 6:58:16 PM (25 years ago)

Author:

ansari

Message:

Ajout documentation - Reza 15/2/2001

File:

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

trunk/SophyaExt/IFFTW/fftwserver.cc

@@ -3 +3 @@
 #include "FFTW/rfftw.h"
 
+/*! 
+  \class SOPHYA::FFTWServer
+  \ingroup IFFTW
+  An implementation of FFTServerInterface based on FFTW, double
+  precision arrays, using FFTW package, availabale from http://www.fftw.org.
+  Refer to FFTServerInterface for details about FFTServer operations.
+
+  \code
+  #include "fftwserver.h"
+  // ...
+  TMatrix<r_8> in(24,32);
+  TMatrix< complex<r_8> > out;
+  in = RandomSequence();
+  FFTWServer ffts;
+  ffts.setNormalize(true);  // To have normalized transforms
+  cout << " FFTServer info string= " << ffts.getInfo() << endl;
+  cout << "in= " << in << endl;
+  cout << " Calling ffts.FFTForward(in, out) : " << endl;
+  ffts.FFTForward(in, out);
+  cout << "out= " << out << endl;
+  \endcode
+
+*/
 
 #define MAXND_FFTW 5
+
+namespace SOPHYA {
 
 class FFTWServerPlan {
@@ -27 +52 @@
    
 };
+
+} // Fin du namespace
 
 FFTWServerPlan::FFTWServerPlan(int n, fftw_direction dir, bool fgreal)
@@ -40 +67 @@
   _n = n;  
   _dir = dir;
-  if (fgreal) rp = rfftw_create_plan(n, dir, FFTW_ESTIMATE);
-  else p = fftw_create_plan(n, dir, FFTW_ESTIMATE);
+  if (fgreal) { 
+    rp = rfftw_create_plan(n, dir, FFTW_ESTIMATE);
+    if (rp == NULL) 
+      throw AllocationError("FFTWServerPlan: failed to create plan (rp) !");
+  }
+  else {
+    p = fftw_create_plan(n, dir, FFTW_ESTIMATE);
+    if (p == NULL) 
+      throw AllocationError("FFTWServerPlan: failed to create plan (p) !");
+  }
+    
 }
 
@@ -62 +98 @@
   for(k=nd; k<MAXND_FFTW; k++) _nxyz[k] = -10;
   _dir = dir;
-  if (fgreal) rpnd = rfftwnd_create_plan(_nd, _nxyz, dir, FFTW_ESTIMATE);
-  else pnd = fftwnd_create_plan(_nd, _nxyz, dir, FFTW_ESTIMATE);
+  if (fgreal) {
+    rpnd = rfftwnd_create_plan(_nd, _nxyz, dir, FFTW_ESTIMATE);
+    if (rpnd == NULL) 
+      throw AllocationError("FFTWServerPlan: failed to create plan (rpnd) !");
+  }
+  else { 
+    pnd = fftwnd_create_plan(_nd, _nxyz, dir, FFTW_ESTIMATE);
+    if (pnd == NULL) 
+      throw AllocationError("FFTWServerPlan: failed to create plan (pnd) !");
+  }
+
 }
 
@@ -86 +131 @@
     fftw_destroy_plan(p);
     p = fftw_create_plan(n, _dir, FFTW_ESTIMATE);
+    if (p == NULL) 
+      throw AllocationError("FFTWServerPlan::Recreate failed to create plan (p) !");
   }
   else {
     rfftw_destroy_plan(rp);
     rp = rfftw_create_plan(n, _dir, FFTW_ESTIMATE);
-  }
+    if (rp == NULL) 
+      throw AllocationError("FFTWServerPlan::Recreate failed to create plan (rp) !");
+  }
+
 }
 
@@ -117 +167 @@
     fftwnd_destroy_plan(pnd);
     pnd = fftwnd_create_plan(_nd, _nxyz, _dir, FFTW_ESTIMATE);
+    if (pnd == NULL) 
+      throw AllocationError("FFTWServerPlan::Recreate failed to create plan (pnd) !");
   }
   else {
     rfftwnd_destroy_plan(rpnd);
     rpnd = rfftwnd_create_plan(_nd, _nxyz, _dir, FFTW_ESTIMATE);
+    if (rpnd == NULL) 
+      throw AllocationError("FFTWServerPlan::Recreate failed to create plan (rpnd) !");
   }