On the Approximate Solution of a Class of Large Discrete Quadratic Programming Problems by Modulation: The Case of Circulant Quadratic Forms