Index: /trunk/BNC/src/bncutils.cpp
===================================================================
--- /trunk/BNC/src/bncutils.cpp	(revision 6163)
+++ /trunk/BNC/src/bncutils.cpp	(revision 6164)
@@ -606,5 +606,5 @@
 //////////////////////////////////////////////////////////////////////////////
 void kalman(const Matrix& AA, const ColumnVector& ll, const DiagonalMatrix& PP, 
-            SymmetricMatrix& QQ, ColumnVector& xx) {
+            SymmetricMatrix& QQ, ColumnVector& dx) {
 
   Tracer tracer("kalman");
@@ -635,5 +635,5 @@
   SymmetricMatrix Hi; Hi << SHi * SHi.t();
 
-  xx += KT.t() * (ll - AA * xx);
+  dx = KT.t() * (ll - AA * xx);
   QQ << (SS.t() * SS);
 }
Index: /trunk/BNC/src/bncutils.h
===================================================================
--- /trunk/BNC/src/bncutils.h	(revision 6163)
+++ /trunk/BNC/src/bncutils.h	(revision 6164)
@@ -97,5 +97,5 @@
 
 void         kalman(const Matrix& AA, const ColumnVector& ll, const DiagonalMatrix& PP, 
-                    SymmetricMatrix& QQ, ColumnVector& xx);
+                    SymmetricMatrix& QQ, ColumnVector& dx);
 
 double       djul(long j1, long m1, double tt);
Index: /trunk/BNC/src/combination/bnccomb.cpp
===================================================================
--- /trunk/BNC/src/combination/bnccomb.cpp	(revision 6163)
+++ /trunk/BNC/src/combination/bnccomb.cpp	(revision 6164)
@@ -612,7 +612,18 @@
     }
 
+    cout << "AA:  " << AA.Nrows() << ' ' << AA.Ncols() << endl;
+    cout << "ll:  " << ll.Nrows() << endl;
+    cout << "PP:  " << PP.Nrows() << ' ' << PP.Ncols() << endl;
+    cout << "_QQ: " << _QQ.Nrows() << ' ' << _QQ.Ncols() << endl;
+    cout << "dx: "  << dx.Nrows() << endl;
+
     dx = 0.0;
     kalman(AA, ll, PP, _QQ, dx);
+
+    cout << "after Kalman" << endl;
+
     ColumnVector vv = ll - AA * dx;
+
+    cout << vv.t() << endl;
 
     int     maxResIndex;
@@ -807,4 +818,11 @@
   int iObs = 0;
 
+  cout << "nObs, nCon, nPar: " << nObs << ' ' << nCon << ' ' << nPar << endl;
+
+  for (int iPar = 1; iPar <= _params.size(); iPar++) {
+    cmbParam* pp = _params[iPar-1];
+    cout << "PAR " << iPar << ' ' << pp->toString().toAscii().data() << endl;
+  }
+
   QVectorIterator<cmbCorr*> itCorr(corrs());
   while (itCorr.hasNext()) {
@@ -816,4 +834,6 @@
     if (corr->_acName == _masterOrbitAC && resCorr.find(prn) == resCorr.end()) {
       resCorr[prn] = new cmbCorr(*corr);
+      cout << "resCor " << corr->_acName.toAscii().data() << ' '
+           << prn.toAscii().data() << endl;
     }
 
@@ -824,4 +844,6 @@
 
     ll(iObs) = corr->_clkCorr._dClk * t_CST::c - DotProduct(AA.Row(iObs), x0);
+    cout << "iObs, AC, prn, ll: " << iObs << ' ' << corr->_acName.toAscii().data() << ' '
+         << corr->_prn.toAscii().data() << ' ' << ll(iObs) << endl;
   }