Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Friends
Eigen::SparseQR< _MatrixType, _OrderingType > Class Template Reference

Sparse left-looking rank-revealing QR factorization. More...

#include <SparseQR.h>

List of all members.

Public Types

typedef MatrixType::Index Index
typedef Matrix< Index, Dynamic, 1 > IndexVector
typedef _MatrixType MatrixType
typedef _OrderingType OrderingType
typedef PermutationMatrix
< Dynamic, Dynamic, Index > 		PermutationType
typedef SparseMatrix< Scalar,
ColMajor, Index > 		QRMatrixType
typedef MatrixType::RealScalar RealScalar
typedef MatrixType::Scalar Scalar
typedef Matrix< Scalar,
Dynamic, 1 > 		ScalarVector

Public Member Functions

template<typename Rhs , typename Dest >
bool _solve (const MatrixBase< Rhs > &B, MatrixBase< Dest > &dest) const
void analyzePattern (const MatrixType &mat)
 Preprocessing step of a QR factorization.
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.
ComputationInfo info () const
 Reports whether previous computation was successful.
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 internal::solve_retval
< SparseQR, Rhs > 		solve (const MatrixBase< Rhs > &B) const
template<typename Rhs >
const
internal::sparse_solve_retval
< SparseQR, Rhs > 		solve (const SparseMatrixBase< Rhs > &B) const
 SparseQR ()
 SparseQR (const MatrixType &mat)

Protected Member Functions

void sort_matrix_Q ()

Protected Attributes

bool m_analysisIsok
IndexVector m_etree
bool m_factorizationIsok
IndexVector m_firstRowElt
ScalarVector m_hcoeffs
ComputationInfo m_info
bool m_isInitialized
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

Friends

struct SparseQR_QProduct
struct SparseQRMatrixQReturnType

Detailed Description

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

Sparse left-looking rank-revealing QR factorization.

This class implements a left-looking rank-revealing QR decomposition of sparse matrices. 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 rank-revealing 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().transpose() to get the transpose. 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:
_MatrixTypeThe type of the sparse matrix A, must be a column-major SparseMatrix<>
_OrderingTypeThe fill-reducing ordering method. See the OrderingMethods module for the list of built-in and external ordering methods.

Definition at line 64 of file SparseQR.h.

Member Typedef Documentation

template<typename _MatrixType , typename _OrderingType >
typedef MatrixType::Index Eigen::SparseQR< _MatrixType, _OrderingType >::Index

Definition at line 71 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
typedef Matrix<Index, Dynamic, 1> Eigen::SparseQR< _MatrixType, _OrderingType >::IndexVector

Definition at line 73 of file SparseQR.h.

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

Definition at line 67 of file SparseQR.h.

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

Definition at line 68 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
typedef PermutationMatrix<Dynamic, Dynamic, Index> Eigen::SparseQR< _MatrixType, _OrderingType >::PermutationType

Definition at line 75 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
typedef SparseMatrix<Scalar,ColMajor,Index> Eigen::SparseQR< _MatrixType, _OrderingType >::QRMatrixType

Definition at line 72 of file SparseQR.h.

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

Definition at line 70 of file SparseQR.h.

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

Definition at line 69 of file SparseQR.h.

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

Definition at line 74 of file SparseQR.h.

Constructor & Destructor Documentation

template<typename _MatrixType , typename _OrderingType >
Eigen::SparseQR< _MatrixType, _OrderingType >::SparseQR ( ) [inline]

Definition at line 77 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
Eigen::SparseQR< _MatrixType, _OrderingType >::SparseQR ( const MatrixType &  mat) [inline]

Definition at line 80 of file SparseQR.h.

Member Function Documentation

template<typename _MatrixType , typename _OrderingType >
template<typename Rhs , typename Dest >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::_solve ( const MatrixBase< Rhs > &  B,
MatrixBase< Dest > &  dest 
) const [inline]

Definition at line 151 of file SparseQR.h.

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

Preprocessing step of a QR factorization.

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 sparcity 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 264 of file SparseQR.h.

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 98 of file SparseQR.h.

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 138 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
void Eigen::SparseQR< _MatrixType, _OrderingType >::compute ( const MatrixType &  mat) [inline]

Definition at line 84 of file SparseQR.h.

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 sparcity pattern than mat.

Parameters:
matThe sparse column-major matrix

Definition at line 300 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
ComputationInfo Eigen::SparseQR< _MatrixType, _OrderingType >::info ( ) const [inline]

Reports whether previous computation was successful.

Returns:
Success if computation was succesful, NumericalIssue if the QR factorization reports a numerical problem InvalidInput if the input matrix is invalid
See also:
iparm()

Definition at line 214 of file SparseQR.h.

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 147 of file SparseQR.h.

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 132 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
const QRMatrixType& Eigen::SparseQR< _MatrixType, _OrderingType >::matrixR ( void  ) const [inline]
Returns:
a const reference to the sparse upper triangular matrix R of the QR factorization.

Definition at line 102 of file SparseQR.h.

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 108 of file SparseQR.h.

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 94 of file SparseQR.h.

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 181 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
template<typename Rhs >
const internal::solve_retval<SparseQR, Rhs> Eigen::SparseQR< _MatrixType, _OrderingType >::solve ( const MatrixBase< Rhs > &  B) const [inline]
Returns:
the solution X of $ A X = B $ using the current decomposition of A.
See also:
compute()

Definition at line 192 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
template<typename Rhs >
const internal::sparse_solve_retval<SparseQR, Rhs> Eigen::SparseQR< _MatrixType, _OrderingType >::solve ( const SparseMatrixBase< Rhs > &  B) const [inline]

Definition at line 199 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
void Eigen::SparseQR< _MatrixType, _OrderingType >::sort_matrix_Q ( ) [inline, protected]

Definition at line 221 of file SparseQR.h.

Friends And Related Function Documentation

template<typename _MatrixType , typename _OrderingType >
friend struct SparseQR_QProduct [friend]

Definition at line 251 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
friend struct SparseQRMatrixQReturnType [friend]

Definition at line 252 of file SparseQR.h.

Member Data Documentation

template<typename _MatrixType , typename _OrderingType >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_analysisIsok [protected]

Definition at line 233 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
IndexVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_etree [protected]

Definition at line 247 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_factorizationIsok [protected]

Definition at line 234 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
IndexVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_firstRowElt [protected]

Definition at line 248 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
ScalarVector Eigen::SparseQR< _MatrixType, _OrderingType >::m_hcoeffs [protected]

Definition at line 240 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
ComputationInfo Eigen::SparseQR< _MatrixType, _OrderingType >::m_info [mutable, protected]

Definition at line 235 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_isInitialized [protected]

Definition at line 232 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_isQSorted [protected]

Definition at line 249 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
std::string Eigen::SparseQR< _MatrixType, _OrderingType >::m_lastError [protected]

Definition at line 236 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
Index Eigen::SparseQR< _MatrixType, _OrderingType >::m_nonzeropivots [protected]

Definition at line 246 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_outputPerm_c [protected]

Definition at line 243 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_perm_c [protected]

Definition at line 241 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
PermutationType Eigen::SparseQR< _MatrixType, _OrderingType >::m_pivotperm [protected]

Definition at line 242 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_pmat [protected]

Definition at line 237 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_Q [protected]

Definition at line 239 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
QRMatrixType Eigen::SparseQR< _MatrixType, _OrderingType >::m_R [protected]

Definition at line 238 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
RealScalar Eigen::SparseQR< _MatrixType, _OrderingType >::m_threshold [protected]

Definition at line 244 of file SparseQR.h.

template<typename _MatrixType , typename _OrderingType >
bool Eigen::SparseQR< _MatrixType, _OrderingType >::m_useDefaultThreshold [protected]

Definition at line 245 of file SparseQR.h.

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

acado
Author(s): Milan Vukov, Rien Quirynen
autogenerated on Thu Aug 27 2015 12:02:03