spacecruft
/
galmon
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

speed up ephemeris calculation by truncating iteration if we are done

hzcorr
bert hubert 2019-09-05 09:31:44 +02:00
parent 4de4d14800
commit dacba30428
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ephemeris.hh
View File

@ -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)