Math.LinearAlgebra.Sparse.Algorithms.Diagonal

Documentation

isDiag :: SparseMatrix α -> Bool

Checks if matrix has diagonal form

toDiag :: Integral α => SparseMatrix α -> (SparseMatrix α, SparseMatrix α, SparseMatrix α)

Transforms matrix to diagonal form and returns also two protocol matrices:

>>> let (d,t,u) = toDiag m  in t × m × (trans u) == d
True

t stores rows transformations and u — columns transformations