Skip to content

Commit 7828019

Browse files
YingboMashashi
andcommitted
Fix sym_lu
Co-authored-by: "Shashi Gowda" <[email protected]>
1 parent 510f597 commit 7828019

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

src/solve.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -19,7 +19,7 @@ function sym_lu(A)
1919
U = copy!(Array{Any}(undef, size(A)),A)
2020
p = BlasInt[1:m;]
2121
for k = 1:m-1
22-
_, i = findmin(map(ii->iszero(U[ii, k]) ? Inf : nterms(U[ii,k]), k:n))
22+
_, i = findmin(map(ii->_iszero(U[ii, k]) ? Inf : nterms(U[ii,k]), k:n))
2323
i += k - 1
2424
# swap
2525
U[k, k:end], U[i, k:end] = U[i, k:end], U[k, k:end]
@@ -28,7 +28,7 @@ function sym_lu(A)
2828

2929
for j = k+1:m
3030
L[j,k] = U[j, k] / U[k, k]
31-
U[j,k:m] .= U[j,k:m] .- L[j,k] .* U[k,k:m]
31+
U[j,k:m] .= U[j,k:m] .- Ref(L[j,k]) .* U[k,k:m]
3232
end
3333
end
3434
for j=1:m

0 commit comments

Comments
 (0)