From dacba304289f5e96d96e0b61f83c930a5d98a893 Mon Sep 17 00:00:00 2001
From: bert hubert <bert.hubert@netherlabs.nl>
Date: Thu, 5 Sep 2019 09:31:44 +0200
Subject: [PATCH] speed up ephemeris calculation by truncating iteration if we
 are done

---
 ephemeris.hh | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/ephemeris.hh b/ephemeris.hh
index 07198b4..39bef97 100644
--- a/ephemeris.hh
+++ b/ephemeris.hh
@@ -79,10 +79,16 @@ void getCoordinates(int wn, double tow, const T& iod, Point* p, bool quiet=true)
   if(!quiet)
     cerr << " M = m0 + n * tk = "<<m0 << " + " << n << " * " << tk <<endl;
   double E = M;
+  double newE;
   for(int k =0 ; k < 10; ++k) {
     if(!quiet)
-      cerr<<"M = "<<M<<", E = "<< E << endl;
-    E = M + e * sin(E);
+      cerr<<"k "<<k<<" M = "<<M<<", E = "<< E << endl;
+    newE = M + e * sin(E);
+    if(fabs(E-newE) < 0.00001) {
+      E = newE;
+      break;
+    }
+    E = newE;
   }
 
   // M = E  - e * sin(E)   -> E = M + e * sin(E)