isdisjoint is not symmetric/commutative

[schillic]

P below is a box-shaped polytope. Q below is a flat polyhedron.

julia> p = HPolytope([HalfSpace([1.0, 0.0], 0.308261),
        HalfSpace([-1.0, 0.0], -0.2),
        HalfSpace([0.0, 1.0], 0.1),
        HalfSpace([0.0, -1.0], 0.106045)]);

julia> q = HPolyhedron([
        HalfSpace([-1.0, 0.0], 0.0),
        HalfSpace([-0.714286, -1.0], 0.0),
        HalfSpace([0.714286, 1.0], 0.0)]);

julia> isdisjoint(p, q)
false

julia> isdisjoint(q, p)
true

(outsourced from #1038)

[mforets]

julia> [is_intersection_empty(q, Hi) for Hi in constraints_list(p)]
4-element Array{Bool,1}:
 false
 false
 false
 false

julia> [is_intersection_empty(p, Hi) for Hi in constraints_list(q)]
3-element Array{Bool,1}:

 false
 false
  true

julia> [-ρ(-Hi.a, q) for Hi in constraints_list(p)]
4-element Array{Float64,1}:
   -0.0
 -Inf
 -Inf
   -0.0

julia> [Hi.b for Hi in constraints_list(p)]
4-element Array{Float64,1}:
  0.308261
 -0.2
  0.1
  0.106045

[mforets]

possibly related to #902 and to #754 (see last comment)