fix #112. MIPVerify uses the bounds computed from the input set to co…

…ntraint the network. So we also need to add the input set as a constraint to the optimization problem.

src/optimization/mipVerify.jl

@@ -34,6 +34,7 @@ function solve(solver::MIPVerify, problem::Problem)
     model = Model(solver)
     neurons = init_neurons(model, problem.network)
     deltas = init_deltas(model, problem.network)
+    add_set_constraint!(model, problem.input, first(neurons))
     add_complementary_set_constraint!(model, problem.output, last(neurons))
     bounds = get_bounds(problem)
     encode_network!(model, problem.network, neurons, deltas, bounds, BoundedMixedIntegerLP())