r334: *** empty log message ***
[ctsim.git] / libctsim / imagefile.cpp
index 0b7eb8192f94690020e16dc0135c5f1e1c3b949c..c21c17f6858facaa7114cd9a92e90b9667dbc35b 100644 (file)
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: imagefile.cpp,v 1.30 2001/01/02 07:18:07 kevin Exp $
+**  $Id: imagefile.cpp,v 1.31 2001/01/02 08:15:02 kevin Exp $
 **
 **  This program is free software; you can redistribute it and/or modify
 **  it under the terms of the GNU General Public License (version 2) as
@@ -866,26 +866,16 @@ ImageFile::fourier (ImageFile& result) const
       complexOut[ix] = new CTSimComplex [m_ny];\r
     \r
     // fourier each x column\r
-#if 1\r
     CTSimComplex* pY = new CTSimComplex [m_ny];\r
     for (ix = 0; ix < m_nx; ix++) {\r
       for (iy = 0; iy < m_ny; iy++) {\r
         double dImag = 0;\r
         if (isComplex())\r
           dImag = vLHSImag[ix][iy];\r
-        pY[iy] = complex<double>(vLHS[ix][iy], dImag);\r
+        pY[iy] = std::complex<double>(vLHS[ix][iy], dImag);\r
       } \r
       ProcessSignal::finiteFourierTransform (pY, complexOut[ix], m_ny,  ProcessSignal::FORWARD);\r
     }\r
-#else\r
-    double* pY = new double [m_ny];\r
-    for (ix = 0; ix < m_nx; ix++) {\r
-      for (iy = 0; iy < m_ny; iy++) {\r
-        pY[iy] = vLHS[ix][iy];\r
-      } \r
-      ProcessSignal::finiteFourierTransform (pY, complexOut[ix], m_ny,  ProcessSignal::FORWARD);\r
-    }\r
-#endif\r
     delete [] pY;\r
     \r
     // fourier each y row\r