(function from pyomo.contrib.gdpopt.create_oa_subproblems)
pyomo.contrib.gdpopt.create_oa_subproblems
Calls the specified transformation (by default bigm) on the original model and removes nonlinear constraints to create a MILP discrete problem.