nyx_space/propagators/rk_methods/
verner.rs

1/*
2    Nyx, blazing fast astrodynamics
3    Copyright (C) 2018-onwards Christopher Rabotin <christopher.rabotin@gmail.com>
4
5    This program is free software: you can redistribute it and/or modify
6    it under the terms of the GNU Affero General Public License as published
7    by the Free Software Foundation, either version 3 of the License, or
8    (at your option) any later version.
9
10    This program is distributed in the hope that it will be useful,
11    but WITHOUT ANY WARRANTY; without even the implied warranty of
12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13    GNU Affero General Public License for more details.
14
15    You should have received a copy of the GNU Affero General Public License
16    along with this program.  If not, see <https://www.gnu.org/licenses/>.
17*/
18
19use super::RK;
20
21/// `Verner56` is an RK Verner integrator of order 5-6.
22///
23/// Coefficients taken from [here (PDF)](http://people.math.sfu.ca/~jverner/classify.1992.ps).
24pub(crate) struct Verner56 {}
25
26impl RK for Verner56 {
27    const ORDER: u8 = 6;
28    const STAGES: usize = 8;
29
30    const A_COEFFS: &'static [f64] = &[
31        1.0 / 6.0,
32        4.0 / 75.0,
33        16.0 / 75.0,
34        5.0 / 6.0,
35        -8.0 / 3.0,
36        5.0 / 2.0,
37        -165.0 / 64.0,
38        55.0 / 6.0,
39        -425.0 / 64.0,
40        85.0 / 96.0,
41        -8_263.0 / 15_000.0,
42        124.0 / 75.0,
43        -643.0 / 680.0,
44        -81.0 / 250.0,
45        2_484.0 / 10_625.0,
46        3_501.0 / 1_720.0,
47        -300.0 / 43.0,
48        297_275.0 / 52_632.0,
49        -319.0 / 2_322.0,
50        24_068.0 / 84_065.0,
51        3_850.0 / 26_703.0,
52        12.0 / 5.0,
53        -8.0,
54        4_015.0 / 612.0,
55        -11.0 / 36.0,
56        88.0 / 255.0,
57        0.0,
58        0.0,
59    ];
60
61    const B_COEFFS: &'static [f64] = &[
62        3.0 / 40.0,
63        0.0,
64        875.0 / 2_244.0,
65        23.0 / 72.0,
66        264.0 / 1_955.0,
67        125.0 / 11_592.0,
68        43.0 / 616.0,
69        0.0,
70        13.0 / 160.0,
71        0.0,
72        2_375.0 / 5_984.0,
73        5.0 / 16.0,
74        12.0 / 85.0,
75        0.0,
76        0.0,
77        3.0 / 44.0,
78    ];
79}
nyx_space/propagators/rk_methods/verner.rs

nyx_space/propagators/rk_methods/
verner.rs
