Changeset 2593 for trunk/src/STGrid.cpp

Timestamp:

07/09/12 16:54:10 (12 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: Yes/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...

New class asapgrid2 is defined. Difference from asapgrid class is
that input and output is a scantable instance, not a table name.

File:

• trunk/src/STGrid.cpp (modified) (6 diffs)

trunk/src/STGrid.cpp

@@ -1185 +1185 @@
 }
 
+void STGrid::table( Table &tab, uInt i )
+{
+  if ( i >= 0 && i < nfile_ )
+    tab = Table( infileList_[i] ) ;
+}
+
 void STGrid::selectData()
 {
@@ -1191 +1197 @@
   tableList_.resize( nfile_ ) ;
   if ( ifno_ == -1 ) {
-    Table taborg( infileList_[0] ) ;
+    //Table taborg( infileList_[0] ) ;
+    Table taborg ;
+    table( taborg, 0 ) ;
     ROScalarColumn<uInt> ifnoCol( taborg, "IFNO" ) ;
     ifno_ = ifnoCol( 0 ) ;
@@ -1198 +1206 @@
   }
   for ( uInt i = 0 ; i < nfile_ ; i++ ) {
-    Table taborg( infileList_[i] ) ;
+    //Table taborg( infileList_[i] ) ;
+    Table taborg ;
+    table( taborg, i ) ;
     TableExprNode node ;
     if ( ifno != -1 || isMultiIF( taborg ) ) {
@@ -1639 +1649 @@
   Table tab ;
   prepareTable( tab, outfile_ ) ;
+  fillTable( tab ) ;
+
+  t1 = mathutil::gettimeofday_sec() ;
+  os << LogIO::DEBUGGING << "saveData: elapsed time is " << t1-t0 << " sec." << LogIO::POST ; 
+
+  return outfile_ ;
+}
+
+void STGrid::prepareTable( Table &tab, String &name ) 
+{
+  Table t( infileList_[0], Table::Old ) ;
+  t.deepCopy( name, Table::New, False, t.endianFormat(), True ) ;
+  tab = Table( name, Table::Update ) ;
+  // 2012/02/13 TN
+  // explicitly copy subtables since no rows including subtables are 
+  // copied by Table::deepCopy with noRows=True
+  TableCopy::copySubTables( tab, t ) ;
+}
+
+void STGrid::fillTable( Table &tab )
+{
   IPosition dshape = data_.shape() ;
   Int nrow = nx_ * ny_ * npol_ ;
@@ -1687 +1718 @@
   data_.freeStorage( data_p, bdata ) ;
 
-  t1 = mathutil::gettimeofday_sec() ;
-  os << LogIO::DEBUGGING << "saveData: elapsed time is " << t1-t0 << " sec." << LogIO::POST ; 
-
   fillMainColumns( tab ) ;
-
-  return outfile_ ;
-}
-
-void STGrid::prepareTable( Table &tab, String &name ) 
-{
-  Table t( infileList_[0], Table::Old ) ;
-  t.deepCopy( name, Table::New, False, t.endianFormat(), True ) ;
-  tab = Table( name, Table::Update ) ;
-  // 2012/02/13 TN
-  // explicitly copy subtables since no rows including subtables are 
-  // copied by Table::deepCopy with noRows=True
-  TableCopy::copySubTables( tab, t ) ;
 }
 
@@ -1779 +1794 @@
 }
 
-}
+// STGrid2
+STGrid2::STGrid2()
+  : STGrid()
+{
+}
+
+STGrid2::STGrid2( const ScantableWrapper &s )
+  : STGrid()
+{
+  setScantable( s ) ;
+}
+
+STGrid2::STGrid2( const vector<ScantableWrapper> &v )
+  : STGrid()
+{
+  setScantableList( v ) ;
+}
+
+void STGrid2::setScantable( const ScantableWrapper &s )
+{
+  nfile_ = 1 ;
+  dataList_.resize( nfile_ ) ;
+  dataList_[0] = s ;
+  infileList_.resize( nfile_ ) ;
+  infileList_[0] = s.getCP()->table().tableName() ;
+}
+
+void STGrid2::setScantableList( const vector<ScantableWrapper> &v )
+{
+  nfile_ = v.size() ;
+  dataList_.resize( nfile_ ) ;
+  infileList_.resize( nfile_ ) ;
+  for ( uInt i = 0 ; i < nfile_ ; i++ ) {
+    dataList_[i] = v[i] ;
+    infileList_[i] = v[i].getCP()->table().tableName() ;
+  }
+}
+
+ScantableWrapper STGrid2::getResultAsScantable()
+{
+  CountedPtr<Scantable> s = new Scantable( Table::Plain ) ;
+  s->setHeader( dataList_[0].getCP()->getHeader() ) ;
+  Table tout, tin ;
+  String subt[] = { "FREQUENCIES", "FOCUS", "WEATHER", 
+                    "TCAL", "MOLECULES", "HISTORY", "FIT" } ;
+  for ( uInt i = 0 ; i < 7 ; i++ ) {
+    tout = s->table().rwKeywordSet().asTable(subt[i]) ;
+    tin = dataList_[0].getCP()->table().rwKeywordSet().asTable(subt[i]) ;
+    TableCopy::copyRows( tout, tin ) ;
+  }
+  fillTable( s->table() ) ;
+  ScantableWrapper sw( s ) ;
+  return sw ;
+}
+
+void STGrid2::table( Table &tab, uInt i ) 
+{
+  if ( i < nfile_ )
+    tab = dataList_[i].getCP()->table() ;
+}
+
+}