A $C^0$ weak Galerkin method with preconditioning for constrained optimal control problems with general tracking

This paper presents a $C^0$ weak Galerkin ($C^0$-WG) method combined with an additive Schwarz preconditioner for solving optimal control problems (OCPs) governed by partial differential equations with general tracking cost functionals and pointwise state constraints. These problems pose significant analytical and numerical challenges due to the presence of fourth-order variational inequalities and the reduced regularity of solutions. Our first contribution is the design of a $C^0$-WG method based on globally continuous quadratic Lagrange elements, enabling efficient elementwise stiffness matrix assembly and parameter-free implementation while maintaining accuracy, as supported by a rigorous error analysis. As a second contribution, we develop an additive Schwarz preconditioner tailored to the $C^0$-WG method to improve solver performance for the resulting ill-conditioned linear systems. Numerical experiments confirm the effectiveness and robustness of the proposed method and preconditioner for both biharmonic and optimal control problems.