Struct ProjectionMethod

Source

pub struct ProjectionMethod<F: Float + FromPrimitive + Debug + ScalarOperand + AddAssign + SubAssign + MulAssign + DivAssign + Display + LowerExp + Sum + IntegrateFloat> {
    pub structure: DAEStructure<F>,
    pub project_after_step: bool,
    pub consistent_initialization: bool,
    pub constraint_tol: F,
}

Expand description

Projection method for index reduction and constraint satisfaction

This method combines index reduction with projection-based stabilization to maintain constraint satisfaction during integration.

Fields§

§structure: DAEStructure<F>

DAE structure information

§project_after_step: bool

Whether to project after each step

§consistent_initialization: bool

Whether to use consistent initialization

§constraint_tol: F

Tolerance for constraint satisfaction

Implementations§

Source §

impl<F: Float + FromPrimitive + Debug + ScalarOperand + AddAssign + SubAssign + MulAssign + DivAssign + Display + LowerExp + Sum + IntegrateFloat> ProjectionMethod<F>

Source

pub fn new(structure: DAEStructure<F>) -> Self

Create a new projection method for constraint satisfaction

Source

pub fn project_solution<GFunc>( &self, t: F, x: &mut Array1<F>, y: &mut Array1<F>, g: &GFunc, ) -> IntegrateResult<()>
where GFunc: Fn(F, ArrayView1<'_, F>, ArrayView1<'_, F>) -> Array1<F>,

Project the solution onto the constraint manifold

This ensures that the solution satisfies the constraints exactly.

Source

pub fn make_consistent<GFunc>( &self, t: F, x: &mut Array1<F>, y: &mut Array1<F>, g: &GFunc, ) -> IntegrateResult<()>
where GFunc: Fn(F, ArrayView1<'_, F>, ArrayView1<'_, F>) -> Array1<F>,

Make the initial conditions consistent with the constraints

This projects the initial state onto the constraint manifold.

Auto Trait Implementations§

§

impl<F> UnwindSafe for ProjectionMethod<F>
where F: UnwindSafe + RefUnwindSafe,

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 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> IntoEither for T

Source §

fn into_either(self, into_left: bool) -> Either<Self, Self>

Converts self into a Left variant of Either<Self, Self> if into_left is true. Converts self into a Right variant of Either<Self, Self> otherwise. Read more

Source §

fn into_either_with<F>(self, into_left: F) -> Either<Self, Self>
where F: FnOnce(&Self) -> bool,

Converts self into a Left variant of Either<Self, Self> if into_left(&self) returns true. Converts self into a Right variant of Either<Self, Self> otherwise. Read more

Source §