Struct DeflationSubspace 

pub struct DeflationSubspace<T: ComplexField> { /* private fields */ }

Pre-computed deflation data for deflated GMRES

Given deflation vectors W = [w₁,…,wᵣ] (n×r, r≪n) and operator A:

• AW = A·W columns
• E = Wᴴ·A·W (r×r dense), factored via LU

Provides projectors:

• Left: P(v) = v - AW·E⁻¹·Wᴴ·v
• Coarse correction: W·E⁻¹·Wᴴ·b
• Recovery: Q(v) = v - W·E⁻¹·Wᴴ·A·v

Implementations

impl<T: ComplexField> DeflationSubspace<T>

pub fn new<A: LinearOperator<T>>( w_columns: Vec<Array1<T>>, operator: &A, ) -> Result<Self, LuError>

Construct deflation subspace from deflation vectors and operator

Computes AW = A·W and factors E = Wᴴ·A·W via LU.

Arguments

• w_columns - Deflation vectors (should be orthonormal or near-orthonormal)
• operator - The linear operator A

pub fn num_vectors(&self) -> usize

Number of deflation vectors

pub fn apply_left_projector(&self, v: &Array1<T>) -> Array1<T>

Apply left deflation projector: P(v) = v - AW·E⁻¹·Wᴴ·v

pub fn coarse_correction(&self, b: &Array1<T>) -> Array1<T>

Compute coarse correction: W·E⁻¹·Wᴴ·b

pub fn apply_recovery<A: LinearOperator<T>>( &self, v: &Array1<T>, operator: &A, ) -> Array1<T>

Apply recovery operator: Q(v) = v - W·E⁻¹·Wᴴ·A·v

Trait Implementations

impl<T: Clone + ComplexField> Clone for DeflationSubspace<T>

fn clone(&self) -> DeflationSubspace<T>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<T: Debug + ComplexField> Debug for DeflationSubspace<T>

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

Formats the value using the given formatter.