On the Solution of Large Quadratic Programming Problems with Bound Constraints