Changeset 2591 for trunk/src

Timestamp:

07/09/12 16:11:50 (12 years ago)

Author:

WataruKawasaki

Message:

New Development: Yes

JIRA Issue: Yes CSV-1869

Ready for Test: Yes

Interface Changes: No

What Interface Changed:

Test Programs:

Put in Release Notes: No

Module(s): asap

Description: added a scantable function to get/set moleculesID columns data.

Location:

trunk/src

Files:

• Scantable.cpp (modified) (16 diffs)
• Scantable.h (modified) (2 diffs)
• ScantableWrapper.h (modified) (1 diff)
• python_Scantable.cpp (modified) (1 diff)

trunk/src/Scantable.cpp

@@ -11 +11 @@
 //
 #include <map>
+#include <sys/time.h>
 
 #include <atnf/PKSIO/SrcType.h>
@@ -2444 +2445 @@
 void Scantable::cubicSplineBaseline(const std::vector<bool>& mask, int nPiece, float thresClip, int nIterClip, bool getResidual, const std::string& progressInfo, const bool outLogger, const std::string& blfile)
 {
+  /*************/
+  //clock_t totTimeStart, totTimeEnd, blTimeStart, blTimeEnd, ioTimeStart, ioTimeEnd;
+  double totTimeStart, totTimeEnd, blTimeStart, blTimeEnd, ioTimeStart, ioTimeEnd, msTimeStart, msTimeEnd, seTimeStart, seTimeEnd, otTimeStart, otTimeEnd, prTimeStart, prTimeEnd;
+  double elapseMs = 0.0;
+  double elapseSe = 0.0;
+  double elapseOt = 0.0;
+  double elapsePr = 0.0;
+  double elapseBl = 0.0;
+  double elapseIo = 0.0;
+  //totTimeStart = clock();
+  totTimeStart = mathutil::gettimeofday_sec();
+  /*************/
+
   try {
     ofstream ofs;
@@ -2474 +2488 @@
     //--------------------------------
     for (int whichrow = 0; whichrow < nRow; ++whichrow) {
+      /******************/
+      //ioTimeStart = clock();
+      ioTimeStart = mathutil::gettimeofday_sec();
+      /**/
+      std::vector<float> sp = getSpectrum(whichrow);
+      /**/
+      //ioTimeEnd = clock();
+      ioTimeEnd = mathutil::gettimeofday_sec();
+      elapseIo += (double)(ioTimeEnd - ioTimeStart);
+
+      //blTimeStart = clock();
+      msTimeStart = mathutil::gettimeofday_sec();
+      /******************/
+
       chanMask = getCompositeChanMask(whichrow, mask);
+
+      /**/
+      msTimeEnd = mathutil::gettimeofday_sec();
+      elapseMs += (double)(msTimeEnd - msTimeStart);
+      blTimeStart = mathutil::gettimeofday_sec();
+      /**/
+
       //fitBaseline(chanMask, whichrow, fitter);
       //setSpectrum((getResidual ? fitter.getResidual() : fitter.getFit()), whichrow);
@@ -2480 +2515 @@
       std::vector<float> params(nPiece*4);
       int nClipped = 0;
-      std::vector<float> res = doCubicSplineFitting(getSpectrum(whichrow), chanMask, nPiece, pieceEdges, params, nClipped, thresClip, nIterClip, getResidual);
+      std::vector<float> res = doCubicSplineFitting(sp, chanMask, nPiece, pieceEdges, params, nClipped, thresClip, nIterClip, getResidual);
+
+      /**/
+      blTimeEnd = mathutil::gettimeofday_sec();
+      elapseBl += (double)(blTimeEnd - blTimeStart);
+      seTimeStart = mathutil::gettimeofday_sec();
+      /**/
+
+
       setSpectrum(res, whichrow);
       //
 
+      /**/
+      seTimeEnd = mathutil::gettimeofday_sec();
+      elapseSe += (double)(seTimeEnd - seTimeStart);
+      otTimeStart = mathutil::gettimeofday_sec();
+      /**/
+
       outputFittingResult(outLogger, outTextFile, chanMask, whichrow, coordInfo, hasSameNchan, ofs, "cubicSplineBaseline()", pieceEdges, params, nClipped);
+
+      /**/
+      otTimeEnd = mathutil::gettimeofday_sec();
+      elapseOt += (double)(otTimeEnd - otTimeStart);
+      prTimeStart = mathutil::gettimeofday_sec();
+      /**/
+
       showProgressOnTerminal(whichrow, nRow, showProgress, minNRow);
+
+      /******************/
+      //blTimeEnd = clock();
+      prTimeEnd = mathutil::gettimeofday_sec();
+      elapsePr += (double)(prTimeEnd - prTimeStart);
+      /******************/
     }
     //--------------------------------
@@ -2494 +2556 @@
     throw;
   }
+  /***************/
+  //totTimeEnd = clock();
+  totTimeEnd = mathutil::gettimeofday_sec();
+  //std::cout << "io    : " << elapseIo/CLOCKS_PER_SEC << " (sec.)" << endl;
+  //std::cout << "bl    : " << elapseBl/CLOCKS_PER_SEC << " (sec.)" << endl;
+  //std::cout << "total : " << (double)(totTimeEnd - totTimeStart)/CLOCKS_PER_SEC << " (sec.)" << endl;
+  std::cout << "io    : " << elapseIo << " (sec.)" << endl;
+  std::cout << "ms    : " << elapseMs << " (sec.)" << endl;
+  std::cout << "bl    : " << elapseBl << " (sec.)" << endl;
+  std::cout << "se    : " << elapseSe << " (sec.)" << endl;
+  std::cout << "ot    : " << elapseOt << " (sec.)" << endl;
+  std::cout << "pr    : " << elapsePr << " (sec.)" << endl;
+  std::cout << "total : " << (double)(totTimeEnd - totTimeStart) << " (sec.)" << endl;
+  /***************/
 }
 
 void Scantable::autoCubicSplineBaseline(const std::vector<bool>& mask, int nPiece, float thresClip, int nIterClip, const std::vector<int>& edge, float threshold, int chanAvgLimit, bool getResidual, const std::string& progressInfo, const bool outLogger, const std::string& blfile)
 {
+  /*************
+  clock_t totTimeStart, totTimeEnd, blTimeStart, blTimeEnd, ioTimeStart, ioTimeEnd;
+  double elapseBl = 0.0;
+  double elapseIo = 0.0;
+  totTimeStart = clock();
+  *************/
+
   try {
     ofstream ofs;
@@ -2530 +2613 @@
 
     for (int whichrow = 0; whichrow < nRow; ++whichrow) {
+      /******************
+      ioTimeStart = clock();
+      */
+      std::vector<float> sp = getSpectrum(whichrow);
+      /*
+      ioTimeEnd = clock();
+      elapseIo += (double)(ioTimeEnd - ioTimeStart);
+
+      blTimeStart = clock();
+      ******************/
 
       //-------------------------------------------------------
@@ -2560 +2653 @@
       std::vector<float> params(nPiece*4);
       int nClipped = 0;
-      std::vector<float> res = doCubicSplineFitting(getSpectrum(whichrow), chanMask, nPiece, pieceEdges, params, nClipped, thresClip, nIterClip, getResidual);
+      std::vector<float> res = doCubicSplineFitting(sp, chanMask, nPiece, pieceEdges, params, nClipped, thresClip, nIterClip, getResidual);
       setSpectrum(res, whichrow);
       //
@@ -2566 +2659 @@
       outputFittingResult(outLogger, outTextFile, chanMask, whichrow, coordInfo, hasSameNchan, ofs, "autoCubicSplineBaseline()", pieceEdges, params, nClipped);
       showProgressOnTerminal(whichrow, nRow, showProgress, minNRow);
+
+      /******************
+      blTimeEnd = clock();
+      elapseBl += (double)(blTimeEnd - blTimeStart);
+      ******************/
+
     }
 
@@ -2573 +2672 @@
     throw;
   }
+
+  /***************
+  totTimeEnd = clock();
+  std::cout << "io    : " << elapseIo/CLOCKS_PER_SEC << " (sec.)" << endl;
+  std::cout << "bl    : " << elapseBl/CLOCKS_PER_SEC << " (sec.)" << endl;
+  std::cout << "total : " << (double)(totTimeEnd - totTimeStart)/CLOCKS_PER_SEC << " (sec.)" << endl;
+  ***************/
 }
 
@@ -3381 +3487 @@
 std::vector<bool> Scantable::getCompositeChanMask(int whichrow, const std::vector<bool>& inMask)
 {
+  /****
+  double ms1TimeStart, ms1TimeEnd, ms2TimeStart, ms2TimeEnd;
+  double elapse1 = 0.0;
+  double elapse2 = 0.0;
+
+  ms1TimeStart = mathutil::gettimeofday_sec();
+  ****/
+
   std::vector<bool> mask = getMask(whichrow);
+
+  /****
+  ms1TimeEnd = mathutil::gettimeofday_sec();
+  elapse1 = ms1TimeEnd - ms1TimeStart;
+  std::cout << "ms1   : " << elapse1 << " (sec.)" << endl;
+  ms2TimeStart = mathutil::gettimeofday_sec();
+  ****/
+
   uInt maskSize = mask.size();
   if (inMask.size() != 0) {
@@ -3391 +3513 @@
     }
   }
+
+  /****
+  ms2TimeEnd = mathutil::gettimeofday_sec();
+  elapse2 = ms2TimeEnd - ms2TimeStart;
+  std::cout << "ms2   : " << elapse2 << " (sec.)" << endl;
+  ****/
 
   return mask;
@@ -3444 +3572 @@
 {
   if (outLogger || outTextFile) {
+  /****
+  double ms1TimeStart, ms1TimeEnd, ms2TimeStart, ms2TimeEnd;
+  double elapse1 = 0.0;
+  double elapse2 = 0.0;
+
+  ms1TimeStart = mathutil::gettimeofday_sec();
+  ****/
+
     float rms = getRms(chanMask, whichrow);
+
+  /****
+  ms1TimeEnd = mathutil::gettimeofday_sec();
+  elapse1 = ms1TimeEnd - ms1TimeStart;
+  std::cout << "ot1   : " << elapse1 << " (sec.)" << endl;
+  ms2TimeStart = mathutil::gettimeofday_sec();
+  ****/
+
     String masklist = getMaskRangeList(chanMask, whichrow, coordInfo, hasSameNchan);
     std::vector<bool> fixed;
@@ -3454 +3598 @@
     }
     if (outTextFile) {
-      ofs << formatPiecewiseBaselineParams(edge, params, fixed, rms, nClipped, masklist, whichrow, true) << flush;
-    }
+      ofs << formatPiecewiseBaselineParams(edge, params, fixed, rms, nClipped, masklist, whichrow, true);// << flush;
+    }
+
+  /****
+  ms2TimeEnd = mathutil::gettimeofday_sec();
+  elapse2 = ms2TimeEnd - ms2TimeStart;
+  std::cout << "ot2   : " << elapse2 << " (sec.)" << endl;
+  ****/
+
   }
 }
@@ -3553 +3704 @@
 float Scantable::getRms(const std::vector<bool>& mask, int whichrow)
 {
+  /****
+  double ms1TimeStart, ms1TimeEnd, ms2TimeStart, ms2TimeEnd;
+  double elapse1 = 0.0;
+  double elapse2 = 0.0;
+
+  ms1TimeStart = mathutil::gettimeofday_sec();
+  ****/
+
   Vector<Float> spec;
+
   specCol_.get(whichrow, spec);
+
+  /****
+  ms1TimeEnd = mathutil::gettimeofday_sec();
+  elapse1 = ms1TimeEnd - ms1TimeStart;
+  std::cout << "rm1   : " << elapse1 << " (sec.)" << endl;
+  ms2TimeStart = mathutil::gettimeofday_sec();
+  ****/
 
   float mean = 0.0;
@@ -3569 +3736 @@
   mean /= (float)n;
   smean /= (float)n;
+
+  /****
+  ms2TimeEnd = mathutil::gettimeofday_sec();
+  elapse2 = ms2TimeEnd - ms2TimeStart;
+  std::cout << "rm2   : " << elapse2 << " (sec.)" << endl;
+  ****/
 
   return sqrt(smean - mean*mean);
@@ -3795 +3968 @@
 }
 
+vector<uint> Scantable::getMoleculeIdColumnData() const
+{
+  Vector<uInt> molIds(mmolidCol_.getColumn());
+  vector<uint> res;
+  molIds.tovector(res);
+  return res;
+}
+
+void Scantable::setMoleculeIdColumnData(const std::vector<uint>& molids)
+{
+  Vector<uInt> molIds(molids);
+  Vector<uInt> arr(mmolidCol_.getColumn());
+  if ( molIds.nelements() != arr.nelements() )
+    throw AipsError("The input data size must be the number of rows.");
+  mmolidCol_.putColumn(molIds);
+}
+
 
 }

trunk/src/Scantable.h

@@ -590 +590 @@
                                             int whichrow, 
                                             bool verbose=false) const;
+  std::vector<uint> getMoleculeIdColumnData() const;
+  void setMoleculeIdColumnData(const std::vector<uint>& molids);
 
 
@@ -766 +768 @@
   void applyChanFlag( casa::uInt whichrow, const std::vector<bool>& msk, casa::uChar flagval);
 
-
 };
 

trunk/src/ScantableWrapper.h

@@ -310 +310 @@
   { return table_->execFFT(whichrow, mask, getRealImag, getAmplitudeOnly); }
 
+  std::vector<uint> getMoleculeIdColumnData() const
+  { return table_->getMoleculeIdColumnData(); }
+  void setMoleculeIdColumnData(const std::vector<uint>& molids)
+  { table_->setMoleculeIdColumnData(molids); }
 
 

trunk/src/python_Scantable.cpp

@@ -161 +161 @@
     //.def("_sspline_baseline", &ScantableWrapper::smoothingSplineBaseline)
     //.def("_test_cin", &ScantableWrapper::testCin)
+    .def("_getmolidcol_list", &ScantableWrapper::getMoleculeIdColumnData)
+    .def("_setmolidcol_list", &ScantableWrapper::setMoleculeIdColumnData)
   ;
 };