Struct nalgebra::linalg::SVD

pub struct SVD<T: ComplexField, R: DimMin<C>, C: Dim> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,  {
    pub u: Option<OMatrix<T, R, DimMinimum<R, C>>>,
    pub v_t: Option<OMatrix<T, DimMinimum<R, C>, C>>,
    pub singular_values: OVector<T::RealField, DimMinimum<R, C>>,
}

Singular Value Decomposition of a general matrix.

Fields

u: Option<OMatrix<T, R, DimMinimum<R, C>>>

The left-singular vectors U of this SVD.

v_t: Option<OMatrix<T, DimMinimum<R, C>, C>>

The right-singular vectors V^t of this SVD.

singular_values: OVector<T::RealField, DimMinimum<R, C>>

The singular values of this SVD.

Implementations

impl<T: ComplexField, R: DimMin<C>, C: Dim> SVD<T, R, C> where
    DimMinimum<R, C>: DimSub<U1>,
    DefaultAllocator: Allocator<T, R, C> + Allocator<T, C> + Allocator<T, R> + Allocator<T, DimDiff<DimMinimum<R, C>, U1>> + Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>> + Allocator<T::RealField, DimDiff<DimMinimum<R, C>, U1>>, 

pub fn new(matrix: OMatrix<T, R, C>, compute_u: bool, compute_v: bool) -> Self

Computes the Singular Value Decomposition of matrix using implicit shift.

pub fn try_new(
    matrix: OMatrix<T, R, C>,
    compute_u: bool,
    compute_v: bool,
    eps: T::RealField,
    max_niter: usize
) -> Option<Self>

Attempts to compute the Singular Value Decomposition of matrix using implicit shift.

Arguments

• compute_u − set this to true to enable the computation of left-singular vectors.
• compute_v − set this to true to enable the computation of right-singular vectors.
• eps − tolerance used to determine when a value converged to 0.
• max_niter − maximum total number of iterations performed by the algorithm. If this number of iteration is exceeded, None is returned. If niter == 0, then the algorithm continues indefinitely until convergence.

pub fn rank(&self, eps: T::RealField) -> usize

Computes the rank of the decomposed matrix, i.e., the number of singular values greater than eps.

pub fn recompose(self) -> Result<OMatrix<T, R, C>, &'static str>

Rebuild the original matrix.

This is useful if some of the singular values have been manually modified. Returns Err if the right- and left- singular vectors have not been computed at construction-time.

pub fn pseudo_inverse(
    self,
    eps: T::RealField
) -> Result<OMatrix<T, C, R>, &'static str> where
    DefaultAllocator: Allocator<T, C, R>, 

Computes the pseudo-inverse of the decomposed matrix.

Any singular value smaller than eps is assumed to be zero. Returns Err if the right- and left- singular vectors have not been computed at construction-time.

pub fn solve<R2: Dim, C2: Dim, S2>(
    &self,
    b: &Matrix<T, R2, C2, S2>,
    eps: T::RealField
) -> Result<OMatrix<T, C, C2>, &'static str> where
    S2: Storage<T, R2, C2>,
    DefaultAllocator: Allocator<T, C, C2> + Allocator<T, DimMinimum<R, C>, C2>,
    ShapeConstraint: SameNumberOfRows<R, R2>, 

Solves the system self * x = b where self is the decomposed matrix and x the unknown.

Any singular value smaller than eps is assumed to be zero. Returns Err if the singular vectors U and V have not been computed.

Trait Implementations

impl<T: Clone + ComplexField, R: Clone + DimMin<C>, C: Clone + Dim> Clone for SVD<T, R, C> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,
    T::RealField: Clone, 

fn clone(&self) -> SVD<T, R, C>

Returns a copy of the value.

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

Performs copy-assignment from source.

impl<T: Debug + ComplexField, R: Debug + DimMin<C>, C: Debug + Dim> Debug for SVD<T, R, C> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,
    T::RealField: Debug, 

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

Formats the value using the given formatter.

impl<'de, T: ComplexField, R: DimMin<C>, C: Dim> Deserialize<'de> for SVD<T, R, C> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,
    DefaultAllocator: Allocator<T::RealField, DimMinimum<R, C>> + Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>>,
    OMatrix<T, R, DimMinimum<R, C>>: Deserialize<'de>,
    OMatrix<T, DimMinimum<R, C>, C>: Deserialize<'de>,
    OVector<T::RealField, DimMinimum<R, C>>: Deserialize<'de>, 

fn deserialize<__D>(__deserializer: __D) -> Result<Self, __D::Error> where
    __D: Deserializer<'de>, 

Deserialize this value from the given Serde deserializer.

impl<T: ComplexField, R: DimMin<C>, C: Dim> Serialize for SVD<T, R, C> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,
    DefaultAllocator: Allocator<T::RealField, DimMinimum<R, C>> + Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>>,
    OMatrix<T, R, DimMinimum<R, C>>: Serialize,
    OMatrix<T, DimMinimum<R, C>, C>: Serialize,
    OVector<T::RealField, DimMinimum<R, C>>: Serialize, 

fn serialize<__S>(&self, __serializer: __S) -> Result<__S::Ok, __S::Error> where
    __S: Serializer, 

Serialize this value into the given Serde serializer.

impl<T: ComplexField, R: DimMin<C>, C: Dim> Copy for SVD<T, R, C> where
    DefaultAllocator: Allocator<T, DimMinimum<R, C>, C> + Allocator<T, R, DimMinimum<R, C>> + Allocator<T::RealField, DimMinimum<R, C>>,
    OMatrix<T, R, DimMinimum<R, C>>: Copy,
    OMatrix<T, DimMinimum<R, C>, C>: Copy,
    OVector<T::RealField, DimMinimum<R, C>>: Copy,