A dual gradient-projection method for large-scale strictly convex quadratic problems