| 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 |
|
} |
| 60 |
< |
A(j, i) = (A(j, i) - sum2) /L(i,i); |
| 60 |
> |
L(j, i) = (A(j, i) - sum2) /L(i,i); |
| 61 |
|
} |
| 62 |
|
} |
| 63 |
|
|
