I have a Quadratic Programming problem with linear constraints. My objective is Quadratic-Convex, the constraint matrix is Totally Unimodular (similar to assignment or network flow problems), and the variables are Binary. I know that if an LP has Totally Unimodular constraint matrix $A$ and integral $b$ matrices, then the Optima will be integral.
Now, while solving my QP with a nonlinear solver (MINOS 5.51) with relaxed variable (ignoring the integrality), I get almost all the solutions as integral. Most of the time, I get all variables as integral. Only 2 in 20,000 variables or 2 in 84 thousands variables are non-integral.
My question is, is there any theory which talks about the implication of Total Unimodularity in QP or MIQP or convex optimization? Is there any explanation why my problem is getting an integral solution? Are there any articles or books which discuss this issue?
Let me assume that the feasible region is full dimensional. There are two issues:
None of the corner points of the feasible region might be optimal (since the objective function is not linear).
A QP is typically solved with a different algorithm (interior point vs [simplex method or interior point followed by the simplex method]). Therefore, if an entire face of the feasible region is optimal (which is possible when the objective function is not strictly convex), the optimization algorithm will find a non-corner point solution.