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

Timestamp:

Apr 12, 2023, 2:39:30 PM (13 months ago)

Author:

stuerze

Message:

minor changes

File:

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

trunk/BNC/src/PPP/pppFilter.cpp

@@ -48 +48 @@
 // Process Single Epoch
 ////////////////////////////////////////////////////////////////////////////
-t_irc t_pppFilter::processEpoch() {
+t_irc t_pppFilter::processEpoch(int num) {
   _numSat = 0;
   const double maxSolGap = 60.0;
+  bool setNeuNoiseToZero = false;
+
+  if (num > 1) {
+    setNeuNoiseToZero = true;
+  }
 
   // Vector of all Observations
@@ -80 +85 @@
     return failure;
   }
-  #ifdef BNC_DEBUG_PPP
+
+#ifdef BNC_DEBUG_PPP
   if (OPT->_obsModelType == OPT->DCMcodeBias ||
       OPT->_obsModelType == OPT->DCMphaseBias) {
@@ -91 +97 @@
   ColumnVector xFltOld = _xFlt;
   SymmetricMatrix QFltOld = _QFlt;
-  bool setNeuNoiseToZero = false;
-  for (int iSys = 0; iSys < usedSystems.size(); iSys++) {
-    char sys = usedSystems[iSys];
-    int num = 0;
-    for (unsigned jj = 0; jj < allObs.size(); jj++) {
-      if (allObs[jj]->prn().system() == sys) {
-        ++num;
-      }
-    }
-//    if (num < OPT->_minObs) {
-//      setNeuNoiseToZero = true;
-//    }
-  }
   setStateVectorAndVarCovMatrix(xFltOld, QFltOld, setNeuNoiseToZero);
 
@@ -687 +680 @@
   // get last epoch
   t_pppObsPool::t_epoch *epoch = _obsPool->lastEpoch();
+#ifdef BNC_DEBUG_PPP
+    LOG << "GET LAST EPOCH" << endl;
+#endif    
   if (!epoch) {
     LOG << "t_pppFilter::datumTransformation: !lastEpoch" << endl;
@@ -692 +688 @@
   }
   _epoTime = epoch->epoTime();
+ 
   LOG.setf(ios::fixed);
-  LOG << string(_epoTime) << "\nDATUM TRANSFORMATION " << endl;
+  LOG  << "DATUM TRANSFORMATION in Epoch  "<<  string(_epoTime) << endl;
 
   vector<t_pppSatObs*> &allObs = epoch->obsVector();
@@ -795 +792 @@
 
 #ifdef BNC_DEBUG_PPP
-  //LOG << "xFltOld:\n" << xFltOld << endl;
-  //LOG << "xFlt   :\n" << _xFlt   << endl;
+//  LOG << "xFltOld:\n" << xFltOld << endl;
+//  LOG << "xFlt   :\n" << _xFlt   << endl;
 #endif
 
@@ -824 +821 @@
   // =================
   _datumTrafo->switchAA();
-
+  
   _obsPool->putEpoch(_epoTime, allObs, pseudoObsIono, refSatMap);
+#ifdef BNC_DEBUG_PPP  
+  LOG << "PUT EPOCH t_pppFilter::datumTransformation " << _epoTime.timestr().c_str() << endl;
+#endif
 
   return success;