diff --git a/src/chordal/decomposition/psd_completion.jl b/src/chordal/decomposition/psd_completion.jl index 74d185ac..c629ca2b 100644 --- a/src/chordal/decomposition/psd_completion.jl +++ b/src/chordal/decomposition/psd_completion.jl @@ -70,13 +70,13 @@ function psd_complete!(A::AbstractMatrix{T}, pattern::SparsityPattern) where {T #sorted using numerical ordering σ(i) α = get_separators(sntree, j) + ν = collect(ν) #unborks indexing below + α = collect(α) #unborks indexing below + # index set containing the row indices of the lower-triangular zeros in # column i (i: representative index) sorted by σ(i) i = ν[1] - ν = collect(ν) #unborks indexing below - α = collect(α) #unborks indexing below - η = collect(i+1:1:N) filter!(x -> !insorted(x, α) && !insorted(x, ν), η) @@ -103,4 +103,3 @@ function psd_complete!(A::AbstractMatrix{T}, pattern::SparsityPattern) where {T A[:, :] = W[ip, ip] end -