Changeset 10031 in ntrip for trunk/BNC/src/PPP/pppFilter.cpp

Timestamp:

Apr 14, 2023, 11:56:49 AM (13 months ago)

Author:

stuerze

Message:

minor changes

File:

• trunk/BNC/src/PPP/pppFilter.cpp (modified) (6 diffs)

trunk/BNC/src/PPP/pppFilter.cpp

@@ -95 +95 @@
   // Status Vector, Variance-Covariance Matrix
   // -----------------------------------------
-  ColumnVector xFltOld = _xFlt;
+  ColumnVector    xFltOld = _xFlt;
   SymmetricMatrix QFltOld = _QFlt;
   setStateVectorAndVarCovMatrix(xFltOld, QFltOld, setNeuNoiseToZero);
@@ -467 +467 @@
         // -----------------------
         else {
-          ColumnVector AA(nPar);
-          for (unsigned iPar = 0; iPar < nPar; iPar++) {
-            const t_pppParam *par = params[iPar];
-            AA[iPar] = par->partial(_epoTime, obs, tLC, refPrn);
-          }
-          double ll = obs->obsValue(tLC) - obs->cmpValue(tLC)
-                    - DotProduct(_x0, AA);
-          double vv = DotProduct(AA, _xFlt) - ll;
-          if (fabs(vv) > SLIP) {
-            LOG << epoTimeStr << " cycle slip detected " << t_lc::toString(tLC)
-                << ' ' << obs->prn().toString() << ' ' << setw(8)
-                << setprecision(4) << vv << endl;
-            if (preProcessing) {
-              _obsPool->setRefSatChangeRequired(sys, true);
-            } else {
+          if (!preProcessing) {
+            ColumnVector AA(nPar);
+            for (unsigned iPar = 0; iPar < nPar; iPar++) {
+              const t_pppParam *par = params[iPar];
+              AA[iPar] = par->partial(_epoTime, obs, tLC, refPrn);
+            }
+            double ll = obs->obsValue(tLC) - obs->cmpValue(tLC)
+                      - DotProduct(_x0, AA);
+            double vv = DotProduct(AA, _xFlt) - ll;
+            if (fabs(vv) > SLIP) {
+              LOG << epoTimeStr << " cycle slip detected " << t_lc::toString(tLC)
+                  << ' ' << obs->prn().toString() << ' ' << setw(8)
+                  << setprecision(4) << vv << endl;
               resetAmb(obs->prn(), obsVector, tLC);
             }
@@ -682 +680 @@
 #ifdef BNC_DEBUG_PPP
     LOG << "GET LAST EPOCH" << endl;
-#endif    
+#endif
   if (!epoch) {
     LOG << "t_pppFilter::datumTransformation: !lastEpoch" << endl;
@@ -688 +686 @@
   }
   _epoTime = epoch->epoTime();
- 
+
   LOG.setf(ios::fixed);
   LOG  << "DATUM TRANSFORMATION in Epoch  "<<  string(_epoTime) << endl;
@@ -785 +783 @@
   }
 
-  ColumnVector xFltOld = _xFlt;
+  ColumnVector    xFltOld = _xFlt;
   SymmetricMatrix QFltOld = _QFlt;
 
@@ -821 +819 @@
   // =================
   _datumTrafo->switchAA();
-  
+
   _obsPool->putEpoch(_epoTime, allObs, pseudoObsIono, refSatMap);
-#ifdef BNC_DEBUG_PPP  
+#ifdef BNC_DEBUG_PPP
   LOG << "PUT EPOCH t_pppFilter::datumTransformation " << _epoTime.timestr().c_str() << endl;
 #endif