One stroke complementarity for poisson-like problems

Taking electrokinetics as a paradigm problem for the sake of simplicity, complementarity originates when an irrotational electric field and a solenoidal current density satisfying boundary conditions are in hand. We first compare three formulations to obtain a solenoidal current density, both in terms of pure computational advantage and in the ability to pursue symmetric energy bounds with respect to the standard electric scalar potential formulation. For these formulations, we devise post-processing techniques that promise to provide bilateral bounds in one stroke, hence requiring the solution of just one linear system.