← Back to branch summary

~georg-zotti/stellarium/gz_planetAxes

« back to all changes in this revision

Viewing changes to src/core/modules/Orbit.cpp

  • Committer: Georg Zotti
  • Date: 2017-01-06 00:15:43 UTC
  • Revision ID: georg.zotti@univie.ac.at-20170106001543-vbpmr4ocoidaunb8
New rotational elements in use. Actual rotation (siderealTime) not properly working yet. A few other codefixes.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/core/modules/Orbit.cpp
99
99
//      const double M=dt*sqrt(GAUSS_GRAV_CONST/(2.0*q*q*q));
100
100
//      const double W=1.5*M;
101
101
        const double W=dt*n;
102
 
        const double Y=cbrt(W+std::sqrt(W*W+1));
 
102
        const double Y=cbrt(W+std::sqrt(W*W+1.));
103
103
        const double tanNu2=Y-1.0/Y; // Heafner (5.5.8) has an error here, writes (Y-1)/Y.
104
104
        rCosNu=q*(1.0-tanNu2*tanNu2);
105
105
        rSinNu=2.0*q*tanNu2;
229
229
          orbitGood(orbitGoodDays)
230
230
{
231
231
        // GZ MAKE SURE THIS IS ALWAYS 0/0/0.
232
 
        qDebug() << "parentRotObliquity" << parentRotObliquity << "parentRotAscendingnode" << parentRotAscendingnode << "parentRotJ2000Longitude" << parentRotJ2000Longitude;
 
232
        //qDebug() << "parentRotObliquity" << parentRotObliquity << "parentRotAscendingnode" << parentRotAscendingnode << "parentRotJ2000Longitude" << parentRotJ2000Longitude;
233
233
        rdot.set(0.0, 0.0, 0.0);
234
234
        const double c_obl = cos(parentRotObliquity);
235
235
        const double s_obl = sin(parentRotObliquity);