Pass correct flatten flag on facet integrals
Add mixed Poisson test with strong boundary conditions