- Timestamp:
- Aug 21, 2014, 9:19:24 AM (10 years ago)
- Location:
- trunk/BNC/src/PPP
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/PPP/pppClient.cpp
r5914 r6020 206 206 ColumnVector& xyzc, bool print) { 207 207 208 t_lc::type tLC = (OPT->dualFreqRequired('G') ? t_lc::cIF : t_lc::c1);208 t_lc::type tLC = t_lc::dummy; 209 209 210 210 while (true) { … … 213 213 for (unsigned ii = 0; ii < obsVector.size(); ii++) { 214 214 const t_pppSatObs* satObs = obsVector.at(ii); 215 if ( satObs->isValid() && satObs->prn().system() == 'G' && 216 (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) { 217 ++iObs; 218 BB[iObs][0] = satObs->xc()[0]; 219 BB[iObs][1] = satObs->xc()[1]; 220 BB[iObs][2] = satObs->xc()[2]; 221 BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC); 215 if (satObs->prn().system() == 'G') { 216 if (tLC == t_lc::dummy) { 217 tLC = satObs->isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1; 218 } 219 if ( satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) { 220 ++iObs; 221 BB[iObs][0] = satObs->xc()[0]; 222 BB[iObs][1] = satObs->xc()[1]; 223 BB[iObs][2] = satObs->xc()[2]; 224 BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC); 225 } 222 226 } 223 227 } … … 276 280 double t_pppClient::cmpOffGG(vector<t_pppSatObs*>& obsVector) { 277 281 278 t_lc::type tLC = (OPT->dualFreqRequired('R') ? t_lc::cIF : t_lc::c1);282 t_lc::type tLC = t_lc::dummy; 279 283 double offGG = 0.0; 280 284 … … 289 293 for (unsigned ii = 0; ii < obsVector.size(); ii++) { 290 294 t_pppSatObs* satObs = obsVector.at(ii); 291 if ( satObs->isValid() && satObs->prn().system() == 'R' && 292 (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) { 293 double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC); 294 ++nObs; 295 offGG += ll; 296 if (fabs(ll) > fabs(maxRes)) { 297 maxRes = ll; 298 maxResIndex = ii; 299 maxResPrn = satObs->prn(); 295 if (satObs->prn().system() == 'R') { 296 if (tLC == t_lc::dummy) { 297 tLC = satObs->isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1; 298 } 299 if (satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle)) { 300 double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC); 301 ++nObs; 302 offGG += ll; 303 if (fabs(ll) > fabs(maxRes)) { 304 maxRes = ll; 305 maxResIndex = ii; 306 maxResPrn = satObs->prn(); 307 } 300 308 } 301 309 } -
trunk/BNC/src/PPP/pppSatObs.h
r6018 r6020 17 17 ~t_pppSatObs(); 18 18 bool isValid() const; 19 bool isValid(t_lc::type tLC) const; 19 20 const t_prn& prn() const {return _prn;} 20 21 const ColumnVector& xc() const {return _xcSat;}
Note:
See TracChangeset
for help on using the changeset viewer.