model.optimize select a qp solver by default

This is a fix for model.optimize. When no solver is passed in, either
gurobi or cplex is used (if available) instead of glpk if a quadratic
objective is given.