“…Preconditioning PDE-constrained optimization problems has been a subject of considerable interest of late (see, e.g., [53,50,60,44,42,43,45,40]), however relatively few such methods have made use of the double saddle-point structure. We refer the reader to recent, related work [37,55,41,8]. We provide theoretical results on preconditioners for double saddle-point systems and, having arranged the linear systems obtained from the active-set approach to this form, we present efficient and flexible approximations to be applied within a Krylov-subspace solver.…”