Sparse left-looking QR factorization with numerical column pivoting. More...

#include <SparseQR.h>

Inheritance diagram for Eigen::SparseQR< _MatrixType, _OrderingType >:

Public Types
enum	{ ColsAtCompileTime = MatrixType::ColsAtCompileTime, MaxColsAtCompileTime = MatrixType::MaxColsAtCompileTime }

typedef Matrix< StorageIndex, Dynamic, 1 >	IndexVector

typedef _MatrixType	MatrixType

typedef _OrderingType	OrderingType

typedef PermutationMatrix< Dynamic, Dynamic, StorageIndex >	PermutationType

typedef SparseMatrix< Scalar, ColMajor, StorageIndex >	QRMatrixType

typedef MatrixType::RealScalar	RealScalar

typedef MatrixType::Scalar	Scalar

typedef Matrix< Scalar, Dynamic, 1 >	ScalarVector

typedef MatrixType::StorageIndex	StorageIndex

Public Member Functions
template<typename Rhs , typename Dest >
bool	_solve_impl (const MatrixBase< Rhs > &B, MatrixBase< Dest > &dest) const

void	_sort_matrix_Q ()

void	analyzePattern (const MatrixType &mat)
	Preprocessing step of a QR factorization. More...

Index	cols () const

const PermutationType &	colsPermutation () const

void	compute (const MatrixType &mat)

void	factorize (const MatrixType &mat)
	Performs the numerical QR factorization of the input matrix. More...

ComputationInfo	info () const
	Reports whether previous computation was successful. More...

std::string	lastErrorMessage () const

SparseQRMatrixQReturnType< SparseQR >	matrixQ () const

const QRMatrixType &	matrixR () const

Index	rank () const

Index	rows () const

void	setPivotThreshold (const RealScalar &threshold)

template<typename Rhs >
const Solve< SparseQR, Rhs >	solve (const MatrixBase< Rhs > &B) const

template<typename Rhs >
const Solve< SparseQR, Rhs >	solve (const SparseMatrixBase< Rhs > &B) const

	SparseQR ()

	SparseQR (const MatrixType &mat)

Public Member Functions inherited from Eigen::SparseSolverBase< SparseQR< _MatrixType, _OrderingType > >
void	_solve_impl (const SparseMatrixBase< Rhs > &b, SparseMatrixBase< Dest > &dest) const

SparseQR< _MatrixType, _OrderingType > &	derived ()

const SparseQR< _MatrixType, _OrderingType > &	derived () const

const Solve< SparseQR< _MatrixType, _OrderingType >, Rhs >	solve (const MatrixBase< Rhs > &b) const

const Solve< SparseQR< _MatrixType, _OrderingType >, Rhs >	solve (const SparseMatrixBase< Rhs > &b) const

	SparseSolverBase ()

	~SparseSolverBase ()

Protected Types
typedef SparseSolverBase< SparseQR< _MatrixType, _OrderingType > >	Base

Protected Attributes
bool	m_analysisIsok

IndexVector	m_etree

bool	m_factorizationIsok

IndexVector	m_firstRowElt

ScalarVector	m_hcoeffs

ComputationInfo	m_info

bool	m_isEtreeOk

bool	m_isQSorted

std::string	m_lastError

Index	m_nonzeropivots

PermutationType	m_outputPerm_c

PermutationType	m_perm_c

PermutationType	m_pivotperm

QRMatrixType	m_pmat

QRMatrixType	m_Q

QRMatrixType	m_R

RealScalar	m_threshold

bool	m_useDefaultThreshold

Protected Attributes inherited from Eigen::SparseSolverBase< SparseQR< _MatrixType, _OrderingType > >
bool	m_isInitialized

Friends
template<typename , typename >
struct	SparseQR_QProduct

Detailed Description

template<typename _MatrixType, typename _OrderingType>
class Eigen::SparseQR< _MatrixType, _OrderingType >

Sparse left-looking QR factorization with numerical column pivoting.

This class implements a left-looking QR decomposition of sparse matrices with numerical column pivoting. When a column has a norm less than a given tolerance it is implicitly permuted to the end. The QR factorization thus obtained is given by A*P = Q*R where R is upper triangular or trapezoidal.

P is the column permutation which is the product of the fill-reducing and the numerical permutations. Use colsPermutation() to get it.

Q is the orthogonal matrix represented as products of Householder reflectors. Use matrixQ() to get an expression and matrixQ().adjoint() to get the adjoint. You can then apply it to a vector.

R is the sparse triangular or trapezoidal matrix. The later occurs when A is rank-deficient. matrixR().topLeftCorner(rank(), rank()) always returns a triangular factor of full rank.

Template Parameters

_MatrixType	The type of the sparse matrix A, must be a column-major SparseMatrix<>
_OrderingType	The fill-reducing ordering method. See the OrderingMethods module for the list of built-in and external ordering methods.

The numerical pivoting strategy and default threshold are the same as in SuiteSparse QR, and detailed in the following paper: Tim Davis, "Algorithm 915, SuiteSparseQR: Multifrontal Multithreaded Rank-Revealing Sparse QR Factorization, ACM Trans. on Math. Soft. 38(1), 2011. Even though it is qualified as "rank-revealing", this strategy might fail for some rank deficient problems. When this class is used to solve linear or least-square problems it is thus strongly recommended to check the accuracy of the computed solution. If it failed, it usually helps to increase the threshold with setPivotThreshold.

Warning: The input sparse matrix A must be in compressed mode (see SparseMatrix::makeCompressed()).; For complex matrices matrixQ().transpose() will actually return the adjoint matrix.

Definition at line 16 of file SparseQR.h.

Member Typedef Documentation

◆ Base

template<typename _MatrixType, typename _OrderingType>

typedef SparseSolverBase<SparseQR<_MatrixType,_OrderingType> > Eigen::SparseQR< _MatrixType, _OrderingType >::Base

protected

Definition at line 87 of file SparseQR.h.

◆ IndexVector

template<typename _MatrixType, typename _OrderingType>

typedef Matrix<StorageIndex, Dynamic, 1> Eigen::SparseQR< _MatrixType, _OrderingType >::IndexVector

Definition at line 97 of file SparseQR.h.

◆ MatrixType

template<typename _MatrixType, typename _OrderingType>

typedef _MatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::MatrixType

Definition at line 91 of file SparseQR.h.

◆ OrderingType

template<typename _MatrixType, typename _OrderingType>

typedef _OrderingType Eigen::SparseQR< _MatrixType, _OrderingType >::OrderingType

Definition at line 92 of file SparseQR.h.

◆ PermutationType

template<typename _MatrixType, typename _OrderingType>

typedef PermutationMatrix<Dynamic, Dynamic, StorageIndex> Eigen::SparseQR< _MatrixType, _OrderingType >::PermutationType

Definition at line 99 of file SparseQR.h.

◆ QRMatrixType

template<typename _MatrixType, typename _OrderingType>

typedef SparseMatrix<Scalar,ColMajor,StorageIndex> Eigen::SparseQR< _MatrixType, _OrderingType >::QRMatrixType

Definition at line 96 of file SparseQR.h.

◆ RealScalar

template<typename _MatrixType, typename _OrderingType>

typedef MatrixType::RealScalar Eigen::SparseQR< _MatrixType, _OrderingType >::RealScalar

Definition at line 94 of file SparseQR.h.

◆ Scalar

template<typename _MatrixType, typename _OrderingType>

typedef MatrixType::Scalar Eigen::SparseQR< _MatrixType, _OrderingType >::Scalar

Definition at line 93 of file SparseQR.h.

◆ ScalarVector

template<typename _MatrixType, typename _OrderingType>

typedef Matrix<Scalar, Dynamic, 1> Eigen::SparseQR< _MatrixType, _OrderingType >::ScalarVector

Definition at line 98 of file SparseQR.h.

◆ StorageIndex

template<typename _MatrixType, typename _OrderingType>

typedef MatrixType::StorageIndex Eigen::SparseQR< _MatrixType, _OrderingType >::StorageIndex

Definition at line 95 of file SparseQR.h.

Member Enumeration Documentation

◆ anonymous enum

template<typename _MatrixType, typename _OrderingType>

anonymous enum

Enumerator
ColsAtCompileTime
MaxColsAtCompileTime

Definition at line 101 of file SparseQR.h.

Constructor & Destructor Documentation

◆ SparseQR() [1/2]

template<typename _MatrixType, typename _OrderingType>

Eigen::SparseQR< _MatrixType, _OrderingType >::SparseQR ( )

inline

Definition at line 107 of file SparseQR.h.

◆ SparseQR() [2/2]

template<typename _MatrixType, typename _OrderingType>

Eigen::SparseQR< _MatrixType, _OrderingType >::SparseQR ( const MatrixType & mat )

inlineexplicit

Construct a QR factorization of the matrix mat.

Warning: The matrix mat must be in compressed mode (see SparseMatrix::makeCompressed()).

See also: compute()

Definition at line 116 of file SparseQR.h.

Member Function Documentation

◆ _solve_impl()

template<typename _MatrixType, typename _OrderingType>

template<typename Rhs , typename Dest >

bool Eigen::SparseQR< _MatrixType, _OrderingType >::_solve_impl	(	const MatrixBase< Rhs > &	B,
		MatrixBase< Dest > &	dest
	)		const

inline

Definition at line 205 of file SparseQR.h.

◆ _sort_matrix_Q()

template<typename _MatrixType, typename _OrderingType>

void Eigen::SparseQR< _MatrixType, _OrderingType >::_sort_matrix_Q ( )

inline

Definition at line 276 of file SparseQR.h.

◆ analyzePattern()

template<typename MatrixType , typename OrderingType >

void Eigen::SparseQR< MatrixType, OrderingType >::analyzePattern ( const MatrixType & mat )

Preprocessing step of a QR factorization.

Warning: The matrix mat must be in compressed mode (see SparseMatrix::makeCompressed()).

In this step, the fill-reducing permutation is computed and applied to the columns of A and the column elimination tree is computed as well. Only the sparsity pattern of mat is exploited.

Note: In this step it is assumed that there is no empty row in the matrix mat.

Definition at line 320 of file SparseQR.h.

◆ cols()

template<typename _MatrixType, typename _OrderingType>

Index Eigen::SparseQR< _MatrixType, _OrderingType >::cols ( void ) const

inline

Returns: the number of columns of the represented matrix.

Definition at line 141 of file SparseQR.h.

◆ colsPermutation()

template<typename _MatrixType, typename _OrderingType>

const PermutationType& Eigen::SparseQR< _MatrixType, _OrderingType >::colsPermutation ( ) const

inline

Returns: a const reference to the column permutation P that was applied to A such that A*P = Q*R It is the combination of the fill-in reducing permutation and numerical column pivoting.

Definition at line 192 of file SparseQR.h.

◆ compute()

template<typename _MatrixType, typename _OrderingType>

void Eigen::SparseQR< _MatrixType, _OrderingType >::compute ( const MatrixType & mat )

inline

Computes the QR factorization of the sparse matrix mat.

Warning: The matrix mat must be in compressed mode (see SparseMatrix::makeCompressed()).

See also: analyzePattern(), factorize()

Definition at line 127 of file SparseQR.h.

◆ factorize()

template<typename MatrixType , typename OrderingType >

void Eigen::SparseQR< MatrixType, OrderingType >::factorize ( const MatrixType & mat )

Performs the numerical QR factorization of the input matrix.

The function SparseQR::analyzePattern(const MatrixType&) must have been called beforehand with a matrix having the same sparsity pattern than mat.

Parameters

mat	The sparse column-major matrix

Definition at line 361 of file SparseQR.h.

◆ info()

template<typename _MatrixType, typename _OrderingType>

ComputationInfo Eigen::SparseQR< _MatrixType, _OrderingType >::info ( ) const

inline

Reports whether previous computation was successful.

Returns: Success if computation was successful, NumericalIssue if the QR factorization reports a numerical problem InvalidInput if the input matrix is invalid

See also: iparm()

Definition at line 268 of file SparseQR.h.

◆ lastErrorMessage()

template<typename _MatrixType, typename _OrderingType>

std::string Eigen::SparseQR< _MatrixType, _OrderingType >::lastErrorMessage ( ) const

inline

Returns: A string describing the type of error. This method is provided to ease debugging, not to handle errors.

Definition at line 201 of file SparseQR.h.

◆ matrixQ()

template<typename _MatrixType, typename _OrderingType>

SparseQRMatrixQReturnType<SparseQR> Eigen::SparseQR< _MatrixType, _OrderingType >::matrixQ ( void ) const

inline

Returns: an expression of the matrix Q as products of sparse Householder reflectors. The common usage of this function is to apply it to a dense matrix or vector
VectorXd B1, B2;
// Initialize B1
B2 = matrixQ() * B1;

To get a plain SparseMatrix representation of Q:

SparseMatrix<double> Q;

Q = SparseQR<SparseMatrix<double> >(A).matrixQ();

Internally, this call simply performs a sparse product between the matrix Q and a sparse identity matrix. However, due to the fact that the sparse reflectors are stored unsorted, two transpositions are needed to sort them before performing the product.

Definition at line 186 of file SparseQR.h.

◆ matrixR()

template<typename _MatrixType, typename _OrderingType>

const QRMatrixType& Eigen::SparseQR< _MatrixType, _OrderingType >::matrixR ( ) const

inline

Returns: a const reference to the sparse upper triangular matrix R of the QR factorization.

Warning: The entries of the returned matrix are not sorted. This means that using it in algorithms expecting sorted entries will fail. This include random coefficient accesses (SpaseMatrix::coeff()), and coefficient-wise operations. Matrix products and triangular solves are fine though.

To sort the entries, you can assign it to a row-major matrix, and if a column-major matrix is required, you can copy it again:

SparseMatrix<double>          R  = qr.matrixR();  // column-major, not sorted!
SparseMatrix<double,RowMajor> Rr = qr.matrixR();  // row-major, sorted
SparseMatrix<double>          Rc = Rr;            // column-major, sorted

Definition at line 156 of file SparseQR.h.

◆ rank()

template<typename _MatrixType, typename _OrderingType>

Index Eigen::SparseQR< _MatrixType, _OrderingType >::rank ( ) const

inline

Returns: the number of non linearly dependent columns as determined by the pivoting threshold.

See also: setPivotThreshold()

Definition at line 162 of file SparseQR.h.

◆ rows()

template<typename _MatrixType, typename _OrderingType>

Index Eigen::SparseQR< _MatrixType, _OrderingType >::rows ( void ) const

inline

Returns: the number of rows of the represented matrix.

Definition at line 137 of file SparseQR.h.

◆ setPivotThreshold()

template<typename _MatrixType, typename _OrderingType>

void Eigen::SparseQR< _MatrixType, _OrderingType >::setPivotThreshold ( const RealScalar & threshold )

inline

Sets the threshold that is used to determine linearly dependent columns during the factorization.

In practice, if during the factorization the norm of the column that has to be eliminated is below this threshold, then the entire column is treated as zero, and it is moved at the end.

Definition at line 235 of file SparseQR.h.

◆ solve() [1/2]

template<typename _MatrixType, typename _OrderingType>

template<typename Rhs >

const Solve<SparseQR, Rhs> Eigen::SparseQR< _MatrixType, _OrderingType >::solve ( const MatrixBase< Rhs > & B ) const

inline

Returns: the solution X of using the current decomposition of A.

See also: compute()

Definition at line 246 of file SparseQR.h.

◆ solve() [2/2]

template<typename _MatrixType, typename _OrderingType>

template<typename Rhs >

const Solve<SparseQR, Rhs> Eigen::SparseQR< _MatrixType, _OrderingType >::solve ( const SparseMatrixBase< Rhs > & B ) const

inline

Definition at line 253 of file SparseQR.h.

Friends And Related Function Documentation

◆ SparseQR_QProduct

template<typename _MatrixType, typename _OrderingType>

template<typename , typename >

friend struct SparseQR_QProduct

friend

Definition at line 306 of file SparseQR.h.

Member Data Documentation

◆ m_analysisIsok

template<typename _MatrixType, typename _OrderingType>

bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_analysisIsok

protected

Definition at line 287 of file SparseQR.h.

◆ m_etree

template<typename _MatrixType, typename _OrderingType>

IndexVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_etree

protected

Definition at line 301 of file SparseQR.h.

◆ m_factorizationIsok

template<typename _MatrixType, typename _OrderingType>

bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_factorizationIsok

protected

Definition at line 288 of file SparseQR.h.

◆ m_firstRowElt

template<typename _MatrixType, typename _OrderingType>

IndexVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_firstRowElt

protected

Definition at line 302 of file SparseQR.h.

◆ m_hcoeffs

template<typename _MatrixType, typename _OrderingType>

ScalarVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_hcoeffs

protected

Definition at line 294 of file SparseQR.h.

◆ m_info

template<typename _MatrixType, typename _OrderingType>

ComputationInfo Eigen::SparseQR< _MatrixType, _OrderingType >::m_info

mutableprotected

Definition at line 289 of file SparseQR.h.

◆ m_isEtreeOk

template<typename _MatrixType, typename _OrderingType>

bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_isEtreeOk

protected

Definition at line 304 of file SparseQR.h.

◆ m_isQSorted

template<typename _MatrixType, typename _OrderingType>

bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_isQSorted

protected

Definition at line 303 of file SparseQR.h.

◆ m_lastError

template<typename _MatrixType, typename _OrderingType>

std::string Eigen::SparseQR< _MatrixType, _OrderingType >::m_lastError

protected

Definition at line 290 of file SparseQR.h.

◆ m_nonzeropivots

template<typename _MatrixType, typename _OrderingType>

Index Eigen::SparseQR< _MatrixType, _OrderingType >::m_nonzeropivots

protected

Definition at line 300 of file SparseQR.h.

◆ m_outputPerm_c

template<typename _MatrixType, typename _OrderingType>

PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_outputPerm_c

protected

Definition at line 297 of file SparseQR.h.

◆ m_perm_c

template<typename _MatrixType, typename _OrderingType>

PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_perm_c

protected

Definition at line 295 of file SparseQR.h.

◆ m_pivotperm

template<typename _MatrixType, typename _OrderingType>

PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_pivotperm

protected

Definition at line 296 of file SparseQR.h.

◆ m_pmat

template<typename _MatrixType, typename _OrderingType>

QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_pmat

protected

Definition at line 291 of file SparseQR.h.

◆ m_Q

template<typename _MatrixType, typename _OrderingType>

QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_Q

protected

Definition at line 293 of file SparseQR.h.

◆ m_R

template<typename _MatrixType, typename _OrderingType>

QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_R

protected

Definition at line 292 of file SparseQR.h.

◆ m_threshold

template<typename _MatrixType, typename _OrderingType>

RealScalar Eigen::SparseQR< _MatrixType, _OrderingType >::m_threshold

protected

Definition at line 298 of file SparseQR.h.

◆ m_useDefaultThreshold

template<typename _MatrixType, typename _OrderingType>

bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_useDefaultThreshold

protected

Definition at line 299 of file SparseQR.h.

The documentation for this class was generated from the following file:

SparseQR.h

Public Types

Public Member Functions

Protected Types

Protected Attributes

Friends

Detailed Description

template<typename _MatrixType, typename _OrderingType> class Eigen::SparseQR< _MatrixType, _OrderingType >

Member Typedef Documentation

◆ Base

◆ IndexVector

◆ MatrixType

◆ OrderingType

◆ PermutationType

◆ QRMatrixType

◆ RealScalar

◆ Scalar

◆ ScalarVector

◆ StorageIndex

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ SparseQR() [1/2]

◆ SparseQR() [2/2]

Member Function Documentation

◆ _solve_impl()

◆ _sort_matrix_Q()

◆ analyzePattern()

◆ cols()

◆ colsPermutation()

◆ compute()

◆ factorize()

◆ info()

◆ lastErrorMessage()

◆ matrixQ()

◆ matrixR()

◆ rank()

◆ rows()

◆ setPivotThreshold()

◆ solve() [1/2]

◆ solve() [2/2]

Friends And Related Function Documentation

◆ SparseQR_QProduct

Member Data Documentation

◆ m_analysisIsok

◆ m_etree

◆ m_factorizationIsok

◆ m_firstRowElt

◆ m_hcoeffs

◆ m_info

◆ m_isEtreeOk

◆ m_isQSorted

◆ m_lastError

◆ m_nonzeropivots

◆ m_outputPerm_c

◆ m_perm_c

◆ m_pivotperm

◆ m_pmat

◆ m_Q

◆ m_R

◆ m_threshold

◆ m_useDefaultThreshold

template<typename _MatrixType, typename _OrderingType>
class Eigen::SparseQR< _MatrixType, _OrderingType >