fly/flightgear
1
0
Fork 0
Code Activity

math domain error fix from Charlie Hotchkiss.

Browse source
This commit is contained in:
curt 1998-12-04 01:25:05 +00:00
parent cfc9e51128
commit 01d6408def
1 changed files with 14 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

16
Astro/celestialBody.cxx
View file

@ -27,7 +27,7 @@
#include "star.hxx" #include "star.hxx"
#include <Debug/logstream.hxx> #include <Debug/logstream.hxx>
#ifdef __BORLANDC__ #ifdef FG_MATH_EXCEPTION_CLASH
# define exception c_exception # define exception c_exception
#endif #endif
#include <math.h> #include <math.h>
@ -90,7 +90,19 @@ void CelestialBody::updatePosition(fgTIME *t, Star *ourSun)
//of the planet //of the planet
R = sqrt (xg*xg + yg*yg + zg*zg); R = sqrt (xg*xg + yg*yg + zg*zg);
s = ourSun->getDistance(); s = ourSun->getDistance();
FV = RAD_TO_DEG * acos( (r*r + R*R - s*s) / (2*r*R));
// It is possible from these calculations for the argument to acos
// to exceed the valid range for acos(). So we do a little extra
// checking.
double tmp = (r*r + R*R - s*s) / (2*r*R);
if ( tmp > 1.0) {
tmp = 1.0;
} else if ( tmp < -1.0) {
tmp = -1.0;
}
FV = RAD_TO_DEG * acos( tmp );
}; };
/**************************************************************************** /****************************************************************************