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

typedef AtomType, in case it would be useful

Browse source
This commit is contained in:
nojhan 2012-01-29 18:19:07 +01:00
commit cfd0579a4a
1 changed files with 7 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

11
cholesky.h
View file

@ -55,6 +55,9 @@ public:
// FIXME check if triangular types behaviour is like having 0
typedef ublas::matrix< T > FactorMat;
//! The type of the scalars
typedef T AtomType;
/** Instanciate without computing anything, you are responsible of
* calling the algorithm and getting the result with operator()
* */
@ -322,8 +325,8 @@ public:
T sum = 0.0;
// for( unsigned int k=1; k<i-1; ++k) {
// for( unsigned int k=1; k<i; ++k) {
for( unsigned int k=0; k<i; ++k) {
// for( unsigned int k=0; k<i-1; ++k) {
// for( unsigned int k=0; k<i; ++k) {
for( unsigned int k=0; k<i-1; ++k) {
sum += this->_L(i, k) * this->_L(i, k);
}
this->_L(i,i) = this->_L(i,i) - sum;
@ -338,8 +341,8 @@ public:
sum = 0.0;
// for( unsigned int k = 1; k < i-1; ++k ) {
// for( unsigned int k = 1; k < i; ++k ) {
for( unsigned int k = 0; k < i; ++k ) {
// for( unsigned int k = 0; k < i-1; ++k ) {
// for( unsigned int k = 0; k < i; ++k ) {
for( unsigned int k = 0; k < i-1; ++k ) {
sum += this->_L(j, k) * this->_L(i, k);
}
this->_L(j,i) = this->_L(j,i) - sum;