| # | Line 47 | Line 47 | int CholeskyDecomposition(MatrixType& A, MatrixType& L | |
|---|---|---|
| 47 | int n = A.getNRow(); | |
| 48 | assert(n == A.getNCol() && n == L.getNRow()&& n==L.getNCol()); | |
| 49 | for(int i = 0; i < n; ++i) { | |
| 50 | < | double sum1 = 0; |
| 50 | > | RealType sum1 = 0; |
| 51 | for (int k = 0; k < i -1; ++k) { | |
| 52 | sum1 +=L(i,k)*L(i,k); | |
| 53 | } | |
| 54 | L(i, i) = sqrt(A(i, i) - sum1); | |
| 55 | for (int j = i+1; j < n; ++j) { | |
| 56 | < | double sum2 = 0; |
| 56 | > | RealType sum2 = 0; |
| 57 | for (int k = 0; k < i-1; ++k) { | |
| 58 | sum2 += L(j ,k)*L(i, k); | |
| 59 | } | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |