sgp4ext.cpp   sgp4ext.cpp 
skipping to change at line 74 skipping to change at line 74
* *
* coupling : * coupling :
* none. * none.
* ------------------------------------------------------------------------- -- */ * ------------------------------------------------------------------------- -- */
double mag double mag
( (
double x[3] double x[3]
) )
{ {
return sqrt(x[0]*x[0] + x[1]*x[1] + x[2]*x[2]); return std::sqrt(x[0]*x[0] + x[1]*x[1] + x[2]*x[2]);
} // end mag } // end mag
/* ------------------------------------------------------------------------ ----- /* ------------------------------------------------------------------------ -----
* *
* procedure cross * procedure cross
* *
* this procedure crosses two vectors. * this procedure crosses two vectors.
* *
* author : david vallado 719-573-2600 1 mar 200 1 * author : david vallado 719-573-2600 1 mar 200 1
* *
skipping to change at line 246 skipping to change at line 246
// --------------------------- circular ------------------------ // --------------------------- circular ------------------------
if ( fabs( ecc ) < sv ) if ( fabs( ecc ) < sv )
{ {
m = nu; m = nu;
e0= nu; e0= nu;
} }
else else
// ---------------------- elliptical ----------------------- // ---------------------- elliptical -----------------------
if ( ecc < 1.0-sv ) if ( ecc < 1.0-sv )
{ {
sine= ( sqrt( 1.0 -ecc*ecc ) * sin(nu) ) / ( 1.0 +ecc*cos(nu) ); sine= ( std::sqrt( 1.0 -ecc*ecc ) * sin(nu) ) / ( 1.0 +ecc*cos( nu) );
cose= ( ecc + cos(nu) ) / ( 1.0 + ecc*cos(nu) ); cose= ( ecc + cos(nu) ) / ( 1.0 + ecc*cos(nu) );
e0 = atan2( sine,cose ); e0 = atan2( sine,cose );
m = e0 - ecc*sin(e0); m = e0 - ecc*sin(e0);
} }
else else
// -------------------- hyperbolic -------------------- // -------------------- hyperbolic --------------------
if ( ecc > 1.0 + sv ) if ( ecc > 1.0 + sv )
{ {
if ((ecc > 1.0 ) && (fabs(nu)+0.00001 < M_PI-acos(1.0 /ecc ))) if ((ecc > 1.0 ) && (fabs(nu)+0.00001 < M_PI-acos(1.0 /ecc )))
{ {
sine= ( sqrt( ecc*ecc-1.0 ) * sin(nu) ) / ( 1.0 + ec c*cos(nu) ); sine= ( std::sqrt( ecc*ecc-1.0 ) * sin(nu) ) / ( 1.0 + ecc*cos(nu) );
#ifdef _MSC_BUILD #ifdef _MSC_BUILD
e0 = StelUtils::asinh( sine ); e0 = StelUtils::asinh( sine );
#else #else
e0 = asinh( sine ); e0 = asinh( sine );
#endif #endif
m = ecc*sinh(e0) - e0; m = ecc*sinh(e0) - e0;
} }
} }
else else
// ----------------- parabolic --------------------- // ----------------- parabolic ---------------------
 End of changes. 3 change blocks. 
3 lines changed or deleted 3 lines changed or added

This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/