Struct peroxide::numerical::ode::BasicODESolver

pub struct BasicODESolver<I: ODEIntegrator> { /* private fields */ }

A basic ODE solver using a specified integrator.

Example

use peroxide::fuga::*;

fn main() -> Result<(), Box<dyn Error>> {
    let rkf = RKF45::new(1e-4, 0.9, 1e-6, 1e-1, 100);
    let basic_ode_solver = BasicODESolver::new(rkf);
    let (t_vec, y_vec) = basic_ode_solver.solve(
        &Test,
        (0f64, 10f64),
        0.01,
    )?;
    let y_vec: Vec<f64> = y_vec.into_iter().flatten().collect();

    Ok(())
}

struct Test;

impl ODEProblem for Test {
    fn initial_conditions(&self) -> Vec<f64> {
        vec![1f64]
    }

    fn rhs(&self, t: f64, y: &[f64], dy: &mut [f64]) -> anyhow::Result<()> {
        dy[0] = (5f64 * t.powi(2) - y[0]) / (t + y[0]).exp();
        Ok(())
    }
}

Implementations

impl<I: ODEIntegrator> BasicODESolver<I>

pub fn new(integrator: I) -> Self

Trait Implementations

impl<I: ODEIntegrator> ODESolver for BasicODESolver<I>

fn solve<P: ODEProblem>( &self, problem: &P, t_span: (f64, f64), dt: f64, ) -> Result<(Vec<f64>, Vec<Vec<f64>>)>