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

Timestamp:

Feb 9, 2001, 6:08:57 PM (25 years ago)

Author:

ansari

Message:

Modifs interface FFTServer , Reza 9/2/2001

File:

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

trunk/SophyaLib/NTools/fftpserver.cc

@@ -28 +28 @@
 FFTPackServer::FFTPackServer()
    : FFTServerInterface("FFTPackServer using extended FFTPack (C-version) package")
-
+  , ckR4(true, true) , ckR8(true, true)
 {
 //the working array and its size for the different
@@ -56 +56 @@
 
 
-void FFTPackServer::FFTForward(TVector< complex<r_8> > const & in, TVector< complex<r_8> > & out)
-{
-  out = in;
-  fftf(out.NElts(), out.Data());
-  if (getNormalize()) out *= (1./(r_8)(in.NElts()));
-}
-
-void FFTPackServer::FFTBackward(TVector< complex<r_8> > const & in, TVector< complex<r_8> > & out)
-{
-  out = in;
-  fftb(out.NElts(), out.Data());
-}
-
-
-
-void FFTPackServer::FFTForward(TVector< complex<r_4> > const & in, TVector< complex<r_4> > & out)
-{
-  out = in;
-  fftf(out.NElts(), out.Data());
-  if (getNormalize()) out *= (1./(r_4)(in.NElts()));
-}
-
-void FFTPackServer::FFTBackward(TVector< complex<r_4> > const & in, TVector< complex<r_4> > & out)
-{
-  out = in;
-  fftb(out.NElts(), out.Data());
-}
-
-void FFTPackServer::FFTForward(TVector< r_4 > const & in, TVector< complex<r_4> > & out)
-{
-  TVector< r_4 > inout(in, false);
-  fftf(inout.NElts(), inout.Data());
-  ReShapetoCompl(inout, out);
-  if (getNormalize()) out *= (1./(r_4)(in.NElts()));
-}
-
-void FFTPackServer::FFTBackward(TVector< complex<r_4> > const & in, TVector< r_4 > & out)
-{
-  ReShapetoReal(in, out);
-  fftb(out.NElts(), out.Data());
-}
-
-
-void FFTPackServer::FFTForward(TVector< r_8 > const & in, TVector< complex<r_8> > & out)
-{
-  TVector< r_8 > inout(in, false);
-  fftf(inout.NElts(), inout.Data());
-  ReShapetoCompl(inout, out);
-  if (getNormalize()) out *= (1./(r_8)(in.NElts()));
-}
-
-void FFTPackServer::FFTBackward(TVector< complex<r_8> > const & in, TVector< r_8 > & out)
-{
-  ReShapetoReal(in, out);
-  fftb(out.NElts(), out.Data());
+void FFTPackServer::FFTForward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out)
+{
+  ckR8.CheckResize(in, out);
+  out = in;
+  fftf(out.Size(), out.Data());
+  if (getNormalize()) out *= (1./(r_8)(in.Size()));
+}
+
+void FFTPackServer::FFTBackward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out)
+{
+  ckR8.CheckResize(in, out);
+  out = in;
+  fftb(out.Size(), out.Data());
+}
+
+
+
+void FFTPackServer::FFTForward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out)
+{
+  ckR4.CheckResize(in, out);  
+  out = in;
+  fftf(out.Size(), out.Data());
+  if (getNormalize()) out *= (1./(r_4)(in.Size()));
+}
+
+void FFTPackServer::FFTBackward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out)
+{
+  ckR4.CheckResize(in, out);  
+  out = in;
+  fftb(out.Size(), out.Data());
+}
+
+void FFTPackServer::FFTForward(TArray< r_4 > const & in, TArray< complex<r_4> > & out)
+{
+  ckR4.CheckResize(in, out);  
+  TArray< r_4 > inout(in, false);
+  fftf(inout.Size(), inout.Data());
+  ckR4.ReShapetoCompl(inout, out);
+  if (getNormalize()) out *= complex<r_4>((1./(r_4)(in.Size())), 0.);
+}
+
+void FFTPackServer::FFTBackward(TArray< complex<r_4> > const & in, TArray< r_4 > & out)
+{
+  ckR4.CheckResize(in, out);    
+  ckR4.ReShapetoReal(in, out);
+  fftb(out.Size(), out.Data());
+}
+
+
+void FFTPackServer::FFTForward(TArray< r_8 > const & in, TArray< complex<r_8> > & out)
+{
+  ckR8.CheckResize(in, out);  
+  TArray< r_8 > inout(in, false);
+  fftf(inout.Size(), inout.Data());
+  ckR8.ReShapetoCompl(inout, out);
+  if (getNormalize()) out *= complex<r_8>((1./(r_8)(in.Size())), 0.);
+}
+
+void FFTPackServer::FFTBackward(TArray< complex<r_8> > const & in, TArray< r_8 > & out)
+{
+  ckR8.CheckResize(in, out);  
+  ckR8.ReShapetoReal(in, out);
+  fftb(out.Size(), out.Data());
 }
 
@@ -208 +216 @@
 }
 
-// Methodes pour reordonner les donnees 
-
-/* --Methode-- */
-void FFTPackServer::ReShapetoReal( TVector< complex<r_8> > const & in, TVector< r_8 >  & out)
-{
-  int n = in.NElts();
-  int ncs = (fabs(in(n-1).imag()) > 1.e-12) ? ncs = 2*n-1 : ncs = n*2-2; 
-  out.ReSize(ncs);
-  int 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(); }
-}
-
-/* --Methode-- */
-void FFTPackServer::ReShapetoReal( TVector< complex<r_4> > const & in, TVector< r_4 >  & out)
-{
-  int n = in.NElts();
-  int ncs = (fabs(in(n-1).imag()) > 1.e-12) ? ncs = 2*n-1 : ncs = n*2-2; 
-  out.ReSize(ncs);
-  int 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(); }
-}
-
-
-/* --Methode-- */
-void FFTPackServer::ReShapetoCompl(TVector< r_8 > const & in, TVector< complex<r_8> > & out)
-{
-  uint_4 n =  in.NElts();
-  uint_4 ncs = n/2+1;
-  uint_4 nc = (n%2 != 0) ? n/2+1 : n/2;
-  out.ReSize(ncs);
-  out(0) = complex<r_8> (in(0),0.);
-  for(int k=1;k<nc;k++) 
-    out(k) =  complex<r_4> (in(2*k-1), in(2*k));
-  if (n%2 == 0) out(ncs-1) = complex<r_8>(in(n-1), 0.);
-
-}
-
-/* --Methode-- */
-void FFTPackServer::ReShapetoCompl(TVector< r_4 > const & in, TVector< complex<r_4> > & out)
-{
-  uint_4 n =  in.NElts();
-  uint_4 ncs = n/2+1;
-  uint_4 nc = (n%2 != 0) ? n/2+1 : n/2;
-  out.ReSize(ncs);
-  out(0) = complex<r_4> (in(0),0.);
-  for(int k=1;k<nc;k++) 
-    out(k) =  complex<r_4> (in(2*k-1), in(2*k));
-  if (n%2 == 0) out(ncs-1) = complex<r_4>(in(n-1), 0.);
-}
+