diff --git a/src/orbits.cpp b/src/orbits.cpp index e50a0e8..da12546 100644 --- a/src/orbits.cpp +++ b/src/orbits.cpp @@ -72,14 +72,12 @@ orbitGetEccentricityVector(glm::dvec3 r, glm::dvec3 v, double mu) return 1 / mu * ((v_mag - mu / r_mag) * r - (glm::dot(r, v) * v)); } -// NOTE: returns position vector in perifocal plane glm::dvec3 orbitGetPositionVector(double r, double theta) { return glm::dvec3(r * cos(theta), r * sin(theta), 0); } -// NOTE: returns velocity vector in perifocal plane glm::dvec3 orbitGetVelocityVector(double mu, double h, double e, double theta) { diff --git a/src/orbits.h b/src/orbits.h index fca8677..156e70e 100644 --- a/src/orbits.h +++ b/src/orbits.h @@ -110,8 +110,10 @@ orbitInit(double a, double e, double iota, double ohm, double omega, double nu); glm::dvec3 orbitGetEccentricityVector(glm::dvec3 r, glm::dvec3 v, double mu); +// NOTE: returns position vector in perifocal plane glm::dvec3 orbitGetPositionVector(double r, double theta); +// NOTE: returns velocity vector in perifocal plane glm::dvec3 orbitGetVelocityVector(double mu, double h, double e, double theta); // NOTE: return transform from perifocal to grav body (IJK) coordinate system