Changeset 1394 in Sophya for trunk/SophyaLib/NTools/fftmserver.cc

Timestamp:

Feb 12, 2001, 6:09:27 PM (25 years ago)

Author:

ansari

Message:

Changement interface FFTServer - Reza 12/2/2001

File:

• trunk/SophyaLib/NTools/fftmserver.cc (modified) (17 diffs)

trunk/SophyaLib/NTools/fftmserver.cc

@@ -7 +7 @@
 FFTMayerServer::FFTMayerServer()
    : FFTServerInterface("FFTMayerServer using extended FFTMayer package")
+  , ckR4("FFTMayerServer: ", true, true) , ckR8("FFTMayerServer: ", true, true)
 
 {
@@ -24 +25 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTForward(TVector< complex<r_8> > const & in, TVector< complex<r_8> > & out)
-{
+void FFTMayerServer::FFTForward(TArray< complex<r_8> > const & ina, TArray< complex<r_8> > & outa)
+{
+ ckR8.CheckResize(ina, outa);
+ TVector< complex<r_8> > in(ina);
+ TVector< complex<r_8> > out(outa);
+
  r_8 a,b,c,d;
  r_8 q,r,s,t;
@@ -55 +60 @@
  fht_r8(inoutim.Data(),inoutim.NElts());
 
- out.ReSize(in.NElts());
+ // out.ReSize(in.NElts());
  r_8 fn = 1./n;
  if (getNormalize()) 
@@ -66 +71 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTBackward(TVector< complex<r_8> > const & in, TVector< complex<r_8> > & out)
-{
+void FFTMayerServer::FFTBackward(TArray< complex<r_8> > const & ina, TArray< complex<r_8> > & outa)
+{
+ ckR8.CheckResize(ina, outa);
+ TVector< complex<r_8> > in(ina);
+ TVector< complex<r_8> > out(outa);
+
  r_8 a,b,c,d;
  r_8 q,r,s,t;
@@ -87 +96 @@
  fht_r8(inoutim.Data(),inoutim.NElts());
 
- out.ReSize(in.NElts());
+ // out.ReSize(in.NElts());
  out(0) = complex<r_8>(inoutre(0), inoutim(0));
  out(n/2) = complex<r_8>(inoutre(n/2), inoutim(n/2));
@@ -102 +111 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTForward(TVector< complex<r_4> > const & in, TVector< complex<r_4> > & out)
-{
+void FFTMayerServer::FFTForward(TArray< complex<r_4> > const & ina, TArray< complex<r_4> > & outa)
+{
+ ckR4.CheckResize(ina, outa);
+ TVector< complex<r_4> > in(ina);
+ TVector< complex<r_4> > out(outa);
+
  r_4 a,b,c,d;
  r_4 q,r,s,t;
@@ -133 +146 @@
  fht_r4(inoutim.Data(),inoutim.NElts());
 
- out.ReSize(in.NElts());
+ // out.ReSize(in.NElts());
  r_4 fn = 1./n;
  if (getNormalize()) 
@@ -145 +158 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTBackward(TVector< complex<r_4> > const & in, TVector< complex<r_4> > & out)
-{
+void FFTMayerServer::FFTBackward(TArray< complex<r_4> > const & ina, TArray< complex<r_4> > & outa)
+{
+ ckR4.CheckResize(ina, outa);
+ TVector< complex<r_4> > in(ina);
+ TVector< complex<r_4> > out(outa);
+
  r_4 a,b,c,d;
  r_4 q,r,s,t;
@@ -167 +184 @@
  fht_r4(inoutim.Data(),inoutim.NElts());
 
- out.ReSize(in.NElts());
+ // out.ReSize(in.NElts());
  out(0) = complex<r_4>(inoutre(0), inoutim(0));
  out(n/2) = complex<r_4>(inoutre(n/2), inoutim(n/2));
@@ -182 +199 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTForward(TVector< r_4 > const & in, TVector< complex<r_4> > & out)
-{
+void FFTMayerServer::FFTForward(TArray< r_4 > const & ina, TArray< complex<r_4> > & outa)
+{
+ ckR4.CheckResize(ina, outa);
+ TVector< r_4 > in(ina);
+ TVector< complex<r_4> > out(outa);
+
  r_4 a,b;
  int i,j,k;
@@ -205 +226 @@
 
   
- out.ReSize(n/2+1);
+ // out.ReSize(n/2+1);
  out(0) = complex<r_4>(inout(0), 0.);
 
@@ -225 +246 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTBackward(TVector< complex<r_4> > const & in, TVector< r_4 > & out)
-{
+void FFTMayerServer::FFTBackward(TArray< complex<r_4> > const & ina, TArray< r_4 > & outa)
+{
+ ckR4.CheckResize(ina, outa);
+ TVector< complex<r_4> > in(ina);
+ TVector< r_4 > out(outa);
+
  r_4 a,b;
  int i,j,k;
@@ -238 +263 @@
  }
 
- out.ReSize(nc);
+ // out.ReSize(nc);
  out(0) = in(0).real();
  if (nc%2 == 0) { // nc pair
@@ -267 +292 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTForward(TVector< r_8 > const & in, TVector< complex<r_8> > & out)
-{
+void FFTMayerServer::FFTForward(TArray< r_8 > const & ina, TArray< complex<r_8> > & outa)
+{
+ ckR8.CheckResize(ina, outa);
+ TVector< r_8 > in(ina);
+ TVector< complex<r_8> > out(outa);
+
  r_8 a,b;
  int i,j,k;
@@ -290 +319 @@
 
   
- out.ReSize(n/2+1);
+ // out.ReSize(n/2+1);
  out(0) = complex<r_8>(inout(0), 0.);
 
@@ -309 +338 @@
 
 /* --Methode-- */
-void FFTMayerServer::FFTBackward(TVector< complex<r_8> > const & in, TVector< r_8 > & out)
-{
+void FFTMayerServer::FFTBackward(TArray< complex<r_8> > const & ina, TArray< r_8 > & outa)
+{
+ ckR8.CheckResize(ina, outa);
+ TVector< r_8 > out(outa);
+ TVector< complex<r_8> > in(ina);
+
  r_8 a,b;
  int i,j,k;
@@ -322 +355 @@
  }
 
- out.ReSize(nc);
+ // out.ReSize(nc);
  out(0) = in(0).real();
  if (nc%2 == 0) { // nc pair