scialg 0.3.0

A collection of scientific algorithms
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

pub mod dormand_prince;
pub mod euler;
pub mod midpoint;
pub mod runge_kutta;

use crate::vector::Vector;

pub enum StepperMethod {
    Euler,
    Midpoint,
    RungeKutta,
    DormandPrince,
}

struct StepperData<const N: usize> {
    h_cur: f64,
    x_cur: f64,
    y_cur: Vector<N>,
    derive: fn(f64, Vector<N>) -> Vector<N>,
}

impl<const N: usize> StepperData<N> {
    fn new(h: f64, p0: Vector<N>, df: fn(f64, Vector<N>) -> Vector<N>) -> Self {
        StepperData {
            h_cur: h,
            x_cur: 0.0,
            y_cur: p0,
            derive: df,
        }
    }
}

pub trait Stepper<const N: usize> {
    fn step(&mut self) -> Vector<N>;
}