@@ -13,13 +13,13 @@ function perf_lucompletepivCopy!(A)
1313 μ, λ = _maxind (As)
1414 μ += k- 1 ; λ += k- 1
1515 rowpiv[k] = μ
16- A[[k,μ], 1 : n] . = A[[μ,k], 1 : n]
16+ A[[k,μ], 1 : n] = A[[μ,k], 1 : n]
1717 colpiv[k] = λ
18- A[1 : n, [k,λ]] . = A[1 : n, [λ,k]]
18+ A[1 : n, [k,λ]] = A[1 : n, [λ,k]]
1919 if A[k,k] ≠ 0
2020 ρ = k+ 1 : n
21- A[ρ, k] ./= A[k, k]
22- A[ρ, ρ] .- = A[ρ, k: k] . * A[k: k, ρ]
21+ A[ρ, k] = A[ρ, k] / A[k, k]
22+ A[ρ, ρ] = A[ρ, ρ] - A[ρ, k: k] * A[k: k, ρ]
2323 end
2424 end
2525 return (A, rowpiv, colpiv)
@@ -34,13 +34,13 @@ function perf_lucompletepivSub!(A)
3434 μ, λ = _maxind (As)
3535 μ += k- 1 ; λ += k- 1
3636 rowpiv[k] = μ
37- A[[k,μ], 1 : n] . = view (A, [μ,k], 1 : n)
37+ A[[k,μ], 1 : n] = view (A, [μ,k], 1 : n)
3838 colpiv[k] = λ
39- A[1 : n, [k,λ]] . = view (A, 1 : n, [λ,k])
39+ A[1 : n, [k,λ]] = view (A, 1 : n, [λ,k])
4040 if A[k,k] ≠ 0
4141 ρ = k+ 1 : n
42- A[ρ, k] ./= A[k, k]
43- A[ρ, ρ] .- = view (A, ρ, k: k) . * view (A, k: k, ρ)
42+ A[ρ, k] = view (A, ρ, k) / A[k, k]
43+ A[ρ, ρ] = view (A, ρ, ρ) - view (A, ρ, k: k) * view (A, k: k, ρ)
4444 end
4545 end
4646 return (A, rowpiv, colpiv)
0 commit comments