Changeset 3441 in ntrip for trunk/BNC

Timestamp:

Sep 22, 2011, 11:22:09 AM (13 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/combination/bnccomb.cpp (modified) (2 diffs)

trunk/BNC/combination/bnccomb.cpp

@@ -399 +399 @@
   }
 
-  Matrix         AA;
-  ColumnVector   ll;
-  DiagonalMatrix PP;
+  SymmetricMatrix QQ_sav = _QQ;
+
+  ColumnVector dx;
   QMap<QString, t_corr*> resCorr;
-
-  if (createAmat(AA, ll, PP, x0, resCorr) != success) {
-    return;
-  }
-
-  ColumnVector dx;
-  SymmetricMatrix QQ_sav = _QQ;
-
+    
   // Update and outlier detection loop
   // ---------------------------------
-  for (int ii = 1; ii < 10; ii++) {
+  while (true) {
+
+    Matrix         AA;
+    ColumnVector   ll;
+    DiagonalMatrix PP;
+
+    if (createAmat(AA, ll, PP, x0, resCorr) != success) {
+      return;
+    }
+
     bncModel::kalman(AA, ll, PP, _QQ, dx);
     ColumnVector vv = ll - AA * dx;
@@ -439 +441 @@
       out << "  Outlier" << endl;
       _QQ = QQ_sav;
-      AA.Row(maxResIndex) = 0.0;
-      ll.Row(maxResIndex) = 0.0;
+      corrs().remove(maxResIndex-1);
     }
     else {