Call scipy.optimize inside pyomo
At the moment (Dec 2017), there is no built-in support for passing a Pyomo model to scipy.optimize.  That said, it would not be a very difficult task to write a reasonably general purpose object that could generate the necessary (value, Jacobian, Hessian) evaluation functions to pass to scipy.optimize.minimize().
Can I integrate scipy.optimize.minimize solver with method=SLSQP inside pyomo? Modeling in pyomo is much faster than in scipy but pyomo documentation does not seem to say explicitly if this is feasible.