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

Timestamp:

Feb 5, 2000, 6:22:55 PM (26 years ago)

Author:

ansari

Message:

Introduction FFTMayer, debug de FFTPack - Reza 5/2/2000

File:

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

trunk/SophyaLib/NTools/fftservintf.cc

@@ -7 +7 @@
 {
   _info = info;
-  _fgnorm = false;
+  _fgnorm = true;
 }
 
@@ -115 +115 @@
 
 
-// Methodes pour reordonner les donnees 
-
-/* --Methode-- */
-void FFTServerInterface::ReShapetoReal( TVector< complex<r_8> > const & in, TVector< r_8 >  & out)
-{
-  int N = in.NElts();
-  out.ReSize(in.NElts()*2);
-  out(0) = in(0).real();
-  if (in(in.NElts()-1).imag() == 0)  {
-    for(int i=1; i<in.NElts(); i++)
-      out(i) = in(i).real();
-
-    for(int i=1; i<in.NElts(); i++)
-      out(i+in.NElts()-1) = in(in.NElts()-i-1).imag();
-
-  }
-  else    {
-    for(int i=1; i<in.NElts(); i++)
-      out(i) = in(i).real();
-    for(int i=1; i<in.NElts(); i++)
-      out(i+in.NElts()-1) = in(in.NElts()-i).imag();
-  }
-}
-/* --Methode-- */
-void FFTServerInterface::ReShapetoReal( TVector< complex<r_4> > const & in, TVector< r_4 >  & out)
-{
-  int N = in.NElts();
-  out.ReSize(in.NElts()*2);
-  out(0) = in(0).real();
-  if (in(in.NElts()-1).imag() == 0)  {
-    for(int i=1; i<in.NElts(); i++)
-      out(i) = in(i).real();
-
-    for(int i=1; i<in.NElts(); i++)
-      out(i+in.NElts()-1) = in(in.NElts()-i-1).imag();
-
-  }
-  else    {
-    for(int i=1; i<in.NElts(); i++)
-      out(i) = in(i).real();
-    for(int i=1; i<in.NElts(); i++)
-      out(i+in.NElts()-1) = in(in.NElts()-i).imag();
-  }
-}
-
-
-/* --Methode-- */
-void FFTServerInterface::ReShapetoCompl(TVector< r_8 > const & in, TVector< complex<r_8> > & out)
-{
-  int N = in.NElts();
-  out.ReSize(N/2+1);
-  out(0) = complex<r_8> (in(0),0.);
-
-  if(in.NElts()%2 !=0) {
-    for(int k=1;k<=N/2+1;k++) {
-      out(k) =  complex<r_8> (in(k),in(N-k));
-    }
-  }
-  else  {
-    for(int k=1;k<N/2;k++)  {
-      out(k) =  complex<r_8> (in(k),in(N-k));
-    }
-    out(N/2) = complex<r_8> (in(N/2),0.);
-  }
-
-}
-
-/* --Methode-- */
-void FFTServerInterface::ReShapetoCompl(TVector< r_4 > const & in, TVector< complex<r_4> > & out)
-{
-  int N = in.NElts();
-  out.ReSize(N/2+1);
-  out(0) = complex<r_4> (in(0),0.);
-
-  if(in.NElts()%2 !=0) {
-    for(int k=1;k<=N/2+1;k++) {
-      out(k) =  complex<r_4> (in(k),in(N-k));
-    }
-  }
-  else  {
-    for(int k=1;k<N/2;k++)  {
-      out(k) =  complex<r_4> (in(k),in(N-k));
-    }
-    out(N/2) = complex<r_4> (in(N/2),0.);
-  }
-
-}
-