Optimization problems constrained by time-dependent Partial Differential Equations (PDEs) are challenging from a computational point of view: even in the simplest case, one needs to solve a system of PDEs coupled globally ...