PointMass.hpp

Go to the documentation of this file.

// Vulcan Point Mass Gravity Model
// Simple inverse-square gravitational acceleration
#pragma once
 
#include <janus/janus.hpp>
#include <vulcan/core/Constants.hpp>
#include <vulcan/core/VulcanTypes.hpp>
 
namespace vulcan::gravity::point_mass {

template <typename Scalar>
Vec3<Scalar> acceleration(const Vec3<Scalar> &r_ecef,
                          double mu = constants::earth::mu) {
    const Scalar r_mag = janus::norm(r_ecef);
    const Scalar r_cubed = r_mag * r_mag * r_mag;
 
    // g = -μ/r³ · r
    return -mu / r_cubed * r_ecef;
}

template <typename Scalar>
Scalar potential(const Vec3<Scalar> &r_ecef, double mu = constants::earth::mu) {
    const Scalar r_mag = janus::norm(r_ecef);
    return -mu / r_mag;
}

template <typename Scalar>
Scalar acceleration_magnitude(const Scalar &r_mag,
                              double mu = constants::earth::mu) {
    return mu / (r_mag * r_mag);
}
 
} // namespace vulcan::gravity::point_mass