From d38598431ac0a62a8bbe65ac72f22b46eed1743c Mon Sep 17 00:00:00 2001 From: Daisuke Oyama Date: Sat, 18 Jan 2025 13:25:59 +0900 Subject: [PATCH] Translate 0-based indices to Julia's 1-based in removevredundancy! --- src/polyhedron.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/polyhedron.jl b/src/polyhedron.jl index effd529..707d286 100644 --- a/src/polyhedron.jl +++ b/src/polyhedron.jl @@ -185,7 +185,7 @@ function Polyhedra.removevredundancy!(p::Polyhedron) detectvlinearity!(p) ext = getext(p) extm = getextm(p, :AlmostFresh) # FIXME does it need to be fresh ? - redset = redund(extm) + redset = BitSet(redund(extm) .+ 1) nonred = setdiff(BitSet(1:size(ext.R, 1)), redset) nonred = collect(setdiff(nonred, ext.linset)) lin = collect(ext.linset)