Struct JacobiOperatorPreconditioner 

pub struct JacobiOperatorPreconditioner { /* private fields */ }

Diagonal (Jacobi) preconditioning built from any LinearOperator: z[i] = r[i] / A[i][i].

Stores the reciprocal diagonal, so applying it is one multiply per row. Cannot break down: a diagonal entry whose magnitude is below 1e-300 is treated as 1.0, matching crate::krylov::JacobiPreconditioner.

Units

The stored reciprocals carry the inverse of whatever units the matrix entries had; treated as dimensionless f64 here.

Implementations

impl JacobiOperatorPreconditioner

pub fn new<A: LinearOperator + ?Sized>(a: &A) -> Self

Extract and invert the diagonal of a.

Cost: one LinearOperator::diagonal call plus n divisions.

pub fn reciprocal_diagonal(&self) -> &[f64]

The stored reciprocal diagonal 1 / A[i][i], length n.

Trait Implementations

impl Clone for JacobiOperatorPreconditioner

fn clone(&self) -> Self

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for JacobiOperatorPreconditioner

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

Formats the value using the given formatter.

impl OperatorPreconditioner for JacobiOperatorPreconditioner

fn apply(&self, r: &[f64], z: &mut [f64])

Write z = M^{-1} r. Both slices have the system length; z is fully overwritten.