Benoît Legat

The `hchebyshevcenter` needs the H-representation so it triggers `computehrep!` which calls the QHull library which seems to be what is segfaulting. Could you try replacing `hchebyshevcenter(my_polyhedron, lp_solver)` by `computehrep!(my_polyhedron)` and...

Arf, this is an issue with QHull.jl (opened in https://github.com/JuliaPolyhedra/QHull.jl/issues/33), replace it with `hrep(my_polyhedron)` then

That's annoying indeed, we could provide a dictionary to help this. The order is the creation order which is the same as their order in `all_variables(model)`

We should probably expose the `MOIU.IndexMap` returned here: https://github.com/JuliaPolyhedra/Polyhedra.jl/blob/8a0d39a1bc8610530e999b29d8acdfa0db17ead8/src/lphrep.jl#L21 We could have a ```julia poly, ref_map = polyhedron_with_map(model, lib) ``` where `ref_map` maps `JuMP.VariableRef` to the index between 1 and...

Indeed, it works with CDD, I am leaving the issue opened as it should be fixed with `DefaultLibrary`.

That should be possible indeed. There is in fact even a quite elegant way to do this now with https://github.com/jump-dev/MutableArithmetics.jl. I'l try to do it today.

Not yet. We need to figure out what it would be. It could be a struct providing a mapping between H-Rep indices and the dual of the corresponding constraint

> The polyhedra shouldn't have white spaces between their neighbor polyhedron Not sure what you mean mean by that. You can try replacing `polyhedron(v)` by `polyhedron(v, CDDLib.Library(:float))` (adding `using CDDLib`...

Which version of `CDDLib`, `Polyhedra` and Julia are you using ? Can you try `CDDLib.Library()` in the REPL to isolate the issue ?

Any update ?