spruce up ephemeris calculation debugging output a bit

berthubert · Sep 11, 2022 · 05a6c6f · 05a6c6f
1 parent 4546811
commit 05a6c6f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/ephemeris.hh b/ephemeris.hh
@@ -89,7 +89,7 @@ double getCoordinates(double tow, const T& iod, Point* p, bool quiet=true)
 
     cerr << "sqrtA = "<< sqrtA << endl;
     cerr << "deltan = "<< deltan << endl;
-    cerr << "t0e = "<< t0e << endl;
+    cerr << "t0e = "<< t0e << "( rel "<<(tow - t0e)<<")"<<endl;
     cerr << "m0 = "<< m0 << " ("<<todeg(m0)<<")"<<endl;
     cerr << "e = "<< e << endl;
     cerr << "omega = " << omega << " ("<<todeg(omega)<<")"<<endl;
@@ -122,7 +122,7 @@ double getCoordinates(double tow, const T& iod, Point* p, bool quiet=true)
 
   double M = m0 + n * tk;
   if(!quiet)
-    cerr << " M = m0 + n * tk = "<<m0 << " + " << n << " * " << tk <<endl;
+    cerr << " M = m0 + n * tk = "<<m0 << " + " << n << " * " << tk << " = " <<M <<endl;
   double E = M;
   double newE;
   for(int k =0 ; k < 10; ++k) {