In this work we construct multigrid preconditioners to accelerate the solution process of a linear-quadratic optimal control problem constrained by the Stokes system. The first order optimality conditions of the control problem form a linear system (the KKT system) connecting the state, adjoint, and control variables. Our approach is to eliminate the state and adjoint variables by essentially solving two Stokes systems, and to construct efficient

There has been a long-standing question of whether certain mesh restrictions are required for a maximum condition to hold for the discrete equations arising from a finite element approximation of an elliptic problem. This is related to knowing whether the discrete Green's function is positive for triangular meshes allowing sufficiently good approximation of $H^1$ functions. We study this question for the Poisson problem in two dimensions discretized