Changeset 10356 in ntrip


Ignore:
Timestamp:
Feb 27, 2024, 10:05:36 AM (7 weeks ago)
Author:
stuerze
Message:

minor changes

Location:
trunk/BNC/src
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/pppClient.cpp

    r10333 r10356  
    255255                                  ColumnVector& xyzc, bool print) {
    256256  t_lc::type tLC = t_lc::dummy;
    257   char sysBancroft = 'G';
    258257  int  numBancroft = obsVector.size();
    259 
    260   _usedSystems['G'] = _usedSystems['R'] = _usedSystems['E'] = _usedSystems['C'] = 0;
    261   for (unsigned jj = 0; jj < obsVector.size(); jj++) {
    262     const t_pppSatObs* satObs = obsVector[jj];
    263     char sys = satObs->prn().system();
    264     _usedSystems[sys]++;
    265   }
    266 
    267   if      ((numBancroft = _usedSystems.value('G')) >= _opt->_minObs) {
    268     sysBancroft = 'G';
    269   }
    270   else if ((numBancroft = _usedSystems.value('E')) >= _opt->_minObs) {
    271     sysBancroft = 'E';
    272   }
    273   else if ((numBancroft = _usedSystems.value('C')) >= _opt->_minObs) {
    274     sysBancroft = 'C';
    275   }
    276   else if ((numBancroft = _usedSystems.value('R')) >= _opt->_minObs) {
    277     sysBancroft = 'R';
    278   }
    279   else {
    280     LOG << "t_pppClient::cmpBancroft not enough observations: " << endl;
    281     return failure;
    282   }
    283258
    284259  while (_running) {
     
    287262    for (unsigned ii = 0; ii < obsVector.size(); ii++) {
    288263      const t_pppSatObs* satObs = obsVector.at(ii);
    289       if (satObs->prn().system() == sysBancroft) {
    290         if (tLC == t_lc::dummy) {
    291           if (satObs->isValid(t_lc::cIF)) {
    292             tLC = t_lc::cIF;
    293           }
    294           else if (satObs->isValid(t_lc::c1)) {
    295             tLC = t_lc::c1;
    296           }
    297           else if (satObs->isValid(t_lc::c2)) {
    298             tLC = t_lc::c2;
    299           }
     264      if (tLC == t_lc::dummy) {
     265        if (satObs->isValid(t_lc::cIF)) {
     266          tLC = t_lc::cIF;
    300267        }
    301         if ( satObs->isValid(tLC) &&
    302             (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
    303           ++iObs;
    304           BB[iObs][0] = satObs->xc()[0];
    305           BB[iObs][1] = satObs->xc()[1];
    306           BB[iObs][2] = satObs->xc()[2];
    307           BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC);
     268        else if (satObs->isValid(t_lc::c1)) {
     269          tLC = t_lc::c1;
    308270        }
     271        else if (satObs->isValid(t_lc::c2)) {
     272          tLC = t_lc::c2;
     273        }
     274      }
     275      if ( satObs->isValid(tLC) &&
     276          (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
     277        ++iObs;
     278        BB[iObs][0] = satObs->xc()[0];
     279        BB[iObs][1] = satObs->xc()[1];
     280        BB[iObs][2] = satObs->xc()[2];
     281        BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC);
    309282      }
    310283    }
     
    327300    for (unsigned ii = 0; ii < obsVector.size(); ii++) {
    328301      const t_pppSatObs* satObs = obsVector.at(ii);
    329       char sys = satObs->prn().system();
    330       if (satObs->isValid() && sys == sysBancroft &&
    331           (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
     302      if (satObs->isValid() && (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
    332303        ColumnVector rr = satObs->xc().Rows(1,3) - xyzc.Rows(1,3);
    333304        double res = rr.NormFrobenius() - satObs->obsValue(tLC)
     
    345316        if (!_epoTimeRover.undef()) LOG << string(_epoTimeRover);
    346317        LOG << "\n---------------------------------------------------------------\n";
    347         LOG << string(epoTime) << " BANCROFT " <<   sysBancroft << ": "
     318        LOG << string(epoTime) << " BANCROFT " << ": "
    348319            << setw(14) << setprecision(3) << xyzc[0] << ' '
    349320            << setw(14) << setprecision(3) << xyzc[1] << ' '
  • trunk/BNC/src/PPP/pppFilter.cpp

    r10335 r10356  
    375375
    376376        // Check Pre-Fit Residuals
    377         /* -----------------------
     377        // -----------------------
    378378        else {
    379379          ColumnVector AA(params.size());
     
    402402            resetAmb(obs->prn(), obsVector, tLC);
    403403          }
    404         }*/
     404        }
    405405      }
    406406    }
  • trunk/BNC/src/PPP/pppParlist.cpp

    r10326 r10356  
    4545  _noise    = 0.0;
    4646  _ambInfo  = 0;
    47   double fac = 1.0;
    4847
    4948  switch (_type) {
     
    6867     break;
    6968   case amb:
    70      if (tLC == t_lc::lIF) {fac = 3.0;}
    7169     _epoSpec = false;
    72      _sigma0  = fac * OPT->_aprSigAmb;
     70     _sigma0  = OPT->_aprSigAmb;
    7371     _ambInfo = new t_ambInfo();
    7472     if (obsVector) {
     
    499497  // GLONASS Clock Offset
    500498  // --------------------
    501   if ( _usedSystems.contains('R')  &&
    502       (_usedSystems.contains('G') || _usedSystems.contains('E') || _usedSystems.contains('C'))) {
     499  if (OPT->useSystem('R')) {
    503500    required.push_back(new t_pppParam(t_pppParam::offGlo, t_prn(), t_lc::dummy));
    504   }
    505   else {
    506     PPP_CLIENT->resetOffGlo();
    507501  }
    508502
    509503  // Galileo Clock Offset
    510504  // --------------------
    511   if (_usedSystems.contains('E') && _usedSystems.contains('G')) {
     505  if (OPT->useSystem('E')) {
    512506    required.push_back(new t_pppParam(t_pppParam::offGal, t_prn(), t_lc::dummy));
    513   }
    514   else {
    515     PPP_CLIENT->resetOffGal();
    516507  }
    517508
    518509  // BDS Clock Offset
    519510  // ----------------
    520   if (_usedSystems.contains('C')  && (_usedSystems.contains('G') || _usedSystems.contains('E'))) {
     511  if (OPT->useSystem('C')) {
    521512    required.push_back(new t_pppParam(t_pppParam::offBds, t_prn(), t_lc::dummy));
    522   }
    523   else {
    524     PPP_CLIENT->resetOffBds();
    525513  }
    526514
  • trunk/BNC/src/pppMain.cpp

    r10327 r10356  
    474474    // Some default values
    475475    // -------------------
    476     opt->_aprSigClk       = 300000.0;
    477     opt->_aprSigClkOff    = 300000.0;
    478     opt->_aprSigAmb       = 10000.0;
    479     opt->_aprSigIon       = 1.0e3;
    480     opt->_noiseIon        = 1.0e3;
     476    opt->_aprSigClk       = 3.0e5;
     477    opt->_aprSigClkOff    = 3.0e5;
     478    opt->_aprSigAmb       = 1.0e4;
     479    opt->_aprSigIon       = 1.0e8;
     480    opt->_noiseIon        = 1.0e8;
    481481    opt->_aprSigCodeBias  = 10000.0;
    482482    opt->_noiseCodeBias   = 10000.0;
Note: See TracChangeset for help on using the changeset viewer.