# | 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 |