Struct InternalSolve 

pub(super) struct InternalSolve {
    pub(super) value: DMatrix<Complex>,
    factor: InternalFactor,
}

The forward solution of the internal system and what its adjoint solve reuses.

Fields

value: DMatrix<Complex>

Solution of the internal system for the forward right-hand sides.

factor: InternalFactor

Implementations

impl InternalSolve

pub(in smatrix) fn state_size( n: usize, columns: usize, ) -> Result<usize>

impl InternalSolve

pub(super) fn new( lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, ) -> Result<Self>

pub(super) fn direct( lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, ) -> Result<Self>

The pivoted LU solve, the branch Self::new takes whenever the Krylov conditions fail.

fn with_threshold( lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, min_rows: usize, ) -> Result<Self>

Self::new with a lower row threshold, which tests use to reach every dispatch branch with small systems.

fn lu(operator: DMatrix<Complex>, rhs: DMatrix<Complex>) -> Result<Self>

Factor the operator and solve for rhs in place.

pub(super) fn solve_adjoint( &self, lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, ) -> Result<AdjointSolve<'_>>

Solve the adjoint system for rhs, reusing the forward solve. lower and upper must be the reflection pair used by that forward solve.

pub(super) fn solve_forward( &self, lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, ) -> Result<DMatrix<Complex>>

Solve a tangent right-hand side with the recorded forward operator.

fn solve_rhs( &self, lower: MatRef<'_, Complex>, upper: MatRef<'_, Complex>, rhs: DMatrix<Complex>, adjoint: bool, ) -> Result<DMatrix<Complex>>

Reuse the recorded solve strategy for either derivative direction.

Trait Implementations

impl Debug for InternalSolve

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

Formats the value using the given formatter.

impl SavedState for InternalSolve

fn save_state(&self) -> Result<Vec<u8>>

Copy the saved numerical state into its fixed byte layout.

fn from_state(bytes: &[u8]) -> Result<Self>

Restore numerical state from bytes produced by this extension.