File tree Expand file tree Collapse file tree 1 file changed +4
-4
lines changed
Expand file tree Collapse file tree 1 file changed +4
-4
lines changed Original file line number Diff line number Diff line change @@ -350,15 +350,15 @@ PMMatrix >> choleskyDecomposition [
350350 1 to: i do: [ :j |
351351 i = j
352352 ifTrue: [
353- rowSum := (1 to: j - 1 ) sum: [ :k |
354- (upperTriangular at: k at: j) squared ].
353+ rowSum := (1 to: j - 1 ) inject: 0 into: [ :sum :k |
354+ sum + (upperTriangular at: k at: j) squared ].
355355
356356 diagonalValue := ((self at: j at: j) - rowSum) sqrt.
357357
358358 upperTriangular at: j at: j put: diagonalValue ]
359359 ifFalse: [
360- partialSum := (1 to: j - 1 ) sum: [ :k |
361- (upperTriangular at: k at: i) * (upperTriangular at: k at: j) ].
360+ partialSum := (1 to: j - 1 ) inject: 0 into: [ :sum :k |
361+ sum + (upperTriangular at: k at: i) * (upperTriangular at: k at: j) ].
362362
363363 factor := upperTriangular at: j at: j.
364364 nonDiagonalValue := ((self at: j at: i) - partialSum) / factor.
You can’t perform that action at this time.
0 commit comments