Public Types | Public Member Functions | Protected Member Functions | Protected Attributes
PastixLDLT< _MatrixType, _UpLo > Class Template Reference

A sparse direct supernodal Cholesky (LLT) factorization and solver based on the PaStiX library. More...

#include <PaStiXSupport.h>

+ Inheritance diagram for PastixLDLT< _MatrixType, _UpLo >:

List of all members.

Public Types

enum  { UpLo }
typedef
internal::pastix_traits
< PastixLDLT< _MatrixType,
_UpLo > >::MatrixType 
_MatrixType
typedef PastixBase< PastixLDLT
< MatrixType, _UpLo > > 
Base
typedef MatrixType::Index Index
typedef _MatrixType MatrixType
typedef MatrixType::RealScalar RealScalar
typedef MatrixType::Scalar Scalar
typedef Matrix< Scalar,
Dynamic, 1 > 
Vector

Public Member Functions

bool _solve (const MatrixBase< Rhs > &b, MatrixBase< Dest > &x) const
void _solve_sparse (const Rhs &b, SparseMatrix< DestScalar, DestOptions, DestIndex > &dest) const
PastixLDLT< _MatrixType, _UpLo > & analyzePattern (MatrixType &mat)
void analyzePattern (const MatrixType &matrix)
Index cols () const
void compute (const MatrixType &matrix)
PastixLDLT< _MatrixType, _UpLo > & derived ()
const PastixLDLT< _MatrixType,
_UpLo > & 
derived () const
Array< RealScalar, IPARM_SIZE, 1 > & dparm ()
double & dparm (int idxparam)
PastixLDLT< _MatrixType, _UpLo > & factorize (MatrixType &mat)
void factorize (const MatrixType &matrix)
ComputationInfo info () const
 Reports whether previous computation was successful.
Array< Index, IPARM_SIZE, 1 > & iparm ()
intiparm (int idxparam)
void PastixInit ()
 PastixLDLT ()
 PastixLDLT (const MatrixType &matrix)
Index rows () const
const internal::solve_retval
< PastixBase, Rhs > 
solve (const MatrixBase< Rhs > &b) const
const
internal::sparse_solve_retval
< PastixBase, Rhs > 
solve (const SparseMatrixBase< Rhs > &b) const

Protected Member Functions

PastixLDLT< _MatrixType, _UpLo > & compute (MatrixType &mat)
void PastixDestroy ()

Protected Attributes

int m_amalgamation
int m_analysisIsOk
int m_comm
Matrix< double, DPARM_SIZE, 1 > m_dparm
int m_factorizationIsOk
bool m_hasTranspose
ComputationInfo m_info
int m_initisOk
Matrix< Index, Dynamic, 1 > m_invp
Matrix< Index, IPARM_SIZE, 1 > m_iparm
bool m_isInitialized
SparseMatrix< Scalar, ColMajorm_mat_null
int m_ordering
pastix_data_t * m_pastixdata
Matrix< Index, Dynamic, 1 > m_perm
int m_size
SparseMatrix< Scalar, ColMajorm_StrMatTrans
Matrix< Scalar, Dynamic, 1 > m_vec_null

Detailed Description

template<typename _MatrixType, int _UpLo>
class Eigen::PastixLDLT< _MatrixType, _UpLo >

A sparse direct supernodal Cholesky (LLT) factorization and solver based on the PaStiX library.

This class is used to solve the linear systems A.X = B via a LDL^T supernodal Cholesky factorization available in the PaStiX library. The matrix A should be symmetric and positive definite WARNING Selfadjoint complex matrices are not supported in the current version of PaStiX The vectors or matrices X and B can be either dense or sparse

Template Parameters:
MatrixTypethe type of the sparse matrix A, it must be a SparseMatrix<>
UpLoThe part of the matrix to use : Lower or Upper. The default is Lower as required by PaStiX
See also:
Solving linear problems

Member Typedef Documentation

typedef internal::pastix_traits<PastixLDLT< _MatrixType, _UpLo > >::MatrixType _MatrixType [inherited]
typedef PastixBase<PastixLDLT<MatrixType, _UpLo> > Base
typedef MatrixType::Index Index
typedef MatrixType::RealScalar RealScalar [inherited]
typedef MatrixType::Scalar Scalar
typedef Matrix<Scalar,Dynamic,1> Vector [inherited]

Member Enumeration Documentation

anonymous enum
Enumerator:
UpLo 

Constructor & Destructor Documentation

PastixLDLT ( ) [inline]
PastixLDLT ( const MatrixType matrix) [inline]

Member Function Documentation

bool _solve ( const MatrixBase< Rhs > &  b,
MatrixBase< Dest > &  x 
) const [inherited]
void _solve_sparse ( const Rhs &  b,
SparseMatrix< DestScalar, DestOptions, DestIndex > &  dest 
) const [inline, inherited]
PastixLDLT< _MatrixType, _UpLo > & analyzePattern ( MatrixType mat) [inherited]
void analyzePattern ( const MatrixType matrix) [inline]

Compute the LDL^T symbolic factorization of matrix using its sparsity pattern The result of this operation can be used with successive matrices having the same pattern as matrix

See also:
factorize()

References PastixBase< Derived >::analyzePattern(), PastixBase< PastixLDLT< _MatrixType, _UpLo > >::m_iparm, and PastixBase< Derived >::PastixInit().

Index cols ( void  ) const [inline, inherited]
PastixLDLT< _MatrixType, _UpLo > & compute ( MatrixType mat) [protected, inherited]
void compute ( const MatrixType matrix) [inline]
PastixLDLT< _MatrixType, _UpLo > & derived ( ) [inline, inherited]
const PastixLDLT< _MatrixType, _UpLo > & derived ( ) const [inline, inherited]
Array<RealScalar,IPARM_SIZE,1>& dparm ( ) [inline, inherited]

Returns a reference to the double vector DPARM of PaStiX parameters The statistics related to the different phases of factorization and solve are saved here as well

See also:
analyzePattern() factorize()
double& dparm ( int  idxparam) [inline, inherited]

Return a reference to a particular index parameter of the DPARM vector

See also:
dparm()
PastixLDLT< _MatrixType, _UpLo > & factorize ( MatrixType mat) [inherited]
void factorize ( const MatrixType matrix) [inline]

Compute the LDL^T supernodal numerical factorization of matrix

References PastixBase< Derived >::factorize(), and PastixBase< PastixLDLT< _MatrixType, _UpLo > >::m_iparm.

ComputationInfo info ( ) const [inline, inherited]

Reports whether previous computation was successful.

Returns:
Success if computation was succesful, NumericalIssue if the PaStiX reports a problem InvalidInput if the input matrix is invalid
See also:
iparm()
Array<Index,IPARM_SIZE,1>& iparm ( ) [inline, inherited]

Returns a reference to the integer vector IPARM of PaStiX parameters to modify the default parameters. The statistics related to the different phases of factorization and solve are saved here as well

See also:
analyzePattern() factorize()
int& iparm ( int  idxparam) [inline, inherited]

Return a reference to a particular index parameter of the IPARM vector

See also:
iparm()
void PastixDestroy ( ) [inline, protected, inherited]
void PastixInit ( ) [inherited]

Initialize the PaStiX data structure. A first call to this function fills iparm and dparm with the default PaStiX parameters

See also:
iparm() dparm()
Index rows ( void  ) const [inline, inherited]
const internal::solve_retval<PastixBase, Rhs> solve ( const MatrixBase< Rhs > &  b) const [inline, inherited]
Returns:
the solution x of $ A x = b $ using the current decomposition of A.
See also:
compute()
const internal::sparse_solve_retval<PastixBase, Rhs> solve ( const SparseMatrixBase< Rhs > &  b) const [inline, inherited]
Returns:
the solution x of $ A x = b $ using the current decomposition of A.
See also:
compute()

Member Data Documentation

int m_amalgamation [mutable, protected, inherited]
int m_analysisIsOk [protected, inherited]
int m_comm [mutable, protected, inherited]
Matrix<double,DPARM_SIZE,1> m_dparm [mutable, protected, inherited]
int m_factorizationIsOk [protected, inherited]
bool m_hasTranspose [mutable, protected, inherited]
ComputationInfo m_info [mutable, protected, inherited]
int m_initisOk [protected, inherited]
Matrix<Index,Dynamic,1> m_invp [mutable, protected, inherited]
Matrix<Index,IPARM_SIZE,1> m_iparm [mutable, protected, inherited]
bool m_isInitialized [protected, inherited]
SparseMatrix<Scalar, ColMajor> m_mat_null [mutable, protected, inherited]
int m_ordering [mutable, protected, inherited]
pastix_data_t* m_pastixdata [mutable, protected, inherited]
Matrix<Index,Dynamic,1> m_perm [mutable, protected, inherited]
int m_size [mutable, protected, inherited]
SparseMatrix<Scalar, ColMajor> m_StrMatTrans [mutable, protected, inherited]
Matrix<Scalar, Dynamic,1> m_vec_null [mutable, protected, inherited]

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