Changeset 2245 in ntrip

Timestamp:

Jan 12, 2010, 5:13:10 PM (15 years ago)

Author:

mervart

Message:

* empty log message *

File:

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

trunk/BNC/bncmodel.cpp

@@ -633 +633 @@
         ll(iObs)      = satData->L3 - rhoCmp;
 
-        cout.setf(ios::fixed);
-        cout << prn.toAscii().data() << " " 
-             << setprecision(3) << rhoCmp << " "
-             << setprecision(3) << satData->P3 << " "
-             << setprecision(3) << satData->L3 << endl;
-
         PP(iObs,iObs) = 1.0 / (sig_L3_GLO * sig_L3_GLO);
         for (int iPar = 1; iPar <= _params.size(); iPar++) {
@@ -648 +642 @@
         }
       
-        //// beg test
-        //double rhoCmp = cmpValue(satData);
-        //cout.setf(ios::fixed);
-        //cout << prn.toAscii().data() << " "  
-        //     << setprecision(3) << rhoCmp      << " "
-        //     << setprecision(3) << satData->P3 << " "
-        //     << setprecision(3) << satData->L3 << " " << endl;
-        //
-        ////// end test
-      }
-    }
-
-    cout << endl;
+        cout.setf(ios::fixed);
+        cout << "iObs = " << iObs << " " << prn.toAscii().data() << " " 
+             << setprecision(3) << rhoCmp << " "
+             << setprecision(3) << satData->P3 << " "
+             << setprecision(3) << satData->L3 << " "
+             << ll(iObs) << endl;
+      }
+    }
 
     // Compute Filter Update
@@ -672 +661 @@
     dx    = _QQ * ATP * ll; 
     vv    = ll - AA * dx;
-    cout << "vv = " << endl;
-    cout << vv.t() << endl;
+
+//    //// beg test
+//    {
+//      cout.setf(ios::fixed);
+//      ColumnVector vv_code(epoData->sizeGPS());
+//      ColumnVector vv_phase(epoData->sizeGPS());
+//      ColumnVector vv_glo(epoData->sizeGlo());
+//
+//      for (int iobs = 1; iobs <= epoData->sizeGPS(); ++iobs) {
+//        vv_code(iobs)  = vv(2*iobs-1);
+//        vv_phase(iobs) = vv(2*iobs);
+//      }
+//      for (int iobs = 1; iobs <= epoData->sizeGlo(); ++iobs) {
+//        vv_glo(iobs)  = vv(2*epoData->sizeGPS()+iobs);
+//      }
+//
+//      cout << "residuals code  " << setprecision(3) << vv_code.t(); 
+//      cout << "residuals phase " << setprecision(3) << vv_phase.t();
+//      cout << "residuals glo   " << setprecision(3) << vv_glo.t();
+//    }
+//    //// end test
 
   } while (outlierDetection(QQsav, vv, epoData->satDataGPS,