- Timestamp:
- 04/16/12 21:18:11 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/STMath.cpp
r2466 r2467 3166 3166 freqrange = ifgfreq[2*i+1] ; 3167 3167 } 3168 3168 3169 3169 3170 3170 // print IF groups … … 3355 3355 3356 3356 // save column values in the vector 3357 vector< vector<uInt> > freqTableIdVec( insize ) ;3358 3357 vector< vector<uInt> > freqIdVec( insize ) ; 3359 3358 vector< vector<uInt> > ifNoVec( insize ) ; 3360 3359 for ( uInt itable = 0 ; itable < insize ; itable++ ) { 3361 ScalarColumn<uInt> freqIDs ;3362 freqIDs.attach( newin[itable]->frequencies().table(), "ID" ) ;3363 3360 ifnoCol.attach( newin[itable]->table(), "IFNO" ) ; 3364 3361 freqIDCol.attach( newin[itable]->table(), "FREQ_ID" ) ; 3365 for ( uInt irow = 0 ; irow < newin[itable]->frequencies().table().nrow() ; irow++ ) {3366 freqTableIdVec[itable].push_back( freqIDs( irow ) ) ;3367 }3368 3362 for ( uInt irow = 0 ; irow < newin[itable]->table().nrow() ; irow++ ) { 3369 3363 freqIdVec[itable].push_back( freqIDCol( irow ) ) ; … … 3380 3374 vector<uInt> freqIdUpdate ; 3381 3375 for ( uInt irow = 0 ; irow < rows ; irow++ ) { 3382 uInt ifno = ifNoVec[itable][irow] ; // IFNO is reset by group index 3376 uInt ifno = ifNoVec[itable][irow] ; // IFNO is reset by group index (IF group index) 3383 3377 double minfreq = ifgfreq[2*ifno] ; 3384 3378 double maxfreq = ifgfreq[2*ifno+1] ; … … 3414 3408 nmaxchan = imax - sigres * (int(cfreq) +( ( cfreq - int(cfreq) >= 0.5 ) ? 1 : 0 )) ; 3415 3409 } 3416 //os << " freqrange = [" << abcissa[nminchan] << ", " << abcissa[nmaxchan] << "]"<< LogIO::POST;3410 //os << "channel range (" << irow << "): [" << nminchan << "," << nmaxchan << "]" << LogIO::POST ; 3417 3411 if ( nmaxchan < nminchan ) { 3418 3412 int tmp = nmaxchan ; … … 3570 3564 // combine frequency group 3571 3565 os << "Combine spectra based on frequency grouping" << LogIO::POST ; 3572 //os << "IFNO is renumbered as frequency group ID (see above)" << LogIO::POST ; 3573 os << "IFNO is renumbered as IF group ID (see above)" << LogIO::POST ; 3566 os << "IFNO is renumbered as frequency group ID (see above)" << LogIO::POST ; 3574 3567 vector<string> coordinfo = tmpout->getCoordInfo() ; 3575 3568 oldinfo[0] = coordinfo[0] ;
Note:
See TracChangeset
for help on using the changeset viewer.