Optimal control problems for the heat equation with pointwisebilateral control-state constraints are considered. A locallysuperlinearly convergent numerical solution algorithm is proposedand its mesh independence is established. Further, for theefficient numerical solution reduced space and Schur complementbased preconditioners are proposed which take into account theactive and inactive set structure of the problem. The paper endsby numerical tests illustrating our theoretical findings andcomparing the efficiency of the proposed preconditioners.
