Changeset 899 in ntrip

Timestamp:

May 9, 2008, 7:53:11 AM (16 years ago)

Author:

mervart

Message:

* empty log message *

File:

• trunk/BNS/bnseph.cpp (modified) (2 diffs)

trunk/BNS/bnseph.cpp

@@ -349 +349 @@
                         ColumnVector& vv) const {
 
-  const static double secPerWeek = 7 * 86400.0;
+  const static double secPerWeek  = 7 * 86400.0;
+  const static double nominalStep = 10.0;
 
   double dtPos = GPSweeks - _tt;
@@ -356 +357 @@
   }
 
+  int nSteps  = int(fabs(dtPos) / nominalStep) + 1;
+  double step = dtPos / nSteps;
+
   cout << _prn.toAscii().data() << "   " 
        << GPSweek  << " " << _GPSweek << "   "
        << GPSweeks << " " <<  _GPSweeks << " " << _tt << "    "
-       << dtPos << endl;
-
-  const static double maxStep = 10.0;
-
-  double tt = 0.0;
-  while (fabs(tt) < fabs(dtPos)) {
-    double step = dtPos > 0.0 ?  maxStep : -maxStep;
-    if (fabs(tt + step) > fabs(dtPos)) {
-      step = dtPos - tt;
-    }
-    _xv = rungeKutta4(tt, _xv, step, glo_deriv);
-    tt += step;
-  }
-
-  // Next Time Start Integration from Current Epoch
-  // ----------------------------------------------
-  _tt += dtPos;
+       << dtPos << " " << nSteps << " " << " " << step << endl;
+
+  for (int ii = 1; ii <= nSteps; ii++) { 
+    _tt += step;
+    _xv = rungeKutta4(_tt, _xv, step, glo_deriv);
+  }
 
   // Position and Velocity