Skip to main content

StiffSolver

Struct StiffSolver 

Source
pub struct StiffSolver {
    pub h: f64,
    pub tol: f64,
}
Expand description

Implicit solver for stiff ODEs using backward Euler (BDF-1) with Newton iteration.

Fields§

§h: f64

Step size h

§tol: f64

Newton iteration tolerance

Implementations§

Source§

impl StiffSolver

Source

pub fn new(h: f64, tol: f64) -> Self

Create with given step size and tolerance.

Source

pub fn step(&self, f: &dyn Fn(f64, f64) -> f64, t: f64, y: f64) -> f64

Perform one backward Euler step via fixed-point iteration.

Source

pub fn solve_to( &self, f: &dyn Fn(f64, f64) -> f64, t0: f64, y0: f64, t_end: f64, ) -> Vec<(f64, f64)>

Solve from t0 to t_end using backward Euler.

Auto Trait Implementations§

Blanket Implementations§

Source§

impl<T> Any for T
where T: 'static + ?Sized,

Source§

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more
Source§

impl<T> Borrow<T> for T
where T: ?Sized,

Source§

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more
Source§

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source§

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more
Source§

impl<T> From<T> for T

Source§

fn from(t: T) -> T

Returns the argument unchanged.

Source§

impl<T, U> Into<U> for T
where U: From<T>,

Source§

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source§

impl<T, U> TryFrom<U> for T
where U: Into<T>,

Source§

type Error = Infallible

The type returned in the event of a conversion error.
Source§

fn try_from(value: U) -> Result<T, <T as TryFrom<U>>::Error>

Performs the conversion.
Source§

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

Source§

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.
Source§

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.