The fast multigrid solution of an optimal control problem governed by a convection–diffusion partial-integro differential equation is investigated. This optimization problem considers a cost functional of tracking type and a constrained distributed control. The optimal control sought is characterized by the solution to the corresponding optimality system, which is approximated by a finite volume and quadrature discretization schemes and solved by multigrid techniques. The proposed multigrid approach combines a multigrid method for the governing model with a fast multigrid integration method. The convergence of this solution procedure is analyzed by local Fourier analysis and validated by results of numerical experiments.