Changeset 1402 in Sophya for trunk/SophyaLib/NTools/fftservintf.cc

Timestamp:

Feb 13, 2001, 7:12:07 PM (25 years ago)

Author:

ansari

Message:

Modifs FFTServerInterface + FFTPackServer (Suite - Fin ?) - Reza 13/2/2001

File:

• trunk/SophyaLib/NTools/fftservintf.cc (modified) (5 diffs)

trunk/SophyaLib/NTools/fftservintf.cc

@@ -41 +41 @@
 
 /* --Methode-- */
-void FFTServerInterface::FFTBackward(TArray< complex<r_8> > const &, TArray< r_8 > &)
+void FFTServerInterface::FFTBackward(TArray< complex<r_8> > const &, TArray< r_8 > &, bool)
 {
   throw NotAvailableOperation("FFTServer::FFTBackward(TArray...) Unsupported operation !");
@@ -68 +68 @@
 
 /* --Methode-- */
-void FFTServerInterface::FFTBackward(TArray< complex<r_4> > const &, TArray< r_4 > &)
+void FFTServerInterface::FFTBackward(TArray< complex<r_4> > const &, TArray< r_4 > &, bool)
 {
   throw NotAvailableOperation("FFTServer::FFTBackward(TArray r_4 ... ) Unsupported operation !");
@@ -183 +183 @@
 /* --Methode-- */
 template <class T>
-int FFTArrayChecker<T>::CheckResize(TArray< complex<T> > const & in, TArray< T > & out)
+int FFTArrayChecker<T>::CheckResize(TArray< complex<T> > const & in, TArray< T > & out, 
+                                    bool usoutsz)
 {
   int k;
@@ -204 +205 @@
       throw(SzMismatchError(msg));
     }
-  sa_size_t sz[BASEARRAY_MAXNDIMS];
-  if (ndg1 > 1) {
-    sz[0] = 2*in.Size(0)-1;
-    for(k=1; k<in.NbDimensions(); k++) 
-      sz[k] = in.Size(k);
+  if (usoutsz) { // We have to use output array size 
+    bool fgerr = false;
+    if (ndg1 > 1) {
+      if (in.Size(0) != out.Size(0)/2+1) fgerr = true;
+    }      
+    else {
+      if (in.Size(in.MaxSizeKA()) != out.Size(in.MaxSizeKA())/2+1) fgerr = true;
+    }
+    if (fgerr) {
+        msg = _msg + "CheckResize(complex in, real out) - Incompatible in-out sizes !";
+        throw(SzMismatchError(msg));
+    }
+  }
+  else {  // We have to resize the output array 
+    sa_size_t sz[BASEARRAY_MAXNDIMS];
+    if (ndg1 > 1) {
+      sz[0] = 2*in.Size(0)-1;
+      for(k=1; k<in.NbDimensions(); k++) 
+        sz[k] = in.Size(k);
     //      sz[k] = in.Size(k)*2-1;
-  }
-  else {
-    for(k=0; k<BASEARRAY_MAXNDIMS; k++)  sz[k] = 1; 
-    T thr = ZeroThreshold();
-    sa_size_t n = in.Size(in.MaxSizeKA());
-    sa_size_t ncs = ( (in[n-1].imag() < -thr) || (in[n-1].imag() > thr) ) ?
-                    ncs = 2*n-1 : ncs = 2*n-2;
-    sz[in.MaxSizeKA()] = ncs;
-  }
-
+    }
+    else {
+      for(k=0; k<BASEARRAY_MAXNDIMS; k++)  sz[k] = 1; 
+      T thr = ZeroThreshold();
+      sa_size_t n = in.Size(in.MaxSizeKA());
+      sa_size_t ncs = ( (in[n-1].imag() < -thr) || (in[n-1].imag() > thr) ) ?
+        ncs = 2*n-1 : ncs = 2*n-2;
+      sz[in.MaxSizeKA()] = ncs;
+    }
   out.ReSize(in.NbDimensions(), sz);
+  }
 
   return(ndg1);
@@ -226 +241 @@
 }
 
-/* --Methode-- 
-template <class T>
-void FFTArrayChecker<T>::ReShapetoReal( TArray< complex<T> > const & in, TArray< T >  & out)
-{
-  sa_size_t n = in.Size();
-  //  int ncs = (fabs(in(n-1).imag()) > 1.e-12) ? ncs = 2*n-1 : ncs = n*2-2;
-  sa_size_t ncs = 2*n-1;
-  sa_size_t k;
-  out[0] = in[0].real();
-  for(k=1;k<n-1;k++) {
-    out[2*k-1] = in[k].real();
-    out[2*k] = in[k].imag();
-  }
-  //  if (ncs == n*2-2)  out[ncs-1] = in[n-1].real();
-  //  else { out[ncs-2] = in[n-1].real();  out[ncs-1] = in[n-1].imag(); }
-  out[ncs-2] = in[n-1].real();  out[ncs-1] = in[n-1].imag(); 
-}
-*/
-
-/* --Methode-- 
-template <class T>
-void FFTArrayChecker<T>::ReShapetoCompl(TArray< T > const & in, TArray< complex<T> > & out)
-{
-  sa_size_t n =  in.Size();
-  sa_size_t ncs = n/2+1;
-  sa_size_t nc = (n%2 != 0) ? n/2+1 : n/2;
-  out[0] = complex<T> (in[0],0.);
-  for(int k=1;k<nc;k++) 
-    out[k] =  complex<T> (in[2*k-1], in[2*k]);
-  if (n%2 == 0) out[ncs-1] = complex<r_8>(in[n-1], 0.);
-}
-*/
 
 #ifdef __CXX_PRAGMA_TEMPLATES__