TwoBody.h

Go to the documentation of this file.

#pragma once
 
#include "objects/geometry/Vector.h"
#include "objects/wrappers/Optional.h"
 
NAMESPACE_SPH_BEGIN
 
namespace Kepler {
 
struct Elements {
    Float a; 
    Float e; 
    Float i; 
 
    Vector L; 
    Vector K; 
 
 
    Float periapsisArgument() const;
 
    Float ascendingNode() const;
 
    Float pericenterDist() const;
 
    Float semiminorAxis() const;
};
 
Optional<Elements> computeOrbitalElements(const Float M, const Float mu, const Vector& r, const Vector& v);
 
Float solveKeplersEquation(const Float M, const Float e, const Size iterCnt = 10);
 
Float eccentricAnomalyToTrueAnomaly(const Float u, const Float e);
 
Float trueAnomalyToEccentricAnomaly(const Float v, const Float e);
 
Vector position(const Float a, const Float e, const Float u);
 
Vector velocity(const Float a, const Float e, const Float u, const Float n);
 
Float meanMotion(const Float a, const Float m_total);
 
} // namespace Kepler
 
NAMESPACE_SPH_END