rust-roche 0.1.0

Rust translation of Tom Marsh's cpp-roche package for modelling Roche distorted stars/binary systems.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42

use std::f64::consts::PI;
use crate::Vec3;

///
/// set_earth_iangle computes the earth vector given an inclination and orbital phase. 
///
/// \param  iangle orbital inclination
/// \param  phase  orbital phase (1 = one orbit)
///
/// \return a Vec3 as the earth vector (unit vector)
///
pub fn set_earth_iangle(iangle: f64, phase: f64) -> Vec3 {
    let iangle_rad: f64 = iangle.to_radians();
    let phase_rad: f64 = 2.*PI*phase;
    let (sini, cosi) = iangle_rad.sin_cos();
    let (sinp, cosp) = phase_rad.sin_cos();
    Vec3{
        x: sini * cosp,
        y: -sini * sinp,
        z: cosi,
    }
}


/// 
/// set_earth computes the earth vector given an inclination and orbital phase. 
///
/// \param  cosi cosine of orbital inclination
/// \param  sini sine of orbital inclination
/// \param  phase  orbital phase (1 = one orbit)
///
/// \return a Vec3 as the earth vector (unit vector)
///
pub  fn set_earth(cosi: f64, sini: f64, phase: f64) -> Vec3 {
    let phase_rad: f64 = 2.*PI*phase;
    let (sinp, cosp) = phase_rad.sin_cos();
    Vec3{
        x: sini*cosp,
        y: -sini*sinp,
        z: cosi
    }
}