Templated class to perform the Cholesky decomposition on a symmetric positive-definite matrix. More...
Public Member Functions | |
LLTMatrix () | |
Construct null. More... | |
LLTMatrix (const SquareMatrix< Type > &M) | |
Construct from a square matrix and perform the decomposition. More... | |
void | decompose (const SquareMatrix< Type > &M) |
Perform the Cholesky decomposition of the matrix. More... | |
void | solve (Field< Type > &x, const Field< Type > &source) const |
Solve the linear system with the given source. More... | |
tmp< Field< Type > > | solve (const Field< Type > &source) const |
Solve the linear system with the given source. More... | |
Public Member Functions inherited from SquareMatrix< Type > | |
SquareMatrix () | |
Null constructor. More... | |
SquareMatrix (const label n) | |
Construct given number of rows/columns. More... | |
template<class MatrixType > | |
SquareMatrix (const ConstMatrixBlock< MatrixType > &) | |
Construct from a block of another matrix. More... | |
template<class MatrixType > | |
SquareMatrix (const MatrixBlock< MatrixType > &) | |
Construct from a block of another matrix. More... | |
SquareMatrix (const label n, const zero) | |
Construct given number of rows/columns. More... | |
SquareMatrix (const label m, const label n, const zero) | |
Construct given number of rows and columns (checked to be equal) More... | |
SquareMatrix (const label n, const Identity< Type >) | |
Construct given number of rows/columns. More... | |
SquareMatrix (const label n, const Type &) | |
Construct with given number of rows and rows. More... | |
SquareMatrix (const RectangularMatrix< Type > &) | |
Construct as copy of a RectangularMatrix. More... | |
SquareMatrix (Istream &) | |
Construct from Istream. More... | |
autoPtr< SquareMatrix< Type > > | clone () const |
Clone. More... | |
void | setSize (const label m) |
Resize the matrix preserving the elements. More... | |
void | shallowResize (const label m) |
Resize the matrix without reallocating storage (unsafe) More... | |
void | operator= (const zero) |
Assignment of all elements to zero. More... | |
void | operator= (const Identity< Type >) |
Assignment elements to the. More... | |
Public Member Functions inherited from Matrix< SquareMatrix< Type >, Type > | |
Matrix () | |
Null constructor. More... | |
Matrix (const label m, const label n) | |
Construct given number of rows and columns. More... | |
Matrix (const label m, const label n, const zero) | |
Construct with given number of rows and columns. More... | |
Matrix (const label m, const label n, const Type &) | |
Construct with given number of rows and columns. More... | |
Matrix (const mType &) | |
Copy constructor. More... | |
Matrix (const Matrix< Form2, Type > &) | |
Copy constructor from matrix of a different form. More... | |
Matrix (const ConstMatrixBlock< MatrixType > &) | |
Construct from a block of another matrix. More... | |
Matrix (const MatrixBlock< MatrixType > &) | |
Construct from a block of another matrix. More... | |
Matrix (Istream &) | |
Construct from Istream. More... | |
autoPtr< mType > | clone () const |
Clone. More... | |
~Matrix () | |
Destructor. More... | |
label | m () const |
Return the number of rows. More... | |
label | n () const |
Return the number of columns. More... | |
label | size () const |
Return the number of elements in matrix (m*n) More... | |
const Type * | v () const |
Return element vector of the constant Matrix. More... | |
Type * | v () |
Return element vector of the Matrix. More... | |
ConstMatrixBlock< mType > | block (const label m, const label n, const label mStart, const label nStart) const |
ConstMatrixBlock< mType > | block (const label mStart, const label nStart) const |
MatrixBlock< mType > | block (const label m, const label n, const label mStart, const label nStart) |
MatrixBlock< mType > | block (const label mStart, const label nStart) |
Foam::ConstMatrixBlock< Foam::Matrix< SquareMatrix< Type >, Type > > | block (const label mStart, const label nStart) const |
Foam::MatrixBlock< Foam::Matrix< SquareMatrix< Type >, Type > > | block (const label mStart, const label nStart) |
ConstMatrixBlock< mType > | col (const label m, const label rowStart) const |
ConstMatrixBlock< mType > | col (const label m, const label mStart, const label nStart) const |
MatrixBlock< mType > | col (const label m, const label rowStart) |
MatrixBlock< mType > | col (const label m, const label mStart, const label nStart) |
void | checki (const label i) const |
Check index i is within valid range (0 ... m-1). More... | |
void | checkj (const label j) const |
Check index j is within valid range (0 ... n-1). More... | |
void | clear () |
Clear the Matrix, i.e. set sizes to zero. More... | |
void | transfer (mType &) |
Transfer the contents of the argument Matrix into this Matrix. More... | |
void | setSize (const label m, const label n) |
Resize the matrix preserving the elements. More... | |
void | shallowResize (const label m, const label n) |
Resize the matrix without reallocating storage (unsafe) More... | |
SquareMatrix< Type > | T () const |
Return the transpose of the matrix. More... | |
Type * | operator[] (const label) |
Return subscript-checked row of Matrix. More... | |
const Type * | operator[] (const label) const |
Return subscript-checked row of constant Matrix. More... | |
const Type & | operator() (const label i, const label j) const |
(i, j) const element access operator More... | |
Type & | operator() (const label i, const label j) |
(i, j) element access operator More... | |
void | operator= (const mType &) |
Assignment operator. Takes linear time. More... | |
void | operator= (const ConstMatrixBlock< MatrixType > &) |
Assignment to a block of another matrix. More... | |
void | operator= (const MatrixBlock< MatrixType > &) |
Assignment to a block of another matrix. More... | |
void | operator= (const zero) |
Assignment of all elements to zero. More... | |
void | operator= (const Type &) |
Assignment of all elements to the given value. More... | |
Additional Inherited Members | |
Public Types inherited from Matrix< SquareMatrix< Type >, Type > | |
typedef Matrix< SquareMatrix< Type >, Type > | mType |
Matrix type. More... | |
typedef Type | cmptType |
Component type. More... | |
Static Public Member Functions inherited from Matrix< SquareMatrix< Type >, Type > | |
static const mType & | null () |
Return a null Matrix. More... | |
Templated class to perform the Cholesky decomposition on a symmetric positive-definite matrix.
Member functions are provided to solve linear systems using the LLT decomposition.
Definition at line 54 of file LLTMatrix.H.
LLTMatrix | ( | ) |
Construct null.
Definition at line 31 of file LLTMatrix.C.
LLTMatrix | ( | const SquareMatrix< Type > & | M | ) |
Construct from a square matrix and perform the decomposition.
Definition at line 36 of file LLTMatrix.C.
void decompose | ( | const SquareMatrix< Type > & | M | ) |
Perform the Cholesky decomposition of the matrix.
Definition at line 45 of file LLTMatrix.C.
References Foam::abort(), Foam::FatalError, FatalErrorInFunction, k, Matrix< SquareMatrix< Type >, Type >::m(), M, LLTMatrix< Type >::solve(), Foam::sqrt(), Foam::sum(), and Foam::Zero.
Solve the linear system with the given source.
and returning the solution in the Field argument x. This function may be called with the same field for x and source.
Definition at line 93 of file LLTMatrix.C.
References Matrix< SquareMatrix< Type >, Type >::m(), and Foam::sum().
Referenced by LLTMatrix< Type >::decompose().
Foam::tmp< Foam::Field< Type > > solve | ( | const Field< Type > & | source | ) | const |
Solve the linear system with the given source.
returning the solution
Definition at line 136 of file LLTMatrix.C.
References tmp< T >::ref(), and solve().