Changeset 3002 in Sophya for trunk/SophyaLib/NTools/toeplitzMatrix.h

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/toeplitzMatrix.h (modified) (3 diffs)

trunk/SophyaLib/NTools/toeplitzMatrix.h

@@ -79 +79 @@
 
 
+// Reza+CMV du 3/7/2006 : on a instancie le FFTPackServer avec le flag preserveinput=true
+// On peut donc faire en principe les const_cast sans danger 
 inline void transformeeFourier(const TVector<complex<double> >& v, TVector< complex<double> >& tv) const
 {
-  fftIntfPtr_-> FFTForward(v, tv);
+  fftIntfPtr_-> FFTForward(const_cast<TVector<complex<double> > & >(v), tv);
 } 
 inline void transformeeFourier(const TVector<double>& v, TVector< complex<double> >& tv) const
 {
-  fftIntfPtr_-> FFTForward(v, tv);
+  fftIntfPtr_-> FFTForward(const_cast<TVector<double> & >(v), tv);
 } 
 
@@ -92 +94 @@
 {
   int n= tv.Size();
-  fftIntfPtr_-> FFTBackward(tv, v);
+  fftIntfPtr_-> FFTBackward(const_cast<TVector<complex<double> > & >(tv), v);
   v/=n;
 } 
@@ -98 +100 @@
 {
   int n= tv.Size();
-  fftIntfPtr_-> FFTBackward(tv, v);
+  fftIntfPtr_-> FFTBackward(const_cast< TVector<complex<double> > & >(tv), v);
   double fac = 1./(2*(n-1));
   v *= fac;