A Method for Global Optimization of Large Systems of Quadratic Constraints