Changeset 9982 in ntrip


Ignore:
Timestamp:
Feb 27, 2023, 5:39:25 PM (4 weeks ago)
Author:
stuerze
Message:

updated to r9974

File:
1 edited

Legend:

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

    r9980 r9982  
    344344      }
    345345      if (preProcessing) {
    346         // for refSats no ambiguity parameters exist
     346        // for refSats no ambiguity parameter exists
    347347        if ((obs->prn() == refPrn)
    348348            && (t_lc::toString(maxOutlierLC) == "l1" ||
     
    354354              << ' ' << setw(8) << setprecision(4) << maxOutlier << endl;
    355355          break;
    356         } //else {obs->setOutlier();}
     356        } else {
     357          obs->setOutlier();
     358        }
    357359      } else {    // fin-processing
    358360        LOG << epoTimeStr << " Outlier " << t_lc::toString(maxOutlierLC) << ' '
     
    360362            << maxOutlier << endl;
    361363        if (par) {
    362           resetAmb(par->prn(), obsVector, &QSav, &xSav);
     364          //if ( par->ambResetCandidate() || (OPT->_obsModelType == OPT->DCMcodeBias ||    OPT->_obsModelType == OPT->DCMphaseBias) ) {
     365            resetAmb(par->prn(), obsVector, &QSav, &xSav);
     366          //}
     367          //else {
     368          //  par->setAmbResetCandidate();
     369          //  obs->setOutlier();
     370          //}
    363371        }
    364372        else {
     
    451459        }
    452460        // Check Pre-Fit Residuals
    453         /* -----------------------
     461        // -----------------------
    454462        else {
    455           if (OPT->_obsModelType == OPT->DCMcodeBias || OPT->_obsModelType == OPT->DCMphaseBias) {
    456             if (!preProcessing || refPrn != t_prn()) {
    457               return success;
    458             }
    459             ColumnVector AA(nPar);
    460             for (unsigned iPar = 0; iPar < nPar; iPar++) {
    461               const t_pppParam *par = params[iPar];
    462               AA[iPar] = par->partial(_epoTime, obs, tLC, refPrn);
    463             }
    464             double ll = obs->obsValue(tLC) - obs->cmpValue(tLC) - DotProduct(_x0, AA);
    465             double vv = DotProduct(AA, _xFlt) - ll;
    466             if (fabs(vv) > SLIP) {
    467               LOG << epoTimeStr << " cycle slip detected " << t_lc::toString(tLC)
    468                   << ' ' << obs->prn().toString() << ' ' << setw(8)
    469                   << setprecision(4) << vv << endl;
     463          if (refPrn != t_prn()) {
     464            return success;
     465          }
     466          ColumnVector AA(nPar);
     467          for (unsigned iPar = 0; iPar < nPar; iPar++) {
     468            const t_pppParam *par = params[iPar];
     469            AA[iPar] = par->partial(_epoTime, obs, tLC, refPrn);
     470          }
     471          double ll = obs->obsValue(tLC) - obs->cmpValue(tLC)
     472                    - DotProduct(_x0, AA);
     473          double vv = DotProduct(AA, _xFlt) - ll;
     474          if (fabs(vv) > SLIP) {
     475            LOG << epoTimeStr << " cycle slip detected " << t_lc::toString(tLC)
     476                << ' ' << obs->prn().toString() << ' ' << setw(8)
     477                << setprecision(4) << vv << endl;
     478            if (preProcessing) {
    470479              _obsPool->setRefSatChangeRequired(sys, true);
     480            } else {
     481              resetAmb(obs->prn(), obsVector);
    471482            }
    472483          }
    473         }*/
     484        }
    474485      }
    475486    }
     
    753764      continue;
    754765    }
    755     _datumTrafo->updateIndices(sys, iObs + 1); LOG << "AA Ncols/Nrows: " << AA.Ncols() << "/" << AA.Nrows() << "  nPar: "  << nPar << endl;    LOG << "AA.SubMatrix(1 .. " << iObs+1 << " , 1 .. " <<  nPar << ")" << endl;
     766    _datumTrafo->updateIndices(sys, iObs + 1); //LOG << "AA Ncols/Nrows: " << AA.Ncols() << "/" << AA.Nrows() << "  nPar: "  << nPar << endl;    //LOG << "AA.SubMatrix(1 .. " << iObs+1 << " , 1 .. " <<  nPar << ")" << endl;
    756767    if (_datumTrafo->prepareAA(AA.SubMatrix(1, iObs + 1, 1, nPar), 2)
    757768        != success) {
Note: See TracChangeset for help on using the changeset viewer.