KJ135
/
sgp4
mirror of http://172.16.0.12:5000/LAMPSARToolSoftware/sgp4.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Start of work on DeepSecular()

feature/19
Daniel Warner 2011-03-30 01:17:05 +01:00
parent a8b25f89c5
commit 909a02b2db
1 changed files with 12 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
SGDP4.cpp
View File

@ -1006,35 +1006,25 @@ void SGDP4::DeepPeriodics(const double& t, double& em,
*/ */
void SGDP4::DeepSecular(const double& t, double& xll, double& omgasm, void SGDP4::DeepSecular(const double& t, double& xll, double& omgasm,
double& xnodes, double& em, double& xinc, double& xn) { double& xnodes, double& em, double& xinc, double& xn) {
#if 0
static const double step = 720.0;
double xldot = 0.0; double xldot = 0.0;
double xndot = 0.0; double xndot = 0.0;
double xnddt = 0.0; double xnddt = 0.0;
xll = xll + d_ssl_ * t; xll += d_ssl_ * t;
omgasm = omgasm + d_ssg_ * t; omgasm += d_ssg_ * t;
xnodes = xnodes + d_ssh_ * t; xnodes += d_ssh_ * t;
em = Eccentricity() + d_sse_ * t; em += d_sse_ * t;
xinc = Inclination() + d_ssi_ * t; xinc += d_ssi_ * t;
double stepp = 720.0; if (!d_resonance_flag_)
double stepn = -720.0;
double step2 = 259200.0;
/*
* check if needed
*/
//if (xinc >= 0.0) {
// xinc = -xinc;
// xnodes = xnodes + Globals::PI();
// omgasm = omgasm - Globals::PI();
//}
if (iresfl == 0)
return; return;
#if 0
if (fabs(d_atime_) < d_stepp_ || if (fabs(d_atime_) < d_stepp_ ||
(d_atime_ < 0.0 && t < d_atime_ - 1.0) || (d_atime_ < 0.0 && t < d_atime_ - 1.0) ||
(d_atime_ > 0.0 && t > d_atime_ + 1.0)) { (d_atime_ > 0.0 && t > d_atime_ + 1.0)) {
@ -1105,6 +1095,7 @@ void SGDP4::DeepSecular(const double& t, double& xll, double& omgasm,
ft = t - d_atime_; ft = t - d_atime_;
} while (fabs(ft) >= d_stepp_); } while (fabs(ft) >= d_stepp_);
} }
#endif
/* /*
* integrator * integrator
*/ */
@ -1113,9 +1104,8 @@ void SGDP4::DeepSecular(const double& t, double& xll, double& omgasm,
double xl = d_xli_ + xldot * ft + xndot * ft * ft * 0.5; double xl = d_xli_ + xldot * ft + xndot * ft * ft * 0.5;
double temp = -xnodes + thgr + t * thdt; double temp = -xnodes + thgr + t * thdt;
if (isynfl == 0) if (!d_synchronous_flag_)
xll = xl + temp + temp; xll = xl + temp + temp;
else else
xll = xl - omgasm + temp; xll = xl - omgasm + temp;
#endif
} }