In this note we present a multigrid preconditioning method for solving quadratic optimization problems constrained by fractional diffusion equation. Multigrid methods within the all-at-once approach to solve first order optimality Karush–Kuhn–Tucker (KKT) systems are widely popular, but their development have relied on underlying being sparse. On other hand, most discretizations, matrix represe...