Index: trunk/BNC/src/PPP_free/pppFilter.cpp
===================================================================
--- trunk/BNC/src/PPP_free/pppFilter.cpp	(revision 6110)
+++ trunk/BNC/src/PPP_free/pppFilter.cpp	(revision 6111)
@@ -1256,2 +1256,35 @@
 }
 
+// 
+////////////////////////////////////////////////////////////////////////////
+void t_pppFilter::cmpDOP(t_epoData* epoData) {
+
+  Tracer tracer("t_pppFilter::cmpDOP");
+
+  _numSat = 0;
+  _pDop   = 0.0;
+
+  if (_params.size() < 4) {
+    return;
+  }
+
+  const unsigned numPar = 4;
+  Matrix AA(epoData->sizeAll(), numPar);
+  QMapIterator<QString, t_satData*> it(epoData->satData);
+  while (it.hasNext()) {
+    it.next();
+    t_satData* satData = it.value();
+    _numSat += 1;
+    for (unsigned iPar = 0; iPar < numPar; iPar++) {
+      AA[_numSat-1][iPar] = _params[iPar]->partial(satData, false);
+    }
+  }
+  if (_numSat < 4) {
+    return;
+  }
+  AA = AA.Rows(1, _numSat);
+  SymmetricMatrix NN; NN << AA.t() * AA;  
+  SymmetricMatrix QQ = NN.i();
+    
+  _pDop = sqrt(QQ(1,1) + QQ(2,2) + QQ(3,3));
+}
Index: trunk/BNC/src/PPP_free/pppFilter.h
===================================================================
--- trunk/BNC/src/PPP_free/pppFilter.h	(revision 6110)
+++ trunk/BNC/src/PPP_free/pppFilter.h	(revision 6111)
@@ -217,4 +217,6 @@
   void bancroft(const Matrix& BBpass, ColumnVector& pos);
 
+  void cmpDOP(t_epoData* epoData);
+
   t_pppClient*          _pppClient;
   bncTime               _time;
@@ -233,4 +235,6 @@
   bncAntex*             _antex;
   t_tides*              _tides;
+  int                   _numSat;
+  double                _pDop;
 };