Struct Householder 

pub struct Householder<A>
where
    A: Scalar,
{ /* private fields */ }

Iterative orthogonalizer using Householder reflection

Implementations

impl<A> Householder<A>

where A: Scalar + Lapack,

pub fn new(dim: usize, tol: <A as Scalar>::Real) -> Householder<A>

Create a new orthogonalizer

pub fn forward_reflection<S>(&self, a: &mut ArrayBase<S, Dim<[usize; 1]>>)

where S: DataMut<Elem = A>,

Take forward reflection P = P_l ... P_1

pub fn backward_reflection<S>(&self, a: &mut ArrayBase<S, Dim<[usize; 1]>>)

where S: DataMut<Elem = A>,

Take backward reflection P = P_1 ... P_l

Trait Implementations

impl<A> Clone for Householder<A>

where A: Clone + Scalar, <A as Scalar>::Real: Clone,

fn clone(&self) -> Householder<A>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<A> Debug for Householder<A>

where A: Debug + Scalar, <A as Scalar>::Real: Debug,

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

Formats the value using the given formatter.

impl<A> Orthogonalizer for Householder<A>

where A: Scalar + Lapack,

type Elem = A

fn dim(&self) -> usize

Dimension of input array

fn len(&self) -> usize

Number of cached basis

fn tolerance(&self) -> <A as Scalar>::Real

fn decompose<S>( &self, a: &mut ArrayBase<S, Dim<[usize; 1]>>, ) -> ArrayBase<OwnedRepr<A>, Dim<[usize; 1]>>

where S: DataMut<Elem = A>,

Decompose given vector into the span of current basis and its tangent space

fn coeff<S>( &self, a: ArrayBase<S, Dim<[usize; 1]>>, ) -> ArrayBase<OwnedRepr<A>, Dim<[usize; 1]>>

where S: Data<Elem = A>,

Calculate the coefficient to the current basis basis

fn div_append<S>( &mut self, a: &mut ArrayBase<S, Dim<[usize; 1]>>, ) -> AppendResult<A>

where S: DataMut<Elem = A>,

Add new vector if the residual is larger than relative tolerance, and return the residual vector

fn append<S>(&mut self, a: ArrayBase<S, Dim<[usize; 1]>>) -> AppendResult<A>

where S: Data<Elem = A>,

Add new vector if the residual is larger than relative tolerance

fn get_q(&self) -> ArrayBase<OwnedRepr<A>, Dim<[usize; 2]>>

Get Q-matrix of generated basis

fn is_full(&self) -> bool

check if the basis spans entire space

fn is_empty(&self) -> bool