Changeset 8956 in ntrip for trunk/BNC/src/PPP/pppFilter.h

Timestamp:

Jun 23, 2020, 11:58:46 AM (4 years ago)

Author:

stuerze

Message:

update regarding PPP

File:

• trunk/BNC/src/PPP/pppFilter.h (modified) (5 diffs)

trunk/BNC/src/PPP/pppFilter.h

@@ -22 +22 @@
   t_irc processEpoch(int num);
 
-  void datumTransformation();
+  void datumTransformation(const ColumnVector& xFltOld, const SymmetricMatrix& QFltOld);
 
   const ColumnVector&    x() const {return _xFlt;}
@@ -75 +75 @@
   };
 
-  class t_datumTrafo{
+  class t_datumTrafo {
   public:
     t_datumTrafo () {initIndices();}
@@ -82 +82 @@
     bool firstSystem() {return _firstSys;}
     void updateIndices(int maxObs) {
-      if (_firstSys) {
+      if (firstSystem()) {
         initIndices();
       }
@@ -94 +94 @@
       _AA2.ReSize(maxObs, numPar); _AA2 = 0.0;
     }
-    void prepareAA(Matrix& AA, int _numEpoProcessing, int nPar) {
-      Matrix& Prep = _AA2;
-      if (_numEpoProcessing == 1) {
-        Prep = _AA1;
+    void prepareAA(const Matrix& AA, int numEpoProcessing, int nPar) {
+      Matrix* Prep = &_AA2;
+      if (numEpoProcessing == 1) {
+        Prep = &_AA1;
       }
-      Prep.SubMatrix(_firstRow, _lastRow, 1, nPar) = AA;
+      Prep->SubMatrix(_firstRow, _lastRow, 1, nPar) << AA;
     }
-    Matrix varCov(const SymmetricMatrix& QFlt) {
+
+    Matrix computeTrafoMatrix() {
       Matrix D21 = (_AA2.t() * _AA2).i() * _AA2.t() * _AA1;
-      return D21 * QFlt * D21.t();
+      return D21;
+    }
+
+    static void printMatrix(const Matrix& X, int nRow, int nCol) {
+      for (int rr = 0; rr < nRow; rr++) {
+        for (int cc = 0; cc < nCol; cc++) {
+          cout << setw(7) << setprecision(4) << X[rr][cc] << " ";
+        }
+        cout << endl;      }
+      cout << endl;
     }
     int     _firstRow;
@@ -130 +140 @@
   void predictCovCrdPart(const SymmetricMatrix& QFltOld);
 
+  t_irc addInfiniteNoise(t_pppParam::e_type para);
+
   bncTime         _epoTime;
   t_pppParlist*   _parlist;