nojhan/eodev
Archived
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

more pertinent names for cholesky factorization functions

Browse source
This commit is contained in:
nojhan 2011-11-12 14:24:54 +01:00
commit d5d8f87a90
1 changed files with 9 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

18
edo/src/edoSamplerNormalMulti.h
View file

@ -61,11 +61,11 @@ public:
typedef ublas::symmetric_matrix< AtomType, ublas::lower > MatrixType; typedef ublas::symmetric_matrix< AtomType, ublas::lower > MatrixType;
enum Method { enum Method {
//! use the standard algorithm, with square root @see factorize //! use the standard algorithm, with square root @see factorize_LLT
standard, standard,
//! use the algorithm using absolute value within the square root @see factorize_abs //! use the algorithm using absolute value within the square root @see factorize_LLT_abs
absolute, absolute,
//! use the robust algorithm, without square root //! use the robust algorithm, without square root @see factorize_LDLT
robust robust
}; };
Method _use; Method _use;
@ -147,11 +147,11 @@ public:
void factorize( const MatrixType& V ) void factorize( const MatrixType& V )
{ {
if( _use == standard ) { if( _use == standard ) {
factorize_std( V ); factorize_LLT( V );
} else if( _use == absolute ) { } else if( _use == absolute ) {
factorize_abs( V ); factorize_LLT_abs( V );
} else if( _use == robust ) { } else if( _use == robust ) {
factorize_robust( V ); factorize_LDLT( V );
} }
} }
@ -162,7 +162,7 @@ public:
* When compiled in debug mode and called on ill-conditionned matrix, * When compiled in debug mode and called on ill-conditionned matrix,
* will raise an assert before calling the square root on a negative number. * will raise an assert before calling the square root on a negative number.
*/ */
void factorize_std( const MatrixType& V) void factorize_LLT( const MatrixType& V)
{ {
unsigned int N = assert_properties( V ); unsigned int N = assert_properties( V );
@ -208,7 +208,7 @@ public:
* Be aware that this increase round-off errors, this is just a ugly * Be aware that this increase round-off errors, this is just a ugly
* hack to avoid crash. * hack to avoid crash.
*/ */
void factorize_abs( const MatrixType & V) void factorize_LLT_abs( const MatrixType & V)
{ {
unsigned int N = assert_properties( V ); unsigned int N = assert_properties( V );
@ -248,7 +248,7 @@ public:
* *
* Computes L and D such as V = L D Lt * Computes L and D such as V = L D Lt
*/ */
void factorize_robust( const MatrixType& V) void factorize_LDLT( const MatrixType& V)
{ {
unsigned int N = assert_properties( V ); unsigned int N = assert_properties( V );