Changeset 2066 in ntrip

Timestamp:

Dec 1, 2009, 2:46:42 PM (14 years ago)

Author:

mervart

Message:

* empty log message *

File:

• trunk/BNC/bncmodel.cpp (modified) (5 diffs)

trunk/BNC/bncmodel.cpp

@@ -185 +185 @@
   double tropDelay = delay_saast(satData->eleSat);
 
-  cout << "tropDelay " << tropDelay << endl;
-
   return satData->rho + _xcBanc(4) - satData->clk + tropDelay;
 }
@@ -200 +198 @@
   double hh =  50.0 * exp(-6.396e-4 * height);
   double ee =  hh / 100.0 * exp(-37.2465 + 0.213166*TT - 0.000256908*TT*TT);
+
+  cout << "Ele: " << 180.0 * Ele / M_PI << " " << height << " " 
+       << pp << " " << TT << " " << hh << " " << ee << endl;
 
   double h_km = height / 1000.0;
@@ -230 +231 @@
   unsigned nObs = epoData->size();
 
+  // Create First-Design Matrix
+  // --------------------------
   _AA.ReSize(nObs, nPar);  // variance-covariance matrix
   _ll.ReSize(nObs);        // tems observed-computed
@@ -251 +254 @@
   }
 
+  // Compute Least-Squares Solution
+  // ------------------------------
   _QQ.ReSize(nPar);
   _QQ << _AA.t() * _AA;
@@ -256 +261 @@
   _dx = _QQ * _AA.t() * _ll;
 
+  // Compute Residuals
+  // -----------------
+  ColumnVector vv = _AA * _dx - _ll;
+
+  cout << setprecision(3) << vv.t(); 
+  cout.flush();
+
+  // Set Solution Vector
+  // -------------------
   _xx.ReSize(nPar);
-
   unsigned iPar = 0;
   QListIterator<bncParam*> itPar(_params);