> On 25 May 2017, at 19:46, Matthew Knepley <knep...@gmail.com> wrote: > > Sounds like DG. I will get out my dead chicken for the incantation
Actually no! Mixed H(div)-L2 for Stokes. Which has facet integrals for partially discontinuous fields. If you do redundant compute for such terms, you need a depth-2 FEM adjacency, which is just grim. Equally we have some strange users who have jump terms in CG formulations. Lawrence