Changeset 1525 for tags/asap2.3.0/apps/test.cpp

Timestamp:

03/13/09 11:55:57 (15 years ago)

Author:

Malte Marquarding

Message:

Tagged 2.3.0 release

Location:

tags/asap2.3.0

Files:

• . (copied) (copied from trunk)
• apps/test.cpp (modified) (1 diff)

tags/asap2.3.0/apps/test.cpp

@@ -10 +10 @@
 
 int main() {
-  const string fname = "/Users/mar637/data/test.asap";
-  //const string fname = "test/data/tid-t002.rpf";
-  //STFillerWrapper stf(fname);
-  //cout << "reading ..." << endl;
-  //stf.read();
-  
-  casa::CountedPtr<Scantable>  st(new Scantable(fname, casa::Table::Plain));
+  const string fname0 = "/tmp/ASAP_problem/c082812231.rpf";
+  ScantableWrapper stw0(0);
+  STFillerWrapper stf0(stw0);
+  cout << "reading ..." << endl;
+  stf0.open(fname0);
+  stf0.read();
+  stf0.close();
+
   STMath stm;
+  vector<bool> msk;
+  vector<casa::CountedPtr<Scantable> > l0;
+  l0.push_back(stw0.getCP());
+  casa::CountedPtr<Scantable> av0 = stm.average(l0,msk, "NONE", "SCAN");
+
+  const string fname1 = "/tmp/ASAP_problem/c082880100.rpf";
+  ScantableWrapper stw1(0);
+  STFillerWrapper stf1(stw1);
+  cout << "reading ..." << endl;
+  stf1.open(fname1);
+  stf1.read();
+  stf1.close();
+  vector<casa::CountedPtr<Scantable> > l1;
+  l1.push_back(stw1.getCP());
+  casa::CountedPtr<Scantable> av1 = stm.average(l1,msk, "NONE", "SCAN");
+
   vector<casa::CountedPtr<Scantable> > on;
-  on.push_back(st);
-  cout << "averaging..." << endl;
-  casa::CountedPtr<Scantable> avst = stm.average(on);
-  /*
-  cout << "quotient" << endl;
-  casa::CountedPtr<Scantable> quot = stm.autoQuotient(avst);
-  cout << quot->summary() << endl;
-  */
+  on.push_back(av0);
+  on.push_back(av1);
+ 
+  casa::CountedPtr<Scantable> avst = stm.merge(on);
+  //avst->makePersistent("t.asap");
   return 0;