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

due to rounding errors arccos can sometimes get fed 1.0000000002, which is bad.

experimental
bert hubert 2019-10-23 10:32:03 +02:00
parent de0faa94b7
commit 4786097822
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
ephemeris.cc
View File

@ -42,7 +42,10 @@ std::pair<double,double> getLongLat(double x, double y, double z)
longitude *= -1;
Vector toUs{core, pos};
double latitude = acos( flat.inner(toUs) / (toUs.length() * flat.length()));
double inp = flat.inner(toUs) / (toUs.length() * flat.length());
if(inp > 1.0 && inp < 1.0000001) // this happens because of rounding errors
inp=1.0;
double latitude = acos( inp);
if(z < 0)
latitude *= -1;