Changeset 9600 in ntrip for trunk/BNC/src/PPP/pppClient.cpp

Timestamp:

Jan 14, 2022, 5:31:00 PM (2 years ago)

Author:

stuerze

Message:

 

File:

• trunk/BNC/src/PPP/pppClient.cpp (modified) (5 diffs)

trunk/BNC/src/PPP/pppClient.cpp

@@ -294 +294 @@
         }
       }
-      if ( satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
+      if ( satObs->isValid(tLC) &&
+          (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
         ++iObs;
         BB[iObs][0] = satObs->xc()[0];
@@ -307 +308 @@
     }
     BB = BB.Rows(1,iObs+1);
-    bancroft(BB, xyzc);
+    if (bancroft(BB, xyzc) != success) {
+      return failure;
+    }
 
     xyzc[3] /= t_CST::c;
@@ -319 +322 @@
       const t_pppSatObs* satObs = obsVector.at(ii);
       if (satObs->isValid() &&
-          satObs->prn().system() == 'G' &&
           (!satObs->modelSet() || satObs->eleSat() >= _opt->_minEle) ) {
         ColumnVector rr = satObs->xc().Rows(1,3) - xyzc.Rows(1,3);
         double res = rr.NormFrobenius() - satObs->obsValue(tLC)
                    - (satObs->xc()[3] - xyzc[3]) * t_CST::c;
-        if (std::isnan(res) || fabs(res) > maxRes) {
-          std::isnan(res) ?
-          maxRes      = res :
-          maxRes      = fabs(res);
+        if (fabs(res) > maxRes) {
+          maxRes = fabs(res);
           maxResIndex = ii;
         }
       }
     }
-    if (!std::isnan(maxRes) && maxRes < BLUNDER) {
+    if (maxRes < BLUNDER) {
       if (print && _numEpoProcessing == 1) {
         LOG.setf(ios::fixed);
@@ -594 +594 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-void t_pppClient::bancroft(const Matrix& BBpass, ColumnVector& pos) {
+t_irc t_pppClient::bancroft(const Matrix& BBpass, ColumnVector& pos) {
 
   if (pos.Nrows() != 4) {
@@ -662 +662 @@
     }
   }
+  if (pos.size() != 4 ||
+      std::isnan(pos(1)) ||
+      std::isnan(pos(2)) ||
+      std::isnan(pos(3)) ||
+      std::isnan(pos(4)))  {
+    return failure;
+  }
+
+  return success;
 }