Backed out changeset: 159f9e7014f8
This commit is contained in:
@@ -126,7 +126,7 @@ bool MultiDiagonalSymmetricMatrix::CreateDiagonal(unsigned int index, unsigned i
|
||||
return true;
|
||||
}
|
||||
|
||||
int MultiDiagonalSymmetricMatrix::FindIndex(unsigned int StartRow) const {
|
||||
int MultiDiagonalSymmetricMatrix::FindIndex(unsigned int StartRow){
|
||||
//There's GOT to be a better way to do this. "Bidirectional map?"
|
||||
for(unsigned int i = 0; i != m; i++)
|
||||
if(StartRows[i] == StartRow)
|
||||
@@ -274,7 +274,7 @@ void MultiDiagonalSymmetricMatrix::KillIncompleteCholeskyFactorization(void){
|
||||
delete IncompleteCholeskyFactorization;
|
||||
}
|
||||
|
||||
void MultiDiagonalSymmetricMatrix::CholeskyBackSolve(float *x, float *b) const {
|
||||
void MultiDiagonalSymmetricMatrix::CholeskyBackSolve(float *x, float *b){
|
||||
//We want to solve L D Lt x = b where D is a diagonal matrix described by Diagonals[0] and L is a unit lower triagular matrix described by the rest of the diagonals.
|
||||
//Let D Lt x = y. Then, first solve L y = b.
|
||||
float *y = new float[n];
|
||||
|
Reference in New Issue
Block a user