Changeset 2385
- Timestamp:
- 12/22/11 19:05:30 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/STGrid.cpp
r2384 r2385 197 197 Float *wgrid_p = gwgt.getStorage( deleteWgtG ) ; 198 198 Float *conv_p = convFunc.getStorage( deleteConv ) ; 199 200 // pass copy of irow to ggridsd since it will be modified in theroutine 201 Int irowCopy = irow ; 199 202 200 203 // call ggridsd … … 208 211 wgt_p, 209 212 &nrow, 210 &irow ,213 &irowCopy, 211 214 grid_p, 212 215 wgrid_p, … … 244 247 { 245 248 LogIO os( LogOrigin("STGrid", "grid", WHERE) ) ; 249 double t0,t1 ; 246 250 247 251 // data selection 252 t0 = mathutil::gettimeofday_sec() ; 248 253 selectData() ; 254 t1 = mathutil::gettimeofday_sec() ; 255 os << "selectData: elapsed time is " << t1-t0 << " sec." << LogIO::POST ; 256 249 257 setupArray() ; 250 258 … … 287 295 // TODO: nchunk_ must be determined from nchan_, npol_, and (nx_,ny_) 288 296 // by considering data size to be allocated for ggridsd input/output 289 nchunk_ = 100 ;297 nchunk_ = 400 ; 290 298 Bool b = nchunk_ >= nrow_ ; 291 299 nchunk_ = min( nchunk_, nrow_ ) ; … … 369 377 t0 = mathutil::gettimeofday_sec() ; 370 378 Int nrow = getDataChunk( spectra, direction, flagtra, rflag, weight ) ; 371 //os << "nrow = " << nrow << LogIO::POST ;372 379 t1 = mathutil::gettimeofday_sec() ; 373 380 eGetData += t1-t0 ; … … 498 505 toPixel( direction, xypos ) ; 499 506 t1 = mathutil::gettimeofday_sec() ; 500 //os << "xypos=" << xypos << LogIO::POST ; 501 os << "toPixel: elapsed time is " << t1-t0 << " sec." << LogIO::POST ; 507 eToPixel += t1-t0 ; 502 508 503 509 // call ggridsd … … 546 552 void STGrid::initPol( Int ipol ) 547 553 { 548 ptab_ = tab_( tab_.col("POLNO") == (uInt)ipol ) ; 554 if ( npol_ == 1 ) 555 ptab_ = tab_ ; 556 else 557 ptab_ = tab_( tab_.col("POLNO") == (uInt)ipol ) ; 549 558 550 559 attach( ptab_ ) ; … … 701 710 Table taborg( infile_ ) ; 702 711 if ( ifno == -1 ) { 703 LogIO os( LogOrigin("STGrid","selectData",WHERE) ) ; 712 LogIO os( LogOrigin("STGrid","selectData",WHERE) ) ; 704 713 // os << LogIO::SEVERE 705 714 // << "Please set IFNO before actual gridding" … … 819 828 Array<Float> &weight ) 820 829 { 830 LogIO os( LogOrigin("STGrid","getDataChunk",WHERE) ) ; 821 831 Int nrow = getDataChunk( spectraF_, direction, flagtraUC_, rflagUI_, weight ) ; 822 832 if ( nrow < nchunk_ ) { … … 864 874 865 875 RefRows rows( nprocessed_, nprocessed_+nrow-1, 1 ) ; 866 os << "rows.nrows()=" << rows.nrows() <<LogIO::POST ;876 os<<LogIO::DEBUGGING<<"nprocessed_="<<nprocessed_<<": rows.nrows()="<<rows.nrows()<<LogIO::POST ; 867 877 spectraCol_.getColumnCells( rows, spectra ) ; 868 878 flagtraCol_.getColumnCells( rows, flagtra ) ; … … 875 885 else 876 886 tsys = tsysTemp[0] ; 877 887 878 888 double t0,t1 ; 879 889 t0 = mathutil::gettimeofday_sec() ;
Note:
See TracChangeset
for help on using the changeset viewer.