On the Use of Elliptic Regularity Theory for the Numerical Solution of Variational Problems