Skip to content

Commit

Permalink
stuff
Browse files Browse the repository at this point in the history
  • Loading branch information
Dzejkop committed Jan 13, 2024
1 parent 5ebca96 commit bf9979a
Show file tree
Hide file tree
Showing 3 changed files with 7 additions and 9 deletions.
4 changes: 1 addition & 3 deletions examples/orbit_viewer.rs
Original file line number Diff line number Diff line change
Expand Up @@ -135,9 +135,7 @@ fn ui(
let h =
orbit.specific_angular_momentum(state.star_mass);
let μ =
standard_gravitational_parameter(
state.star_mass,
);
standard_gravitational_parameter(state.star_mass);

let r = (h.powi(2) / μ) / (1.0 + e * v.cos());
ui.label(format!("r: {r}"));
Expand Down
6 changes: 5 additions & 1 deletion src/astro.rs
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
use crate::constants::G;
use crate::constants::{G, TWO_PI};
use crate::Num;

pub mod elliptic;
Expand All @@ -13,3 +13,7 @@ pub fn standard_gravitational_parameter(mass: Num) -> Num {
pub fn soi(r: Num, m1: Num, m2: Num) -> Num {
r * (m1 / m2).powf(2.0 / 5.0)
}

pub fn period(a: Num, mass: Num) -> Num {
TWO_PI * (a.powi(3) / standard_gravitational_parameter(mass)).sqrt()
}
6 changes: 1 addition & 5 deletions src/elements.rs
Original file line number Diff line number Diff line change
Expand Up @@ -60,11 +60,7 @@ impl KeplerianElements {

/// https://en.wikipedia.org/wiki/Orbital_period
pub fn period(&self, mass: Num) -> Num {
Self::period_static(self.semi_major_axis, mass)
}

pub fn period_static(a: Num, mass: Num) -> Num {
TWO_PI * (a.powi(3) / standard_gravitational_parameter(mass)).sqrt()
astro::period(self.semi_major_axis, mass)
}

/// https://en.wikipedia.org/wiki/Mean_anomaly
Expand Down

0 comments on commit bf9979a

Please sign in to comment.