The covariant formulation with P4estMesh as implemented in #47 does not support MPI nor adaptivity. We will eventually want to add support for these features, but it may be best to wait to do this with t8code, as that will allow for more general element types.