SOFA API  c76874b7
Open source framework for multi-physics simuation
sofa::component::linearsolver::CholeskySolver< TMatrix, TVector > Class Template Reference

#include <CholeskySolver.h>

Direct linear solver based on Cholesky factorization, for dense matrices. More...

Inheritance diagram for sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >:

Detailed Description

template<class TMatrix, class TVector>
class sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >

Direct linear solver based on Cholesky factorization, for dense matrices.

Public Attributes

Data< boolf_verbose
 Dump system state at each iteration. More...
 

Public Member Functions

 SOFA_CLASS (SOFA_TEMPLATE2(CholeskySolver, TMatrix, TVector), SOFA_TEMPLATE2(sofa::component::linearsolver::MatrixLinearSolver, TMatrix, TVector))
 
 CholeskySolver ()
 
void solve (Matrix &M, Vector &x, Vector &b) override
 Compute x such as Mx=b. M is not used, it must have been factored before using method invert(Matrix& M) More...
 
void invert (Matrix &M) override
 Factors the matrix. Must be done before solving. More...
 

Attribute details

template<class TMatrix , class TVector >
Data<bool> sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >::f_verbose

Dump system state at each iteration.

Constructor details

template<class TMatrix , class TVector >
sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >::CholeskySolver ( )

Function details

template<class TMatrix , class TVector >
void sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >::invert ( Matrix M)
override

Factors the matrix. Must be done before solving.

template<class TMatrix , class TVector >
sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >::SOFA_CLASS ( SOFA_TEMPLATE2(CholeskySolver< TMatrix, TVector >, TMatrix, TVector)  ,
SOFA_TEMPLATE2(sofa::component::linearsolver::MatrixLinearSolver, TMatrix, TVector)   
)
template<class TMatrix , class TVector >
void sofa::component::linearsolver::CholeskySolver< TMatrix, TVector >::solve ( Matrix M,
Vector z,
Vector r 
)
override

Compute x such as Mx=b. M is not used, it must have been factored before using method invert(Matrix& M)

Factorisation : A = LL^t A x = b <=> LL^t x = b <=> L u = b , L^t x = u