Changeset 1525 for tags/asap2.3.0/apps
- Timestamp:
- 03/13/09 11:55:57 (16 years ago)
- Location:
- tags/asap2.3.0
- Files:
-
- 1 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
tags/asap2.3.0/apps/test.cpp
r1407 r1525 10 10 11 11 int main() { 12 const string fname = "/Users/mar637/data/test.asap"; 13 //const string fname = "test/data/tid-t002.rpf"; 14 //STFillerWrapper stf(fname); 15 //cout << "reading ..." << endl; 16 //stf.read(); 17 18 casa::CountedPtr<Scantable> st(new Scantable(fname, casa::Table::Plain)); 12 const string fname0 = "/tmp/ASAP_problem/c082812231.rpf"; 13 ScantableWrapper stw0(0); 14 STFillerWrapper stf0(stw0); 15 cout << "reading ..." << endl; 16 stf0.open(fname0); 17 stf0.read(); 18 stf0.close(); 19 19 20 STMath stm; 21 vector<bool> msk; 22 vector<casa::CountedPtr<Scantable> > l0; 23 l0.push_back(stw0.getCP()); 24 casa::CountedPtr<Scantable> av0 = stm.average(l0,msk, "NONE", "SCAN"); 25 26 const string fname1 = "/tmp/ASAP_problem/c082880100.rpf"; 27 ScantableWrapper stw1(0); 28 STFillerWrapper stf1(stw1); 29 cout << "reading ..." << endl; 30 stf1.open(fname1); 31 stf1.read(); 32 stf1.close(); 33 vector<casa::CountedPtr<Scantable> > l1; 34 l1.push_back(stw1.getCP()); 35 casa::CountedPtr<Scantable> av1 = stm.average(l1,msk, "NONE", "SCAN"); 36 20 37 vector<casa::CountedPtr<Scantable> > on; 21 on.push_back(st); 22 cout << "averaging..." << endl; 23 casa::CountedPtr<Scantable> avst = stm.average(on); 24 /* 25 cout << "quotient" << endl; 26 casa::CountedPtr<Scantable> quot = stm.autoQuotient(avst); 27 cout << quot->summary() << endl; 28 */ 38 on.push_back(av0); 39 on.push_back(av1); 40 41 casa::CountedPtr<Scantable> avst = stm.merge(on); 42 //avst->makePersistent("t.asap"); 29 43 return 0; 30 44
Note:
See TracChangeset
for help on using the changeset viewer.