Unexpected paving #149
Description
using IntervalConstraintProgramming
julia> G1 = @constraint(1.0 + x0*t ≥ 1.2, [x0, t]);
julia> G2 = @constraint(1.0 + x0*t ≤ 1.4, [x0, t]);
julia> S = G1 ∩ G2
Separator:
- variables: x0, t
- expression: (1.0 + x0 * t ∈ [1.19999, ∞]) ∩ (1.0 + x0 * t ∈ [-∞, 1.40001])
julia> D = IntervalBox(-1..1, 0..0.4);
julia> inner, outer = S(D)
([0.499999, 1] × [0, 0.400001], [-1, 1] × [0, 0.400001])The result [0, 0.400001] in the "inner" part for t is unexpected. Compare to the following plots that help to check that the valid contraction for t is [0.2, 0.4] by changing the eps in:
plot(pave(S, D, eps).inner, lab="inner")
plot!(pave(S, D, eps).boundary, lab="boundary")The transition between empty inner and non-empty inner happens for eps at <= 0.5.