Eliminated conditional constraint matrix since we will not be able to have Dirichlet boundaries with the skew symmetrized stability term version of the problem (only because deal.ii can not support it).