Struct BogackiShampine 

pub struct BogackiShampine {
    pub abs_tol: Scalar,
    pub rel_tol: Scalar,
    pub dt_beta: Scalar,
    pub dt_expn: Scalar,
    pub dt_cut: Scalar,
    pub dt_min: Scalar,
}

Explicit, three-stage, third-order, variable-step, Runge-Kutta method.¹

\frac{dy}{dt} = f(t, y)

t_{n+1} = t_n + h

k_1 = f(t_n, y_n)

k_2 = f(t_n + \tfrac{1}{2} h, y_n + \tfrac{1}{2} h k_1)

k_3 = f(t_n + \tfrac{3}{4} h, y_n + \tfrac{3}{4} h k_2)

y_{n+1} = y_n + \frac{h}{9}\left(2k_1 + 3k_2 + 4k_3\right)

k_4 = f(t_{n+1}, y_{n+1})

e_{n+1} = \frac{h}{72}\left(-5k_1 + 6k_2 + 8k_3 - 9k_4\right)

---

1. P. Bogacki and L.F. Shampine, Appl. Math. Lett. 2, 321 (1989). ↩

Fields

abs_tol: Scalar

Absolute error tolerance.

rel_tol: Scalar

Relative error tolerance.

dt_beta: Scalar

Multiplier for adaptive time steps.

dt_expn: Scalar

Exponent for adaptive time steps.

dt_cut: Scalar

Cut back factor for the time step.

dt_min: Scalar

Minimum value for the time step.

Trait Implementations

impl Debug for BogackiShampine

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Default for BogackiShampine

fn default() -> Self

Returns the “default value” for a type.

impl<Y, U> Explicit<Y, U> for BogackiShampine

where Self: OdeSolver<Y, U>, Y: Tensor, for<'a> &'a Y: Mul<Scalar, Output = Y> + Sub<&'a Y, Output = Y>, U: TensorVec<Item = Y>,

const SLOPES: usize = 4usize

fn dt_beta(&self) -> Scalar

Returns the multiplier for adaptive time steps.

fn dt_expn(&self) -> Scalar

Returns the exponent for adaptive time steps.

fn integrate( &self, function: impl FnMut(Scalar, &Y) -> Result<Y, String>, time: &[Scalar], initial_condition: Y, ) -> Result<(Vector, U, U), IntegrationError>

Solves an initial value problem by explicitly integrating a system of ordinary differential equations.

fn time_step(&self, error: Scalar, dt: &mut Scalar)

Provides the adaptive time step as a function of the error.

impl<Y, Z, U, V> ExplicitIV<Y, Z, U, V> for BogackiShampine

where Self: OdeSolver<Y, U>, Y: Tensor, Z: Tensor, for<'a> &'a Y: Mul<Scalar, Output = Y> + Sub<&'a Y, Output = Y>, U: TensorVec<Item = Y>, V: TensorVec<Item = Z>,

const SLOPES: usize = 4usize

fn integrate( &self, function: impl FnMut(Scalar, &Y, &Z) -> Result<Y, String>, evaluate: impl FnMut(Scalar, &Y, &Z) -> Result<Z, String>, time: &[Scalar], initial_condition: Y, initial_evaluation: Z, ) -> Result<(Vector, U, U, V), IntegrationError>

Solves an initial value problem by explicitly integrating a system of ordinary differential equations with internal variables.

impl<Y, U> InterpolateSolution<Y, U> for BogackiShampine

where Y: Tensor, for<'a> &'a Y: Mul<Scalar, Output = Y> + Sub<&'a Y, Output = Y>, U: TensorVec<Item = Y>,

fn interpolate( &self, time: &Vector, tp: &Vector, yp: &U, function: impl FnMut(Scalar, &Y) -> Result<Y, String>, ) -> Result<(U, U), IntegrationError>

Solution interpolation.

impl<Y, Z, U, V> InterpolateSolutionIV<Y, Z, U, V> for BogackiShampine

where Y: Tensor, Z: Tensor, for<'a> &'a Y: Mul<Scalar, Output = Y> + Sub<&'a Y, Output = Y>, U: TensorVec<Item = Y>, V: TensorVec<Item = Z>,

fn interpolate( &self, time: &Vector, tp: &Vector, yp: &U, zp: &V, function: impl FnMut(Scalar, &Y, &Z) -> Result<Y, String>, evaluate: impl FnMut(Scalar, &Y, &Z) -> Result<Z, String>, ) -> Result<(U, U, V), IntegrationError>

Solution interpolation with internal variables.

impl<Y, U> OdeSolver<Y, U> for BogackiShampine

where Y: Tensor, U: TensorVec<Item = Y>,

fn abs_tol(&self) -> Scalar

Returns the absolute error tolerance.

fn dt_cut(&self) -> Scalar

Returns the cut back factor for function errors.

fn dt_min(&self) -> Scalar

Returns the minimum value for the time step.