Changeset 6112 in ntrip

Timestamp:

Sep 8, 2014, 4:32:01 PM (10 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC/src/PPP_free

Files:

• pppClient.cpp (modified) (1 diff)
• pppFilter.cpp (modified) (3 diffs)
• pppFilter.h (modified) (2 diffs)

trunk/BNC/src/PPP_free/pppClient.cpp

@@ -139 +139 @@
     output->_xyzRover[1] = _filter->y();
     output->_xyzRover[2] = _filter->z();
-    // copy(&_filter->Q().data()[0], &_filter->Q().data()[6], _output->_covMatrix);
+    copy(&_filter->Q().data()[0], &_filter->Q().data()[6], output->_covMatrix);
     output->_neu[0]      = 0.0;
     output->_neu[1]      = 0.0;

trunk/BNC/src/PPP_free/pppFilter.cpp

@@ -172 +172 @@
   // ---------------------------------------------
   _epoData_sav = new t_epoData();
+
+  // Some statistics
+  // ---------------
+  _neu.ReSize(3); _neu = 0.0;
+  _numSat = 0;
+  _pDop   = 0.0;
 }
 
@@ -624 +630 @@
   LOG << endl << endl;
 
+  // Compute dilution of precision
+  // -----------------------------
+  cmpDOP(epoData);
+
   // Final Message (both log file and screen)
   // ----------------------------------------
@@ -645 +655 @@
     double ellRef[3];
     xyz2ell(OPT->_xyzAprRover.data(), ellRef);
-    double neu[3];
-    xyz2neu(ellRef, xyz, neu);
+    xyz2neu(ellRef, xyz, _neu.data());
 
     LOG << "  NEU "
-        << setw(8) << setprecision(3) << neu[0] << " "
-        << setw(8) << setprecision(3) << neu[1] << " "
-        << setw(8) << setprecision(3) << neu[2] << endl << endl;
+        << setw(8) << setprecision(3) << _neu[0] << " "
+        << setw(8) << setprecision(3) << _neu[1] << " "
+        << setw(8) << setprecision(3) << _neu[2] << endl << endl;
   }
   else {

trunk/BNC/src/PPP_free/pppFilter.h

@@ -148 +148 @@
   t_irc update(t_epoData* epoData);
   bncTime time()  const {return _time;}
+  const SymmetricMatrix& Q() const {return _QQ;}
+  const ColumnVector& neu() const {return _neu;}
+  int    numSat() const {return _numSat;}
+  double PDOP()   const {return _pDop;}
   double x()      const {return _params[0]->xx;}
   double y()      const {return _params[1]->xx;}
@@ -235 +239 @@
   bncAntex*             _antex;
   t_tides*              _tides;
+  ColumnVector          _neu;
   int                   _numSat;
   double                _pDop;