Changeset 3002 in Sophya for trunk/SophyaLib/NTools/fftpserver.cc

Timestamp:

Jul 3, 2006, 1:04:51 PM (19 years ago)

Author:

ansari

Message:

Suppression const des arguments FFTForward/Backward + adaptation de Toeplitz, cmv+Reza 03/07/2006

File:

• trunk/SophyaLib/NTools/fftpserver.cc (modified) (9 diffs)

trunk/SophyaLib/NTools/fftpserver.cc

@@ -40 +40 @@
 
 
-FFTPackServer::FFTPackServer()
+//! Constructor - If preserve_input==true, input arrays will not be overwritten.
+FFTPackServer::FFTPackServer(bool preserve_input)
    : FFTServerInterface("FFTPackServer using extended FFTPack (C-version) package")
-  , ckR4("FFTPackServer: ", true, true) , ckR8("FFTPackServer: ", true, true)
+  , ckR4("FFTPackServer: ", true, true) , ckR8("FFTPackServer: ", true, true),
+  _preserve_input(preserve_input)
 {
 //the working array and its size for the different
@@ -66 +68 @@
 FFTServerInterface * FFTPackServer::Clone()
 {
-  return (new FFTPackServer);
-}
-
-
-void FFTPackServer::FFTForward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out)
+  return (new FFTPackServer(_preserve_input) );
+}
+
+
+void FFTPackServer::FFTForward(TArray< complex<r_8> > & in, TArray< complex<r_8> > & out)
 {
   ckR8.CheckResize(in, out);
@@ -78 +80 @@
 }
 
-void FFTPackServer::FFTBackward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out)
+void FFTPackServer::FFTBackward(TArray< complex<r_8> > & in, TArray< complex<r_8> > & out)
 {
   ckR8.CheckResize(in, out);
@@ -86 +88 @@
 
 
-void FFTPackServer::FFTForward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out)
+void FFTPackServer::FFTForward(TArray< complex<r_4> > & in, TArray< complex<r_4> > & out)
 {
   ckR4.CheckResize(in, out);  
@@ -94 +96 @@
 }
 
-void FFTPackServer::FFTBackward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out)
+void FFTPackServer::FFTBackward(TArray< complex<r_4> > & in, TArray< complex<r_4> > & out)
 {
   ckR4.CheckResize(in, out);  
@@ -101 +103 @@
 }
 
-void FFTPackServer::FFTForward(TArray< r_4 > const & in, TArray< complex<r_4> > & out)
-{
-  ckR4.CheckResize(in, out);  
-  TArray< r_4 > inout(in, false);
+void FFTPackServer::FFTForward(TArray< r_4 > & in, TArray< complex<r_4> > & out)
+{
+  ckR4.CheckResize(in, out);
+  bool share = (_preserve_input) ? false : true;
+  TArray< r_4 > inout(in, share);
   fftf(inout.Size(), inout.Data());
   ReShapetoCompl(inout, out);
@@ -110 +113 @@
 }
 
-void FFTPackServer::FFTBackward(TArray< complex<r_4> > const & in, TArray< r_4 > & out,
+void FFTPackServer::FFTBackward(TArray< complex<r_4> > & in, TArray< r_4 > & out,
                                 bool usoutsz)
 {
@@ -119 +122 @@
 
 
-void FFTPackServer::FFTForward(TArray< r_8 > const & in, TArray< complex<r_8> > & out)
+void FFTPackServer::FFTForward(TArray< r_8 > & in, TArray< complex<r_8> > & out)
 {
   ckR8.CheckResize(in, out);  
-  TArray< r_8 > inout(in, false);
+  bool share = (_preserve_input) ? false : true;
+  TArray< r_8 > inout(in, share);
   fftf(inout.Size(), inout.Data());
   ReShapetoCompl(inout, out);
@@ -128 +132 @@
 }
 
-void FFTPackServer::FFTBackward(TArray< complex<r_8> > const & in, TArray< r_8 > & out,
+void FFTPackServer::FFTBackward(TArray< complex<r_8> > & in, TArray< r_8 > & out,
                                 bool usoutsz)
 {