Struct Ilu0 

pub struct Ilu0<I, T> { /* private fields */ }

Numeric ILU(0) factor of a CSC matrix A.

Holds the L and U value buffers together with the symbolic structure and the dense workspace used during refactorisation. After construction the factor can be refactored against any matrix with the same sparsity pattern — Ilu0::refactorize performs zero heap allocations.

Implementations

impl<I: Index, T: ComplexField> Ilu0<I, T>

pub fn new_with_symbolic(symbolic: SymbolicIlu0<I>) -> Self

Allocate the value buffers and dense workspace for a given symbolic factor.

The returned factor is not yet populated — call Ilu0::refactorize to fill it with values from a matrix matching symbolic’s pattern.

pub fn try_new(a: SparseColMatRef<'_, I, T>) -> Result<Self, Ilu0Error>

Build a fully populated ILU(0) factor from a CSC matrix A in one shot.

Allocates internally for the symbolic structure, the L/U value buffers, and the dense workspace.

pub fn dim(&self) -> usize

Dimension n of the factored matrix.

pub fn symbolic(&self) -> &SymbolicIlu0<I>

Borrow the symbolic factor.

pub fn refactorize( &mut self, a: SparseColMatRef<'_, I, T>, ) -> Result<(), Ilu0Error>

Refactorise against a new matrix with the same sparsity pattern.

Performs zero heap allocations.

Errors

• Ilu0Error::PatternMismatch if a’s shape or column lengths disagree with the symbolic structure, or if the diagonal entries land at unexpected positions.
• Ilu0Error::ZeroPivot if a pivot of magnitude zero is encountered during elimination.

pub fn l_view(&self) -> SparseColMatRef<'_, I, T>

Construct a SparseColMatRef view over the L factor (unit lower triangular, diagonal stored first in each column with value 1).

pub fn u_view(&self) -> SparseColMatRef<'_, I, T>

Construct a SparseColMatRef view over the U factor (upper triangular, diagonal stored last in each column).

Trait Implementations

impl<I, T> BiLinOp<T> for Ilu0<I, T>

where I: Index, T: ComplexField + Debug + Sync,

fn transpose_apply_scratch(&self, _rhs_ncols: usize, _par: Par) -> StackReq

computes the workspace layout required to apply the transpose or adjoint o self to a matrix with rhs_ncols columns

fn transpose_apply( &self, out: MatMut<'_, T>, rhs: MatRef<'_, T>, par: Par, _stack: &mut MemStack, )

applies the transpose of self to rhs, and stores the result in out

fn adjoint_apply( &self, out: MatMut<'_, T>, rhs: MatRef<'_, T>, par: Par, _stack: &mut MemStack, )

applies the adjoint of self to rhs, and stores the result in out

impl<I, T> BiPrecond<T> for Ilu0<I, T>

where I: Index, T: ComplexField + Debug + Sync,

fn transpose_apply_in_place_scratch( &self, _rhs_ncols: usize, _par: Par, ) -> StackReq

computes the workspace layout required to apply the transpose or adjoint of self to a matrix with rhs_ncols columns in place

fn transpose_apply_in_place( &self, rhs: MatMut<'_, T>, par: Par, _stack: &mut MemStack, )

applies the transpose of self to rhs, and stores the result in rhs

fn adjoint_apply_in_place( &self, rhs: MatMut<'_, T>, par: Par, _stack: &mut MemStack, )

applies the adjoint of self to rhs, and stores the result in rhs

impl<I: Clone, T: Clone> Clone for Ilu0<I, T>

fn clone(&self) -> Ilu0<I, T>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<I: Debug, T: Debug> Debug for Ilu0<I, T>

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

Formats the value using the given formatter.

impl<I, T> LinOp<T> for Ilu0<I, T>

where I: Index, T: ComplexField + Debug + Sync,

fn apply_scratch(&self, _rhs_ncols: usize, _par: Par) -> StackReq

computes the workspace layout required to apply self or the conjugate o self to a matrix with rhs_ncols columns

fn nrows(&self) -> usize

output dimension of the operator

fn ncols(&self) -> usize

input dimension of the operator

fn apply( &self, out: MatMut<'_, T>, rhs: MatRef<'_, T>, par: Par, _stack: &mut MemStack, )

applies self to rhs, and stores the result in out

fn conj_apply( &self, out: MatMut<'_, T>, rhs: MatRef<'_, T>, par: Par, _stack: &mut MemStack, )

applies the conjugate of self to rhs, and stores the result in out

impl<I, T> Precond<T> for Ilu0<I, T>

where I: Index, T: ComplexField + Debug + Sync,

fn apply_in_place_scratch(&self, _rhs_ncols: usize, _par: Par) -> StackReq

computes the workspace layout required to apply self or the conjugate of self to a matrix with rhs_ncols columns in place

fn apply_in_place(&self, rhs: MatMut<'_, T>, par: Par, _stack: &mut MemStack)

applies self to rhs, and stores the result in rhs

fn conj_apply_in_place( &self, rhs: MatMut<'_, T>, par: Par, _stack: &mut MemStack, )

applies the conjugate of self to rhs, and stores the result in rhs