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

Timestamp:

Nov 11, 2026, 8:57:49 AM (less than one hour ago)

Author:

stuerze

Message:

some PPP updates

File:

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

trunk/BNC/src/PPP/pppFilter.cpp

@@ -37 +37 @@
 ////////////////////////////////////////////////////////////////////////////
 t_pppFilter::t_pppFilter() {
-  _parlist = 0;
+  _parlist     = 0;
+  _numSat      = 0;
+  _fixRatio    = 0.0;
+  _resRmsCode  = 0.0;
+  _resRmsPhase = 0.0;
 }
 
@@ -50 +54 @@
 t_irc t_pppFilter::processEpoch(t_pppObsPool* obsPool) {
 
-  _numSat = 0;
+  _numSat      = 0;
+  _resRmsCode  = 0.0;
+  _resRmsPhase = 0.0;
   const double maxSolGap = 60.0;
   bool setNeuNoiseToZero = false;
@@ -112 +118 @@
   // ----------------------
   cmpDOP(allObs);
+
+  // Post-fit Residuals: RMS per epoch
+  // ---------------------------------
+  cmpResiduals(allObs);
 
   // Ambiguity Resolution
@@ -135 +145 @@
 
   return success;
+}
+
+// RMS of the post-fit residuals of the epoch (set in processSystem)
+////////////////////////////////////////////////////////////////////////////
+void t_pppFilter::cmpResiduals(const vector<t_pppSatObs*>& allObs) {
+
+  // RMS of code and phase residuals
+  // -------------------------------
+  double   sumCode  = 0.0;
+  double   sumPhase = 0.0;
+  unsigned numCode  = 0;
+  unsigned numPhase = 0;
+  for (const t_pppSatObs* obs : allObs) {
+    const vector<t_lc>& LCs = OPT->LCs(obs->prn().system());
+    for (const t_lc& LC : LCs) {
+      if (!obs->hasRes(LC)) {
+        continue;
+      }
+      double res = obs->getRes(LC);
+      if      (LC._type == t_lc::code  || LC._type == t_lc::codeIF) {
+        sumCode  += res * res;
+        numCode  += 1;
+      }
+      else if (LC._type == t_lc::phase || LC._type == t_lc::phaseIF) {
+        sumPhase += res * res;
+        numPhase += 1;
+      }
+    }
+  }
+  _resRmsCode  = (numCode  > 0) ? sqrt(sumCode  / numCode)  : 0.0;
+  _resRmsPhase = (numPhase > 0) ? sqrt(sumPhase / numPhase) : 0.0;
+
+  LOG << string(_epoTime) << " RMS code "
+      << setw(8) << setprecision(4) << _resRmsCode  << " (" << numCode  << ") phase "
+      << setw(8) << setprecision(4) << _resRmsPhase << " (" << numPhase << ")" << endl;
 }