Changeset 3003 in Sophya

Timestamp:

Jul 3, 2006, 1:05:40 PM (19 years ago)

Author:

ansari

Message:

Adaptation SphericalTransformServer a uppression const des arguments FFTForward/Backward cmv+Reza 03/07/2006

Location:

trunk/SophyaLib/Samba

Files:

• sphericaltransformserver.cc (modified) (2 diffs)
• sphericaltransformserver.h (modified) (2 diffs)

trunk/SophyaLib/Samba/sphericaltransformserver.cc

@@ -567 +567 @@
     }
   TVector<complex<T> > transformedData(nph);
-  fftIntfPtr_-> FFTForward(datain, transformedData);
+  // Il faut avoir instancie le serveur de FFT avec l'option preserveinput=true
+  fftIntfPtr_-> FFTForward(const_cast<TVector< complex<T> > &>(datain), transformedData);
 
   TVector< complex<T> > dataout(nmmax+1);
@@ -616 +617 @@
   // du vecteur reel entre est nph)
   //   cout << " longueur de datain  = " << nph << endl;
-  fftIntfPtr_-> FFTForward(datain, transformedData);
+  // Il faut avoir instancie le serveur de FFT avec l'option preserveinput=true
+  fftIntfPtr_-> FFTForward(const_cast< TVector<T> &>(datain), transformedData);
   //  cout <<  " taille de la transformee " << transformedData.Size() << endl;
   //  TVector< complex<T> > dataout(nmmax+1);

trunk/SophyaLib/Samba/sphericaltransformserver.h

@@ -20 +20 @@
   SphericalTransformServer()  
     {
-      fftIntfPtr_=new FFTPackServer;
+      fftIntfPtr_=new FFTPackServer(true); // preserveinput = true
       fftIntfPtr_->setNormalize(false);
     };
@@ -27 +27 @@
  /*!
  Set a fft server. The constructor sets a default fft server (fft-pack). So it is not necessary to call this method for a standard use.
+ \warning The FFTServerInterface object should NOT overwrite the input arrays
 */
   void SetFFTServer(FFTServerInterface* srv=NULL)