Is there a way to query if a solver supports indicator constraints? I’d like to write a function that produces a model that uses indicator constraints if they are supported and reverts to a big-M constraints if not given an optimizer.
If you try from JuMP, with SCIP solver, you can do:
using JuMP, SCIP
const MOI = MathOptInterface
model = Model(SCIP.Optimizer)