Changeset 1616 for branches/alma/src

Timestamp:

08/07/09 15:49:27 (15 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-729, CAS-1147

Ready to Release: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

Use LogIO instead of std::cout and std::cerr.

Location:

branches/alma/src

Files:

• STFiller.cpp (modified) (13 diffs)
• STFitter.cpp (modified) (2 diffs)
• STHeader.cpp (modified) (3 diffs)
• STMath.cpp (modified) (43 diffs)
• Scantable.cpp (modified) (13 diffs)

branches/alma/src/STFiller.cpp

@@ -35 +35 @@
 #include <casa/System/ProgressMeter.h>
 #include <atnf/PKSIO/NROReader.h>
+#include <casa/Logging/LogIO.h>
 
 #include <time.h>
@@ -470 +471 @@
   time( &t0 ) ;
   tm *ttm = localtime( &t0 ) ;
-  
-  cout << "STFiller::openNRO()  Start time = " << t0 
-       << " (" 
-       << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
-       << " " 
-       << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
-       << ")" << endl ;
+  LogIO os( LogOrigin( "STFiller", "openNRO()", WHERE ) ) ;
+//   cout << "STFiller::openNRO()  Start time = " << t0 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+  os << "Start time = " << t0 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << LogIO::POST ;
 
   // fill STHeader
@@ -498 +505 @@
                                 header_->epoch,
                                 header_->poltype ) ) {
-    cout << "STFiller::openNRO()  Failed to get header information." << endl ;
-    return ;
+//     cout << "STFiller::openNRO()  Failed to get header information." << endl ;
+//     return ;
+    throw( AipsError("Failed to get header information.") ) ;
   }
 
@@ -526 +534 @@
   }
 
-  // DEBUG
-  //cout << "STFiller::openNRO()  nIF " << endl ;
-  //
-
   beamOffset_ = 0;
   vector<Bool> beams = nreader_->getBeams() ;
@@ -549 +553 @@
   }
 
-  // DEBUG
-  //cout << "STFiller::openNRO()  nBeam " << endl ;
-  //
-
   header_->nbeam = nBeam_ ;
   header_->nif = nIF_ ;
@@ -558 +558 @@
   // set header
   table_->setHeader( *header_ ) ;
-
-  // DEBUG
-  //cout << "STFiller::openNRO() Velocity Definition = " << nheader->getVDEF() << endl ;
 
   // DEBUG
@@ -566 +563 @@
   time( &t1 ) ;
   ttm = localtime( &t1 ) ;
-  cout << "STFiller::openNRO()  End time = " << t1 
-       << " (" 
-       << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
-       << " " 
-       << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
-       << ")" << endl ;
-  cout << "STFiller::openNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+//   cout << "STFiller::openNRO()  End time = " << t1 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+//   cout << "STFiller::openNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os << "End time = " << t1 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << endl ;
+  os << "Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os.post() ;
   //
 
@@ -584 +589 @@
   time( &t0 ) ;
   tm *ttm = localtime( &t0 ) ;
-  cout << "STFiller::readNRO()  Start time = " << t0 
-       << " (" 
-       << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
-       << " " 
-       << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
-       << ")" << endl ;
+  LogIO os( LogOrigin( "STFiller", "readNRO()", WHERE ) ) ;
+//   cout << "STFiller::readNRO()  Start time = " << t0 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+  os << "Start time = " << t0 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << LogIO::POST ;
   //
 
@@ -630 +642 @@
   Array<Double> scanrate ;
   for ( i = 0 ; i < imax ; i++ ) {
-//     if( nreader_->getDataset()->getRecord( i ) == NULL ) {
-//       cerr << "STFiller::readNRO()  error while reading row " << i << endl ;
-//       return -1 ;
-//     }
-
     string scanType = nreader_->getScanType( i ) ;
     Int srcType = -1 ;
     if ( scanType.compare( 0, 2, "ON") == 0 ) {
-      // cout << "ON srcType: " << i << endl ;
+      // os << "ON srcType: " << i << LogIO::POST ;
       srcType = 0 ;
     }
     else if ( scanType.compare( 0, 3, "OFF" ) == 0 ) {
-      //cout << "OFF srcType: " << i << endl ;
+      //os << "OFF srcType: " << i << LogIO::POST ;
       srcType = 1 ;
     }
     else if ( scanType.compare( 0, 4, "ZERO" ) == 0 ) {
-      //cout << "ZERO srcType: " << i << endl ;
+      //os << "ZERO srcType: " << i << LogIO::POST ;
       srcType = 2 ;
     }
     else {
-      //cout << "Undefined srcType: " << i << endl ;
+      //os << "Undefined srcType: " << i << LogIO::POST ;
       srcType = 3 ;
     }
@@ -691 +698 @@
                                   srcdir,
                                   scanrate ) ) {
-        cerr << "STFiller::readNRO()  Failed to get scan information." << endl ;
-        return 1 ;
+//         cerr << "STFiller::readNRO()  Failed to get scan information." << endl ;
+//         return 1 ;
+        throw( AipsError("Failed to get scan information.") ) ;
       }
 
@@ -716 +724 @@
         int iadd = -1 ;
         for ( uInt iif = 0 ; iif < freqs.size() ; iif++ ) {
-          //cout << "STFiller::readNRO()  freqs[" << iif << "][1] = " << freqs[iif][1] << endl ;
+          //os << "freqs[" << iif << "][1] = " << freqs[iif][1] << LogIO::POST ;
           double fdiff = abs( freqs[iif][1] - fqs[1] ) / freqs[iif][1] ;
-          //cout << "STFiller::readNRO()  fdiff = " << fdiff << endl ;
+          //os << "fdiff = " << fdiff << LogIO::POST ;
           if ( fdiff < 1.0e-8 ) {
             iadd = iif ;
@@ -798 +806 @@
     // DEBUG
     //int rownum = nreader_->getRowNum() ;
-    //cout << "STFiller::readNRO() Finished row " << i << "/" << rownum << endl ;
+    //os << "Finished row " << i << "/" << rownum << LogIO::POST ;
     //
   }
@@ -806 +814 @@
   time( &t1 ) ;
   ttm = localtime( &t1 ) ;
-  cout << "STFiller::readNRO()  Processed " << i << " rows" << endl ;
-  cout << "STFiller::readNRO()  Added " << i - count << " rows (ignored " 
-       << count << " \"ZERO\" scans)" << endl ;
-  cout << "STFiller::readNRO()  End time = " << t1 
-       << " (" 
-       << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
-       << " " 
-       << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
-       << ")" << endl ;
-  cout << "STFiller::readNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+//   cout << "STFiller::readNRO()  Processed " << i << " rows" << endl ;
+//   cout << "STFiller::readNRO()  Added " << i - count << " rows (ignored " 
+//        << count << " \"ZERO\" scans)" << endl ;
+//   cout << "STFiller::readNRO()  End time = " << t1 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+//   cout << "STFiller::readNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os << "Processed " << i << " rows" << endl ;
+  os << "Added " << i - count << " rows (ignored " 
+     << count << " \"ZERO\" scans)" << endl ;
+  os.post() ;
+  os << "End time = " << t1 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << endl ;
+  os << "Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os.post() ;
   //
 

branches/alma/src/STFitter.cpp

@@ -32 +32 @@
 #include <casa/Arrays/ArrayMath.h>
 #include <casa/Arrays/ArrayLogical.h>
+#include <casa/Logging/LogIO.h>
 #include <scimath/Fitting.h>
 #include <scimath/Fitting/LinearFit.h>
@@ -147 +148 @@
     funcs_[0] = new Polynomial<Float>(ncomp);
   } else {
-    cerr << " compiled functions not yet implemented" << endl;
+    //cerr << " compiled functions not yet implemented" << endl;
+    LogIO os( LogOrigin( "Fitter", "setExpression()", WHERE ) ) ;
+    os << LogIO::WARN << " compiled functions not yet implemented" << LogIO::POST;
     //funcs_.resize(1);
     //funcs_[0] = new CompiledFunction<Float>();

branches/alma/src/STHeader.cpp

@@ -37 +37 @@
 #include <casa/Arrays/IPosition.h>
 #include <casa/Quanta/MVTime.h>
+#include <casa/Logging/LogIO.h>
 
-
+#include <sstream>
 
 #include "STDefs.h"
@@ -79 +80 @@
   MVTime mvt(this->utc);
   mvt.setFormat(MVTime::YMD);
-  cout << "Observer: " << this->observer << endl
-       << "Project: " << this->project << endl
-       << "Obstype: " << this->obstype << endl
-       << "Antenna: " << this->antennaname << endl
-       << "Ant. Position: " << this->antennaposition << endl
-       << "Equinox: " << this->equinox << endl
-       << "Freq. ref.: " << this->freqref << endl
-       << "Ref. frequency: " << this->reffreq << endl
-       << "Bandwidth: "  << this->bandwidth << endl
-       << "Time (utc): "
-       << mvt
-       << endl;
+//   cout << "Observer: " << this->observer << endl
+//        << "Project: " << this->project << endl
+//        << "Obstype: " << this->obstype << endl
+//        << "Antenna: " << this->antennaname << endl
+//        << "Ant. Position: " << this->antennaposition << endl
+//        << "Equinox: " << this->equinox << endl
+//        << "Freq. ref.: " << this->freqref << endl
+//        << "Ref. frequency: " << this->reffreq << endl
+//        << "Bandwidth: "  << this->bandwidth << endl
+//        << "Time (utc): "
+//        << mvt
+//        << endl;
+  LogIO os( LogOrigin( "STHeader", "print()", WHERE ) ) ;
+  os << "Observer: " << this->observer << endl
+     << "Project: " << this->project << endl
+     << "Obstype: " << this->obstype << endl
+     << "Antenna: " << this->antennaname << endl
+     << "Ant. Position: " << this->antennaposition << endl
+     << "Equinox: " << this->equinox << endl
+     << "Freq. ref.: " << this->freqref << endl
+     << "Ref. frequency: " << this->reffreq << endl
+     << "Bandwidth: "  << this->bandwidth << endl
+     << "Time (utc): "
+     << mvt
+     << LogIO::POST ;
   //setprecision(10) << this->utc << endl;
 }
@@ -129 +143 @@
 {
    if (n_>0) {
-      cerr << "Source    ID" << endl;
-      for (uInt i=0; i<n_; i++) {
-         cout << setw(11) << source_(i) << ID_(i) << endl;
-      }
+//       cerr << "Source    ID" << endl;
+//       for (uInt i=0; i<n_; i++) {
+//          cout << setw(11) << source_(i) << ID_(i) << endl;
+     LogIO os( LogOrigin( "SDDataDesc", "summary()", WHERE ) ) ;
+     ostringstream oss ;
+     oss << "Source    ID" << endl;
+     for (uInt i=0; i<n_; i++) {
+       oss << setw(11) << source_(i) << ID_(i) << endl;
+     }
+     os << oss.str() << LogIO::POST ;
    }
 }

branches/alma/src/STMath.cpp

@@ -44 +44 @@
 #include <scimath/Functionals/Polynomial.h>
 
+#include <casa/Logging/LogIO.h>
+#include <sstream>
+
 #include "MathUtils.h"
 #include "RowAccumulator.h"
@@ -74 +77 @@
                  const std::string& avmode)
 {
+  LogIO os( LogOrigin( "STMath", "average()", WHERE ) ) ;
   if ( avmode == "SCAN" && in.size() != 1 )
     throw(AipsError("Can't perform 'SCAN' averaging on multiple tables.\n"
@@ -81 +85 @@
   // check if OTF observation
   String obstype = in[0]->getHeader().obstype ;
-  //bool otfscan = false ;
   Double tol = 0.0 ;
   if ( obstype.find( "OTF" ) != String::npos ) {
-    //cout << "OTF scan" << endl ;
-    //otfscan = true ;
     tol = TOL_OTF ;
   }
@@ -156 +157 @@
     for ( Int i = 0 ; i < length ; i++ ) {
       Vector<Double> t = dircol(i).getAngle(Unit(String("rad"))).getValue() ;
-      //cout << setw( 3 ) << count++ << ": " ;
-      //cout << "[" << t[0] << "," << t[1] << "]" <<  endl ;
+      //os << << count++ << ": " ;
+      //os << "[" << t[0] << "," << t[1] << "]" << LogIO::POST ;
       bool adddir = true ;
       for ( uInt j = 0 ; j < dirs.size() ; j++ ) {
@@ -176 +177 @@
     }
     uInt rowNum = dirs.size() ;
-    //cout << "dirs.size()=" << dirs.size() << endl ;
     tout.addRow( rowNum ) ;
     for ( uInt i = 0 ; i < rowNum ; i++ ) {
@@ -233 +233 @@
         }
       }
-      //cout << "removeRows.size()=" << removeRows.size() << endl ;
       if ( removeRows.size() != 0 ) {
-        //cout << "[" ;
-        //for ( uInt irow=0 ; irow<removeRows.size()-1 ; irow++ )
-        //cout << removeRows[irow] << "," ;
-        //cout << removeRows[removeRows.size()-1] << "]" << endl ;
         subt.removeRow( removeRows ) ;
       }
@@ -298 +293 @@
   }
   if (rowstodelete.nelements() > 0) {
-    cout << rowstodelete << endl;
+    //cout << rowstodelete << endl;
+    os << rowstodelete << LogIO::POST ;
     tout.removeRow(rowstodelete);
     if (tout.nrow() == 0) {
@@ -314 +310 @@
   // check if OTF observation
   String obstype = in->getHeader().obstype ;
-  //bool otfscan = false ;
   Double tol = 0.0 ;
   if ( obstype.find( "OTF" ) != String::npos ) {
-    //cout << "OTF scan" << endl ;
-    //otfscan = true ;
     tol = TOL_OTF ;
   }
@@ -428 +421 @@
       Vector<uChar> flags = flagCol.getColumn(Slicer(Slice(0)));
       // mask spectra for different DIRECTION
-      //cout << "irow=" << outrowCount+irow << ": flagged [" ;
       for ( uInt jrow = 0 ; jrow < subt.nrow() ; jrow++ ) {
         Vector<Double> direction = dircol(jrow).getAngle(Unit(String("rad"))).getValue() ;
@@ -437 +429 @@
         //if ( !allNearAbs( t, direction, tol ) ) {
         if ( dd > tol ) {
-          //cout << jrow << " " ;
           flags[jrow] = userflag ;
         }
       }
-      //cout << "]" << endl ;
-      //cout << "flags=" << flags << endl ;
       Vector<Float> outspec(nchan);
       Vector<uChar> outflag(nchan,0);
@@ -835 +824 @@
          //Debug
          //if(noff!=ndiff) cerr<<"noff and ndiff is not equal"<<endl;
+         //LogIO os( LogOrigin( "STMath", "dototalpower()", WHERE ) ) ; 
+         //if(noff!=ndiff) os<<"noff and ndiff is not equal"<<LogIO::POST;
          meanoff = sum(spoff)/noff;
          meandiff = sum(spdiff)/ndiff;
@@ -954 +945 @@
       //Debug
       //cerr<<"Tsys used="<<tsysrefscalar<<endl;
+      //LogIO os( LogOrigin( "STMath", "dosigref", WHERE ) ) ;
+      //os<<"Tsys used="<<tsysrefscalar<<LogIO::POST;
       // fill the result, replay signal tsys by reference tsys
       outintCol.put(i, resint);
@@ -1120 +1113 @@
       //Array<Float> avtsys =  Float(0.5) * (tsys1 + tsys2);
       // cerr<< "Tsys1="<<tsys1<<" Tsys2="<<tsys2<<endl;
+      // LogIO os( LogOrigin( "STMath", "donod", WHERE ) ) ;
+      // os<< "Tsys1="<<tsys1<<" Tsys2="<<tsys2<<LogIO::POST;
       tsys1[0] = sqrt(tsyssq1 + tsyssq2);
       Array<Float> avtsys =  tsys1;
@@ -1185 +1180 @@
   out2->frequencies().getEntry(rp2, rv2, inc2, freqid2);
   //cerr << out1->frequencies().table().nrow() << " " << out2->frequencies().table().nrow() << endl ;
+  //LogIO os( LogOrigin( "STMath", "dofs()", WHERE ) ) ;
+  //os << out1->frequencies().table().nrow() << " " << out2->frequencies().table().nrow() << LogIO::POST ;
   if (rp1==rp2) {
     Double foffset = rv1 - rv2;
     uInt choffset = static_cast<uInt>(foffset/abs(inc2));
     if (choffset >= nchan) {
-      cerr<<"out-band frequency switching, no folding"<<endl;
+      //cerr<<"out-band frequency switching, no folding"<<endl;
+      LogIO os( LogOrigin( "STMath", "dofs()", WHERE ) ) ;
+      os<<"out-band frequency switching, no folding"<<LogIO::POST;
       nofold = True;
     }
@@ -1230 +1229 @@
   // check
   if ( choffset == 0 ) {
-    cerr << "channel offset is zero, no folding" << endl ;
+    LogIO os( LogOrigin( "STMath", "dofold()", WHERE ) ) ;
+    os << "channel offset is zero, no folding" << LogIO::POST ;
     return out ;
   }
   int nchan = ref->nchan() ;
   if ( abs(choffset) >= nchan ) {
-    cerr << "out-band frequency switching, no folding" << endl ;
+    LogIO os( LogOrigin( "STMath", "dofold()", WHERE ) ) ;
+    os << "out-band frequency switching, no folding" << LogIO::POST ;
     return out ;
   }
@@ -2384 +2385 @@
   throw ( casa::AipsError )
 {
+  LogIO os( LogOrigin( "STMath", "new_average()", WHERE ) ) ;
   if ( avmode == "SCAN" && in.size() != 1 )
     throw(AipsError("Can't perform 'SCAN' averaging on multiple tables.\n"
@@ -2390 +2392 @@
   // check if OTF observation
   String obstype = in[0]->getHeader().obstype ;
-  //bool otfscan = false ;
   Double tol = 0.0 ;
   if ( obstype.find( "OTF" ) != String::npos ) {
-    //cout << "OTF scan" << endl ;
-    //otfscan = true ;
     tol = TOL_OTF ;
   }
@@ -2407 +2406 @@
 
     // TEST: do normal average in each table before IF grouping
-    cout << "Do preliminary averaging" << endl ;
+    os << "Do preliminary averaging" << LogIO::POST ;
     vector< CountedPtr<Scantable> > tmpin( insize ) ;
     for ( uInt itable = 0 ; itable < insize ; itable++ ) {
@@ -2415 +2414 @@
 
     // warning
-    cout << "Average spectra with different spectral resolution" << endl ;
-    cout << endl ;
+    os << "Average spectra with different spectral resolution" << LogIO::POST ;
 
     // temporarily set coordinfo
@@ -2447 +2445 @@
     //          ...
     //          [min_idn0, max_idn0, min_idn1, max_idn1, ...]]
-    //cout << "Check IF settings in each table" << endl ;
+    //os << "Check IF settings in each table" << LogIO::POST ;
     vector< vector<uInt> > freqid( insize );
     vector< vector<double> > iffreq( insize ) ;
@@ -2462 +2460 @@
           uInt id = freqIDCol( irow ) ;
           if ( freqid[itable].size() == 0 || count( freqid[itable].begin(), freqid[itable].end(), id ) == 0 ) {
-            //cout << "itable = " << itable << ": IF " << id << " is included in the list" << endl ;
+            //os << "itable = " << itable << ": IF " << id << " is included in the list" << LogIO::POST ;
             vector<double> abcissa = tmpin[itable]->getAbcissa( irow ) ;
             freqid[itable].push_back( id ) ;
@@ -2473 +2471 @@
 
     // debug
-    //cout << "IF settings summary:" << endl ;
+    //os << "IF settings summary:" << endl ;
     //for ( uInt i = 0 ; i < freqid.size() ; i++ ) {
-    //cout << "   Table" << i << endl ;
+    //os << "   Table" << i << endl ;
     //for ( uInt j = 0 ; j < freqid[i].size() ; j++ ) {
-    //cout << "      id = " << freqid[i][j] << " (min,max) = (" << iffreq[i][2*j] << "," << iffreq[i][2*j+1] << ")" << endl ;
+    //os << "      id = " << freqid[i][j] << " (min,max) = (" << iffreq[i][2*j] << "," << iffreq[i][2*j+1] << ")" << endl ;
     //}
     //}
-    //cout << endl ;
+    //os << endl ;
+    //os.post() ;
 
     // IF grouping based on their frequency coverage
@@ -2492 +2491 @@
     // ifgfreq[numgrp*2] 
     // ifgfreq: [min0_grp0, max0_grp0, min1_grp1, max1_grp1, ...]
-    //cout << "IF grouping based on their frequency coverage" << endl ;
+    //os << "IF grouping based on their frequency coverage" << LogIO::POST ;
     vector< vector<uInt> > ifgrp ;
     vector<double> ifgfreq ;
@@ -2576 +2575 @@
 
     // print IF groups
-    cout << "IF Group summary: " << endl ;
-    cout << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: (TABLE_ID, FREQ_ID)" << endl ;
+    ostringstream oss ;
+    oss << "IF Group summary: " << endl ;
+    oss << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: (TABLE_ID, FREQ_ID)" << endl ;
     for ( uInt i = 0 ; i < ifgrp.size() ; i++ ) {
-      cout << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*i] << "," << ifgfreq[2*i+1] << "]: " ;
+      oss << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*i] << "," << ifgfreq[2*i+1] << "]: " ;
       for ( uInt j = 0 ; j < ifgrp[i].size()/2 ; j++ ) {
-        cout << "(" << ifgrp[i][2*j] << "," << ifgrp[i][2*j+1] << ") " ; 
-      }
-      cout << endl ;
-    }
-    cout << endl ;
+        oss << "(" << ifgrp[i][2*j] << "," << ifgrp[i][2*j+1] << ") " ; 
+      }
+      oss << endl ;
+    }
+    oss << endl ;
+    os << oss.str() << LogIO::POST ;
     
     // print frequency group
-    cout << "Frequency Group summary: " << endl ;
-    cout << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: IF_GROUP_ID" << endl ;
+    oss.str("") ;
+    oss << "Frequency Group summary: " << endl ;
+    oss << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: IF_GROUP_ID" << endl ;
     for ( uInt i = 0 ; i < freqgrp.size() ; i++ ) {
-      cout << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*freqgrp[i][0]] << "," << ifgfreq[2*freqgrp[i][freqgrp[i].size()-1]+1] << "]: " ;
+      oss << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*freqgrp[i][0]] << "," << ifgfreq[2*freqgrp[i][freqgrp[i].size()-1]+1] << "]: " ;
       for ( uInt j = 0 ; j < freqgrp[i].size() ; j++ ) {
-        cout << freqgrp[i][j] << " " ;
-      }
-      cout << endl ;
-    }
-    cout << endl ;
+        oss << freqgrp[i][j] << " " ;
+      }
+      oss << endl ;
+    }
+    oss << endl ;
+    os << oss.str() << LogIO::POST ;
 
     // membership check
@@ -2623 +2626 @@
 
     // print membership
+    //oss.str("") ;
     //for ( uInt i = 0 ; i < insize ; i++ ) {
-    //cout << "Table " << i << endl ;
+    //oss << "Table " << i << endl ;
     //for ( uInt j = 0 ; j < groups[i].size() ; j++ ) {
-    //cout << "   FREQ_ID " <<  setw( 2 ) << freqid[i][j] << ": " ;
+    //oss << "   FREQ_ID " <<  setw( 2 ) << freqid[i][j] << ": " ;
     //for ( uInt k = 0 ; k < groups[i][j].size() ; k++ ) {
-    //cout << setw( 2 ) << groups[i][j][k] << " " ;
+    //oss << setw( 2 ) << groups[i][j][k] << " " ;
     //}
-    //cout << endl ;
+    //oss << endl ;
     //}
     //}
+    //os << oss.str() << LogIO::POST ;
 
     // set back coordinfo
@@ -2649 +2654 @@
     vector<uInt> newifids( insize, 0 ) ;
     for ( uInt itable = 0 ; itable < insize ; itable++ ) {
-      //cout << "Table " << setw(2) << itable << ": " ;
+      //os << "Table " << itable << ": " ;
       for ( uInt ifrow = 0 ; ifrow < groups[itable].size() ; ifrow++ ) {
         addrow[itable].push_back( groups[itable][ifrow].size()-1 ) ;
-        //cout << addrow[itable][ifrow] << " " ;
+        //os << addrow[itable][ifrow] << " " ;
       }
       addtable[itable] = *max_element( addrow[itable].begin(), addrow[itable].end() ) ;
-      //cout << "(" << addtable[itable] << ")" << endl ;
+      //os << "(" << addtable[itable] << ")" << LogIO::POST ;
     }
     newin.resize( insize ) ;
@@ -2681 +2686 @@
         }
         string taql = taqlstream.str() ;
-        //cout << "taql = " << taql << endl ;
+        //os << "taql = " << taql << LogIO::POST ;
         STSelector selector = STSelector() ;
         selector.setTaQL( taql ) ;
@@ -2708 +2713 @@
 
     // print IF groups again for debug
-    //cout << "IF Group summary: " << endl ;
-    //cout << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: (TABLE_ID, FREQ_ID)" << endl ;
+    //oss.str( "" ) ;
+    //oss << "IF Group summary: " << endl ;
+    //oss << "   GROUP_ID [FREQ_MIN, FREQ_MAX]: (TABLE_ID, FREQ_ID)" << endl ;
     //for ( uInt i = 0 ; i < ifgrp.size() ; i++ ) {
-    //cout << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*i] << "," << ifgfreq[2*i+1] << "]: " ;
+    //oss << "   GROUP " << setw( 2 ) << i << " [" << ifgfreq[2*i] << "," << ifgfreq[2*i+1] << "]: " ;
     //for ( uInt j = 0 ; j < ifgrp[i].size()/2 ; j++ ) {
-    //cout << "(" << ifgrp[i][2*j] << "," << ifgrp[i][2*j+1] << ") " ; 
+    //oss << "(" << ifgrp[i][2*j] << "," << ifgrp[i][2*j+1] << ") " ; 
     //}
-    //cout << endl ;
+    //oss << endl ;
     //}
-    //cout << endl ;
+    //oss << endl ;
+    //os << oss.str() << LogIO::POST ;
 
     // reset SCANNO and IFNO/FREQ_ID: IF is reset by the result of sortation 
-    cout << "All scan number is set to 0" << endl ;
-    //cout << "All IF number is set to IF group index" << endl ;
-    cout << endl ;
+    os << "All scan number is set to 0" << LogIO::POST ;
+    //os << "All IF number is set to IF group index" << LogIO::POST ;
     insize = newin.size() ;
     for ( uInt itable = 0 ; itable < insize ; itable++ ) {
@@ -2773 +2779 @@
 
     // reset spectra and flagtra: pick up common part of frequency coverage
-    //cout << "Pick common frequency range and align resolution" << endl ;
+    //os << "Pick common frequency range and align resolution" << LogIO::POST ;
     for ( uInt itable = 0 ; itable < insize ; itable++ ) {
       uInt rows = newin[itable]->nrow() ;
@@ -2783 +2789 @@
         double minfreq = ifgfreq[2*ifno] ;
         double maxfreq = ifgfreq[2*ifno+1] ;
-        //cout << "frequency range: [" << minfreq << "," << maxfreq << "]" << endl ;
+        //os << "frequency range: [" << minfreq << "," << maxfreq << "]" << LogIO::POST ;
         vector<double> abcissa = newin[itable]->getAbcissa( irow ) ;
         int nchan = abcissa.size() ;
         double resol = abcissa[1] - abcissa[0] ;
-        //cout << "abcissa range  : [" << abcissa[0] << "," << abcissa[nchan-1] << "]" << endl ;
+        //os << "abcissa range  : [" << abcissa[0] << "," << abcissa[nchan-1] << "]" << LogIO::POST ;
         if ( minfreq <= abcissa[0] ) 
           nminchan = 0 ;
@@ -2802 +2808 @@
           nmaxchan = abcissa.size() - 1 - int(cfreq) - ( ( cfreq - int(cfreq) >= 0.5 ) ? 1 : 0 ) ;
         }
-        //cout << "channel range (" << irow << "): [" << nminchan << "," << nmaxchan << "]" << endl ;
+        //os << "channel range (" << irow << "): [" << nminchan << "," << nmaxchan << "]" << LogIO::POST ;
         if ( nmaxchan > nminchan ) {
           newin[itable]->reshapeSpectrum( nminchan, nmaxchan, irow ) ;
@@ -2829 +2835 @@
     
     // reset spectra and flagtra: align spectral resolution
-    //cout << "Align spectral resolution" << endl ;
+    //os << "Align spectral resolution" << LogIO::POST ;
     // gmaxdnu: the coarsest frequency resolution in the frequency group
     // gmemid: member index that have a resolution equal to gmaxdnu
@@ -2858 +2864 @@
             int nchan = abcissa.size() ;
             double dnu = abcissa[1] - abcissa[0] ;
-            //cout << "GROUP " << igrp << " (" << tableid << "," << rowid << "): nchan = " << nchan << " (minchan = " << minchan << ")" << endl ; 
+            //os << "GROUP " << igrp << " (" << tableid << "," << rowid << "): nchan = " << nchan << " (minchan = " << minchan << ")" << LogIO::POST ; 
             if ( nchan < minchan ) {
               minchan = nchan ;
@@ -2869 +2875 @@
         }
         // regrid spectra in each group
-        cout << "GROUP " << igrp << endl ;
-        cout << "   Channel number is adjusted to " << minchan << endl ;
-        cout << "   Corresponding frequency resolution is " << maxdnu << "Hz" << endl ;
+        os << "GROUP " << igrp << endl ;
+        os << "   Channel number is adjusted to " << minchan << endl ;
+        os << "   Corresponding frequency resolution is " << maxdnu << "Hz" << LogIO::POST ;
         for ( uInt imem = 0 ; imem < ifgrp[igrp].size()/2 ; imem++ ) {
           uInt tableid = ifgrp[igrp][2*imem] ;
           uInt rowid = ifgrp[igrp][2*imem+1] ;
           freqIDCol.attach( newin[tableid]->table(), "FREQ_ID" ) ;
-          //cout << "tableid = " << tableid << " rowid = " << rowid << ": " << endl ;
-          //cout << "   regridChannel applied to " ;
+          //os << "tableid = " << tableid << " rowid = " << rowid << ": " << LogIO::POST ;
+          //os << "   regridChannel applied to " ;
           if ( tableid != reftable ) 
             refreqid = newin[tableid]->frequencies().addEntry( refpixref, refvalref, refinc ) ;
@@ -2883 +2889 @@
             uInt tfreqid = freqIdVec[tableid][irow] ;
             if ( tfreqid == rowid ) {     
-              //cout << irow << " " ;
+              //os << irow << " " ;
               newin[tableid]->regridChannel( minchan, maxdnu, irow ) ;
               freqIDCol.put( irow, refreqid ) ;
@@ -2889 +2895 @@
             }
           }
-          //cout << endl ;
+          //os << LogIO::POST ;
         } 
       }
@@ -2913 +2919 @@
       }
     }
-    cout << endl ;
 
     // set back coordinfo
@@ -2938 +2943 @@
 
     // combine frequency group
-    cout << "Combine spectra based on frequency grouping" << endl ;
-    cout << "IFNO is renumbered as frequency group ID (see above)" << endl ;
+    os << "Combine spectra based on frequency grouping" << LogIO::POST ;
+    os << "IFNO is renumbered as frequency group ID (see above)" << LogIO::POST ;
     vector<string> coordinfo = tmpout->getCoordInfo() ;
     oldinfo[0] = coordinfo[0] ;
@@ -2955 +2960 @@
     }
     string taql = taqlstream.str() ;
-    //cout << "taql = " << taql << endl ;
+    //os << "taql = " << taql << LogIO::POST ;
     STSelector selector = STSelector() ;
     selector.setTaQL( taql ) ;
@@ -3009 +3014 @@
       dircol.attach( iter.table(), "DIRECTION" ) ;
       uInt polno = polnos( 0 ) ;
-      //cout << "POLNO iteration: " << polno << endl ;
+      //os << "POLNO iteration: " << polno << LogIO::POST ;
 //       for ( uInt igrp = 0 ; igrp < freqgrp.size() ; igrp++ ) {
 //      sizes[igrp].resize( freqgrp[igrp].size() ) ;
@@ -3022 +3027 @@
 //            vector<uChar> fvec ;
 //            flag.tovector( fvec ) ;
-//            //cout << "spec.size() = " << svec.size() << " fvec.size() = " << fvec.size() << endl ;
+//            //os << "spec.size() = " << svec.size() << " fvec.size() = " << fvec.size() << LogIO::POST ;
 //            specout[igrp].insert( specout[igrp].end(), svec.begin(), svec.end() ) ;
 //            flagout[igrp].insert( flagout[igrp].end(), fvec.begin(), fvec.end() ) ;
-//            //cout << "specout[" << igrp << "].size() = " << specout[igrp].size() << endl ;
+//            //os << "specout[" << igrp << "].size() = " << specout[igrp].size() << LogIO::POST ;
 //            sizes[igrp][imem] = spec.nelements() ;
 //          }
@@ -3090 +3095 @@
                 vector<uChar> fvec ;
                 flag.tovector( fvec ) ;
-                //cout << "spec.size() = " << svec.size() << " fvec.size() = " << fvec.size() << endl ;
+                //os << "spec.size() = " << svec.size() << " fvec.size() = " << fvec.size() << LogIO::POST ;
                 specout[igrp].insert( specout[igrp].end(), svec.begin(), svec.end() ) ;
                 flagout[igrp].insert( flagout[igrp].end(), fvec.begin(), fvec.end() ) ;
-                //cout << "specout[" << igrp << "].size() = " << specout[igrp].size() << endl ;
+                //os << "specout[" << igrp << "].size() = " << specout[igrp].size() << LogIO::POST ;
               }
             }
@@ -3119 +3124 @@
         if ( ifnoColOut( irow ) == gmemid[igrp] && !updated[igrp] ) {
           uInt freqidOut = freqidColOut( irow ) ;
-          //cout << "freqgrp " << igrp << " freqidOut = " << freqidOut << endl ;
+          //os << "freqgrp " << igrp << " freqidOut = " << freqidOut << LogIO::POST ;
           double refpix ;
           double refval ;

branches/alma/src/Scantable.cpp

@@ -31 +31 @@
 #include <casa/Containers/RecordField.h>
 #include <casa/Utilities/GenSort.h>
+#include <casa/Logging/LogIO.h>
 
 #include <tables/Tables/TableParse.h>
@@ -1067 +1068 @@
 void Scantable::addFit( const STFitEntry& fit, int row )
 {
-  cout << mfitidCol_(uInt(row)) << endl;
+  //cout << mfitidCol_(uInt(row)) << endl;
+  LogIO os( LogOrigin( "Scantable", "addFit()", WHERE ) ) ;
+  os << mfitidCol_(uInt(row)) << LogIO::POST ;
   uInt id = fitTable_.addEntry(fit, mfitidCol_(uInt(row)));
   mfitidCol_.put(uInt(row), id);
@@ -1101 +1104 @@
       Table subt = t( t.col("SCAN") == scanlist[i]+1 );
       if (subt.nrow()==0) {
-        cerr <<"Scan "<<scanlist[i]<<" cannot be found in the scantable."<<endl;
+        //cerr <<"Scan "<<scanlist[i]<<" cannot be found in the scantable."<<endl;
+        LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+        os <<LogIO::WARN<<"Scan "<<scanlist[i]<<" cannot be found in the scantable."<<LogIO::POST;
         ret = 1;
         break;
@@ -1113 +1118 @@
           Table subt2 = t( t.col("SCAN") == scanlist[i+1]+1 );
           if ( subt2.nrow() == 0) {
-            cerr<<"Scan "<<scanlist[i+1]<<" cannot be found in the scantable."<<endl;
+            LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+
+            //cerr<<"Scan "<<scanlist[i+1]<<" cannot be found in the scantable."<<endl;
+            os<<LogIO::WARN<<"Scan "<<scanlist[i+1]<<" cannot be found in the scantable."<<LogIO::POST;
             ret = 1;
             break;
@@ -1123 +1131 @@
           if (scan1seqn == 1 && scan2seqn == 2) {
             if (laston1 == laston2) {
-              cerr<<"A valid scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<endl;
+              LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+              //cerr<<"A valid scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<endl;
+              os<<"A valid scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<LogIO::POST;
               i +=1;
             }
             else {
-              cerr<<"Incorrect scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<endl;
+              LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+              //cerr<<"Incorrect scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<endl;
+              os<<LogIO::WARN<<"Incorrect scan pair ["<<scanlist[i]<<","<<scanlist[i+1]<<"]"<<LogIO::POST;
             }
           }
           else if (scan1seqn==2 && scan2seqn == 1) {
             if (laston1 == laston2) {
-              cerr<<"["<<scanlist[i]<<","<<scanlist[i+1]<<"] is a valid scan pair but in incorrect order."<<endl;
+              LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+              //cerr<<"["<<scanlist[i]<<","<<scanlist[i+1]<<"] is a valid scan pair but in incorrect order."<<endl;
+              os<<LogIO::WARN<<"["<<scanlist[i]<<","<<scanlist[i+1]<<"] is a valid scan pair but in incorrect order."<<LogIO::POST;
               ret = 1;
               break;
@@ -1138 +1152 @@
           }
           else {
-            cerr<<"The other scan for  "<<scanlist[i]<<" appears to be missing. Check the input scan numbers."<<endl;
+            LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+            //cerr<<"The other scan for  "<<scanlist[i]<<" appears to be missing. Check the input scan numbers."<<endl;
+            os<<LogIO::WARN<<"The other scan for  "<<scanlist[i]<<" appears to be missing. Check the input scan numbers."<<LogIO::POST;
             ret = 1;
             break;
@@ -1145 +1161 @@
       }
       else {
-        cerr<<"The scan does not appear to be standard obsevation."<<endl;
+        LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+        //cerr<<"The scan does not appear to be standard obsevation."<<endl;
+        os<<LogIO::WARN<<"The scan does not appear to be standard obsevation."<<LogIO::POST;
       }
     //if ( i >= nscan ) break;
@@ -1151 +1169 @@
   }
   else {
-    cerr<<"No reference to GBT_GO table."<<endl;
+    LogIO os( LogOrigin( "Scantable", "checkScanInfo()", WHERE ) ) ;
+    //cerr<<"No reference to GBT_GO table."<<endl;
+    os<<LogIO::WARN<<"No reference to GBT_GO table."<<LogIO::POST;
     ret = 1;
   }
@@ -1185 +1205 @@
     nmax = nmin ;
     nmin = tmp ;
-    cout << "Swap values. Applied range is [" 
-         << nmin << ", " << nmax << "]" << endl ;
+    LogIO os( LogOrigin( "Scantable", "reshapeSpectrum()", WHERE ) ) ;
+    os << "Swap values. Applied range is [" 
+       << nmin << ", " << nmax << "]" << LogIO::POST ;
   }
   
@@ -1198 +1219 @@
     if ( nmin == 0 ) {
       // nothing to do
-      cout << "Whole range is selected. Nothing to do." << endl ;
+      LogIO os( LogOrigin( "Scantable", "reshapeSpectrum()", WHERE ) ) ;
+      os << "Whole range is selected. Nothing to do." << LogIO::POST ;
       return ;
     }
     else {
-      cout << "Specified maximum exceeds nChan. Applied range is ["
-           << nmin << ", " << nChan-1 << "]." << endl ;
+      LogIO os( LogOrigin( "Scantable", "reshapeSpectrum()", WHERE ) ) ;
+      os << "Specified maximum exceeds nChan. Applied range is ["
+         << nmin << ", " << nChan-1 << "]." << LogIO::POST ;
       nmax = nChan - 1 ;
     }
@@ -1256 +1279 @@
 void asap::Scantable::regridChannel( int nChan, double dnu )
 {
-  cout << "Regrid abcissa with channel number " << nChan << " and spectral resoultion " << dnu << "Hz." << endl ;
+  LogIO os( LogOrigin( "Scantable", "regridChannel()", WHERE ) ) ;
+  os << "Regrid abcissa with channel number " << nChan << " and spectral resoultion " << dnu << "Hz." << LogIO::POST ;
   // assumed that all rows have same nChan
   Vector<Float> arr = specCol_( 0 ) ;
@@ -1263 +1287 @@
   // if oldsize == nChan, nothing to do
   if ( oldsize == nChan ) {
-    cout << "Specified channel number is same as current one. Nothing to do." << endl ;
+    os << "Specified channel number is same as current one. Nothing to do." << LogIO::POST ;
     return ;
   }
@@ -1269 +1293 @@
   // if oldChan < nChan, unphysical operation
   if ( oldsize < nChan ) {
-    cout << "Unphysical operation. Nothing to do." << endl ;
+    os << "Unphysical operation. Nothing to do." << LogIO::POST ;
     return ;
   }