Public Types | Public Member Functions | Protected Member Functions | Protected Attributes
PastixLLT< _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 PastixLLT< _MatrixType, _UpLo >:

List of all members.

Public Types

enum  { UpLo }
typedef
internal::pastix_traits
< PastixLLT< _MatrixType,
_UpLo > >::MatrixType 
_MatrixType
typedef PastixBase< PastixLLT
< MatrixType, _UpLo > > 
Base
typedef Base::ColSpMatrix ColSpMatrix
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
void analyzePattern (const MatrixType &matrix)
Index cols () const
void compute (const MatrixType &matrix)
PastixLLT< _MatrixType, _UpLo > & derived ()
const PastixLLT< _MatrixType,
_UpLo > & 
derived () const
Array< RealScalar, IPARM_SIZE, 1 > & dparm ()
double & dparm (int idxparam)
void factorize (const MatrixType &matrix)
ComputationInfo info () const
 Reports whether previous computation was successful.
Array< Index, IPARM_SIZE, 1 > & iparm ()
intiparm (int idxparam)
 PastixLLT ()
 PastixLLT (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

void analyzePattern (ColSpMatrix &mat)
void clean ()
void compute (ColSpMatrix &mat)
void factorize (ColSpMatrix &mat)
void grabMatrix (const MatrixType &matrix, ColSpMatrix &out)
void init ()

Protected Attributes

int m_analysisIsOk
int m_comm
Matrix< double, DPARM_SIZE, 1 > m_dparm
int m_factorizationIsOk
ComputationInfo m_info
int m_initisOk
Matrix< Index, Dynamic, 1 > m_invp
Matrix< int, IPARM_SIZE, 1 > m_iparm
bool m_isInitialized
pastix_data_t * m_pastixdata
Matrix< Index, Dynamic, 1 > m_perm
int m_size

Detailed Description

template<typename _MatrixType, int _UpLo>
class Eigen::PastixLLT< _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 LL^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<PastixLLT< _MatrixType, _UpLo > >::MatrixType _MatrixType
inherited
typedef PastixBase<PastixLLT<MatrixType, _UpLo> > Base
typedef MatrixType::Index Index
inherited
typedef MatrixType::RealScalar RealScalar
inherited
typedef MatrixType::Scalar Scalar
inherited
typedef Matrix<Scalar,Dynamic,1> Vector
inherited

Member Enumeration Documentation

anonymous enum
Enumerator:
UpLo 

Constructor & Destructor Documentation

PastixLLT ( )
inline
PastixLLT ( 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
inlineinherited
void analyzePattern ( ColSpMatrix mat)
protectedinherited
void analyzePattern ( const MatrixType matrix)
inline

Compute the LL^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(), and PastixLLT< _MatrixType, _UpLo >::grabMatrix().

void clean ( )
inlineprotectedinherited
Index cols ( void  ) const
inlineinherited
void compute ( ColSpMatrix mat)
protectedinherited
void compute ( const MatrixType matrix)
inline

Compute the L factor of the LL^T supernodal factorization of matrix

See also:
analyzePattern() factorize()

References PastixBase< Derived >::compute(), and PastixLLT< _MatrixType, _UpLo >::grabMatrix().

Referenced by PastixLLT< _MatrixType, _UpLo >::PastixLLT().

PastixLLT< _MatrixType, _UpLo > & derived ( )
inlineinherited
const PastixLLT< _MatrixType, _UpLo > & derived ( ) const
inlineinherited
Array<RealScalar,IPARM_SIZE,1>& dparm ( )
inlineinherited

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)
inlineinherited

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

See also:
dparm()
void factorize ( ColSpMatrix mat)
protectedinherited
void factorize ( const MatrixType matrix)
inline

Compute the LL^T supernodal numerical factorization of matrix

See also:
analyzePattern()

References PastixBase< Derived >::factorize(), and PastixLLT< _MatrixType, _UpLo >::grabMatrix().

void grabMatrix ( const MatrixType matrix,
ColSpMatrix out 
)
inlineprotected
ComputationInfo info ( ) const
inlineinherited

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()
void init ( )
inlineprotected
Array<Index,IPARM_SIZE,1>& iparm ( )
inlineinherited

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)
inlineinherited

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

See also:
iparm()
Index rows ( void  ) const
inlineinherited
const internal::solve_retval<PastixBase, Rhs> solve ( const MatrixBase< Rhs > &  b) const
inlineinherited
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
inlineinherited
Returns:
the solution x of $ A x = b $ using the current decomposition of A.
See also:
compute()

Member Data Documentation

int m_analysisIsOk
protectedinherited
int m_comm
mutableprotectedinherited
Matrix<double,DPARM_SIZE,1> m_dparm
mutableprotectedinherited
int m_factorizationIsOk
protectedinherited
ComputationInfo m_info
mutableprotectedinherited
int m_initisOk
protectedinherited
Matrix<Index,Dynamic,1> m_invp
mutableprotectedinherited
Matrix<int,IPARM_SIZE,1> m_iparm
mutableprotectedinherited
bool m_isInitialized
protectedinherited
pastix_data_t* m_pastixdata
mutableprotectedinherited
Matrix<Index,Dynamic,1> m_perm
mutableprotectedinherited
int m_size
mutableprotectedinherited

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