Changeset 5800 in ntrip for trunk/BNC

Timestamp:

Aug 4, 2014, 6:26:53 PM (10 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/src/PPP/filter.cpp (modified) (1 diff)

trunk/BNC/src/PPP/filter.cpp

@@ -207 +207 @@
     // ------------------
     //// beg test
-    DiagonalMatrix PP(iObs+1);
-    for (int ii = 1; ii < iObs+1; ii++) {
-      PP(ii,ii) = Sl(ii,ii) * Sl(ii,ii);
-    }
-    bncModel::kalman(AA, ll, PP, _QFlt, _xFlt);
+    DiagonalMatrix PP(iObs+1); PP = 0.0;
+    for (int ii = 1; ii <= iObs+1; ii++) {
+      PP(ii,ii) = 1.0 / (Sl(ii,ii) * Sl(ii,ii));
+    }
+    ColumnVector dx;
+    bncModel::kalman(AA, ll, PP, _QFlt, dx);
+    _xFlt += dx;
     //// end test