Changeset 2263 for branches/parallel/src/NROFiller.cpp

Timestamp:

08/04/11 15:42:43 (13 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: No

Ready for Test: No

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...

Tuning NROFiller.
Filler is about 6 times faster than before (~660sec->~110sec).

File:

• branches/parallel/src/NROFiller.cpp (modified) (6 diffs)

branches/parallel/src/NROFiller.cpp

@@ -100 +100 @@
   String srcname ;
   String fieldname ;
-  Array<Float> spectra ;
-  Array<uChar> flagtra ;
-  Array<Float> tsys ;
-  Array<Double> direction ;
+  Vector<Float> spectra ;
+  Vector<uChar> flagtra ;
+  Vector<Float> tsys ;
+  Vector<Double> direction ;
   Float azimuth ;
   Float elevation ;
-  Float parangle ;
+  Float parangle = 0.0 ;
   Float opacity ;
   uInt tcalid ;
@@ -117 +117 @@
   Float winddir ;
   Double srcvel ;
-  Array<Double> propermotion ;
+  Vector<Double> propermotion( 2, 0.0 ) ;
   Vector<Double> srcdir ;
-  Array<Double> scanrate ;
+  Vector<Double> scanrate( 2, 0.0 ) ;
   Int rowCount = 0 ;
 
@@ -126 +126 @@
   Vector<Float> defaultTcal( 1, 1.0 ) ;
   String tcalTime = MVTime( header.utc ).string( MVTime::YMD ) ;
+
+  // TCAL subtable rows
+  setTcal( tcalTime, defaultTcal ) ;
+
+  // FOCUS subtable rows
+  setFocus( parangle ) ;
+
   for ( Int irow = 0 ; irow < (Int)nRow ; irow++ ) {
     // check scan intent
@@ -195 +202 @@
     setMolecule( restfreq ) ;
 
-    // FOCUS subtable row
-    setFocus( parangle ) ;
-
     // WEATHER subtable row 
-    setWeather( temperature, pressure, humidity, windvel, winddir ) ;
-
-    // TCAL subtable row
-    // use default since NRO input is calibrated data
-    setTcal( tcalTime, defaultTcal ) ;
-    
+    setWeather2( temperature, pressure, humidity, windvel, winddir ) ;
 
     // set row attributes
     // SPECTRA, FLAGTRA, and TSYS
-    Vector<Float> spectrum( spectra );
-    Vector<uChar> flags( flagtra ) ;
-    Vector<Float> Tsys( tsys ) ;
-    setSpectrum( spectrum, flags, Tsys ) ;
+    setSpectrum( spectra, flagtra, tsys ) ;
 
     // SCANNO, CYCLENO, IFNO, POLNO, and BEAMNO
-    //uInt ifno = table_->frequencies().addEntry( (Double)fqs[0], (Double)fqs[1], (Double)fqs[2] ) ;
     setIndex( scanno, cycleno, ifno, polno, beamno ) ;
 
@@ -221 +216 @@
 
     // DIRECTION
-    Vector<Double> dir( direction ) ;
-    setDirection(dir, azimuth, elevation ) ;
+    setDirection( direction, azimuth, elevation ) ;
 
     // TIME and INTERVAL
@@ -228 +222 @@
 
     // SRCNAME, SRCTYPE, FIELDNAME, SRCDIRECTION, SRCPROPERMOTION, and SRCVELOCITY
-    Vector<Double> propermot( propermotion ) ;
-    setSource( srcname, srcType, fieldname, srcdir, propermot, srcvel ) ;
+    setSource( srcname, srcType, fieldname, srcdir, propermotion, srcvel ) ;
 
     // SCANRATE
-    Vector<Double> srate( scanrate ) ;
-    setScanRate( srate ) ;
+    setScanRate( scanrate ) ;
 
     // OPACITY