Struct DoPri5 

pub struct DoPri5;

Dormand-Prince 5(4) solver.

Example

use numra_ode::{OdeProblem, DoPri5, Solver, SolverOptions};

// Exponential decay: dy/dt = -y
let problem = OdeProblem::new(
    |_t, y: &[f64], dydt: &mut [f64]| { dydt[0] = -y[0]; },
    0.0,
    5.0,
    vec![1.0]
);

let options = SolverOptions::default();
let result = DoPri5::solve(&problem, 0.0, 5.0, &[1.0], &options).unwrap();

// y(5) ≈ e^(-5) ≈ 0.00674
let y_final = result.y_final().unwrap();
assert!((y_final[0] - (-5.0_f64).exp()).abs() < 1e-5);

Implementations

impl DoPri5

pub fn new() -> Self

Create a new DoPri5 solver.

Trait Implementations

impl Clone for DoPri5

fn clone(&self) -> DoPri5

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for DoPri5

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

Formats the value using the given formatter.

impl Default for DoPri5

fn default() -> DoPri5

Returns the “default value” for a type.

impl<S: Scalar> Solver<S> for DoPri5

fn solve<Sys: OdeSystem<S>>( problem: &Sys, t0: S, tf: S, y0: &[S], options: &SolverOptions<S>, ) -> Result<SolverResult<S>, SolverError>

Solve the ODE problem.