sgp4unit.cpp   sgp4unit.cpp 
skipping to change at line 1095 skipping to change at line 1095
2.0 * (d4410 * cos(x2omi + x2li - g44) + 2.0 * (d4410 * cos(x2omi + x2li - g44) +
d4422 * cos(x2li - g44) + d5421 * cos(xomi + x2li - g54) + d4422 * cos(x2li - g44) + d5421 * cos(xomi + x2li - g54) +
d5433 * cos(-xomi + x2li - g54)); d5433 * cos(-xomi + x2li - g54));
xnddt = xnddt * xldot; xnddt = xnddt * xldot;
} }
/* ----------------------- integrator ------------------- */ /* ----------------------- integrator ------------------- */
// sgp4fix move end checks to end of routine // sgp4fix move end checks to end of routine
if (fabs(t - atime) >= stepp) if (fabs(t - atime) >= stepp)
{ {
iret = 0; // NOTE: Never read!
// iret = 0;
iretn = 381; iretn = 381;
} }
else // exit here else // exit here
{ {
ft = t - atime; ft = t - atime;
iretn = 0; iretn = 0;
} }
if (iretn == 381) if (iretn == 381)
{ {
skipping to change at line 1706 skipping to change at line 1707
------------------------------------------------------------------------- ---*/ ------------------------------------------------------------------------- ---*/
bool sgp4 bool sgp4
( (
gravconsttype whichconst, elsetrec& satrec, double tsince, gravconsttype whichconst, elsetrec& satrec, double tsince,
double r[3], double v[3] double r[3], double v[3]
) )
{ {
double am , axnl , aynl , betal , cosim , cnod , double am , axnl , aynl , betal , cosim , cnod ,
cos2u, coseo1=0.0, cosi , cosip , cosisq, cossu , cosu, cos2u, coseo1=0.0, cosi , cosip , cosisq, cossu , cosu,
delm , delomg, em , emsq , ecose , el2 , eo1 , delm , delomg, em , ecose , el2 , eo1 ,
ep , esine , argpm, argpp , argpdf, pl, mrt = 0.0, ep , esine , argpm, argpp , argpdf, pl, mrt = 0.0,
mvt , rdotl , rl , rvdot , rvdotl, sinim , mvt , rdotl , rl , rvdot , rvdotl, sinim ,
sin2u, sineo1=0.0, sini , sinip , sinsu , sinu , sin2u, sineo1=0.0, sini , sinip , sinsu , sinu ,
snod , su , t2 , t3 , t4 , tem5 , temp, snod , su , t2 , t3 , t4 , tem5 , temp,
temp1, temp2 , tempa, tempe , templ , u , ux , temp1, temp2 , tempa, tempe , templ , u , ux ,
uy , uz , vx , vy , vz , inclm , mm , uy , uz , vx , vy , vz , inclm , mm ,
nm , nodem, xinc , xincp , xl , xlm , mp , nm , nodem, xinc , xincp , xl , xlm , mp ,
xmdf , xmx , xmy , nodedf, xnode , nodep, tc , dndt, xmdf , xmx , xmy , nodedf, xnode , nodep, tc , dndt,
twopi, x2o3 , j2=0 , j3 , tumin, j4=0 , xke=0 , j3oj2=0, r adiusearthkm=0, twopi, x2o3 , j2=0 , j3 , tumin, j4=0 , xke=0 , j3oj2=0, r adiusearthkm=0,
mu, vkmpersec; mu, vkmpersec;
skipping to change at line 1813 skipping to change at line 1814
// printf("# error em %f\n", em); // printf("# error em %f\n", em);
satrec.error = 1; satrec.error = 1;
// sgp4fix to return if there is an error in eccentricity // sgp4fix to return if there is an error in eccentricity
return false; return false;
} }
// sgp4fix fix tolerance to avoid a divide by zero // sgp4fix fix tolerance to avoid a divide by zero
if (em < 1.0e-6) if (em < 1.0e-6)
em = 1.0e-6; em = 1.0e-6;
mm = mm + satrec.no * templ; mm = mm + satrec.no * templ;
xlm = mm + argpm + nodem; xlm = mm + argpm + nodem;
emsq = em * em; // NOTE: Never read!
temp = 1.0 - emsq; // emsq = em * em;
// temp = 1.0 - emsq;
nodem = fmod(nodem, twopi); nodem = fmod(nodem, twopi);
argpm = fmod(argpm, twopi); argpm = fmod(argpm, twopi);
xlm = fmod(xlm, twopi); xlm = fmod(xlm, twopi);
mm = fmod(xlm - argpm - nodem, twopi); mm = fmod(xlm - argpm - nodem, twopi);
/* ----------------- compute extra mean quantities ------------- */ /* ----------------- compute extra mean quantities ------------- */
sinim = sin(inclm); sinim = sin(inclm);
cosim = cos(inclm); cosim = cos(inclm);
 End of changes. 3 change blocks. 
4 lines changed or deleted 6 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/