rebound-rs 4.6.0-alpha.1

Rust wrapper for the REBOUND N-body simulation library.
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95

use rebound_bind as rb;

#[derive(Debug, Clone, Copy, PartialEq, Eq)]
pub enum AdaptiveMode {
    Individual = rb::reb_integrator_ias15_REB_IAS15_INDIVIDUAL as isize,
    Global = rb::reb_integrator_ias15_REB_IAS15_GLOBAL as isize,
    Prs23 = rb::reb_integrator_ias15_REB_IAS15_PRS23 as isize,
    Aarseth85 = rb::reb_integrator_ias15_REB_IAS15_AARSETH85 as isize,
}

impl From<AdaptiveMode> for rb::reb_integrator_ias15__bindgen_ty_1 {
    fn from(value: AdaptiveMode) -> Self {
        match value {
            AdaptiveMode::Individual => rb::reb_integrator_ias15_REB_IAS15_INDIVIDUAL,
            AdaptiveMode::Global => rb::reb_integrator_ias15_REB_IAS15_GLOBAL,
            AdaptiveMode::Prs23 => rb::reb_integrator_ias15_REB_IAS15_PRS23,
            AdaptiveMode::Aarseth85 => rb::reb_integrator_ias15_REB_IAS15_AARSETH85,
        }
    }
}

impl AdaptiveMode {
    fn from_raw(value: rb::reb_integrator_ias15__bindgen_ty_1) -> Option<Self> {
        match value {
            rb::reb_integrator_ias15_REB_IAS15_INDIVIDUAL => Some(Self::Individual),
            rb::reb_integrator_ias15_REB_IAS15_GLOBAL => Some(Self::Global),
            rb::reb_integrator_ias15_REB_IAS15_PRS23 => Some(Self::Prs23),
            rb::reb_integrator_ias15_REB_IAS15_AARSETH85 => Some(Self::Aarseth85),
            _ => None,
        }
    }
}

pub struct IntegratorIas15<'a> {
    pub(crate) sim: *mut rb::reb_simulation,
    pub(crate) inner: *mut rb::reb_integrator_ias15,
    pub(crate) _marker: core::marker::PhantomData<&'a mut rb::reb_simulation>,
}

impl<'a> IntegratorIas15<'a> {
    pub fn set_epsilon(self, epsilon: f64) -> Self {
        unsafe {
            (*self.inner).epsilon = epsilon;
        }
        self
    }

    pub fn set_min_dt(self, min_dt: f64) -> Self {
        unsafe {
            (*self.inner).min_dt = min_dt;
        }
        self
    }

    pub fn set_adaptive_mode(self, adaptive_mode: AdaptiveMode) -> Self {
        unsafe {
            (*self.inner).adaptive_mode = adaptive_mode.into();
        }
        self
    }

    pub fn epsilon(&self) -> f64 {
        unsafe { (*self.inner).epsilon }
    }

    pub fn min_dt(&self) -> f64 {
        unsafe { (*self.inner).min_dt }
    }

    pub fn adaptive_mode(&self) -> Option<AdaptiveMode> {
        unsafe { AdaptiveMode::from_raw((*self.inner).adaptive_mode) }
    }

    pub fn iterations_max_exceeded(&self) -> u64 {
        unsafe { (*self.inner).iterations_max_exceeded }
    }

    pub fn timescale(&mut self) -> f64 {
        unsafe { rb::reb_integrator_ias15_timescale(self.sim) }
    }

    pub fn step(&mut self) -> &mut Self {
        unsafe {
            rb::reb_integrator_ias15_step(self.sim);
        }
        self
    }

    pub fn reset(&mut self) -> &mut Self {
        unsafe {
            rb::reb_integrator_ias15_reset(self.sim);
        }
        self
    }
}