Implement a QkPkDG finite element map
Summary
The finite element maps of PDELab are missing the discontinuous counterpart of PkQkLocalFiniteElementMap
.
Proposal
Create a finite element map that implements discontinuous function spaces on rectangles and simplices, a PkQkDGLocalFiniteElementMap
.
Use case
I solve a DG discretization on a locally adaptive UG grid and use either PkLocal...
or QkDGLocalFiniteElementMap
, depending on the geometry type (rectangles or simplices). For rectangles, UG implements an 'automatic closure' feature to avoid hanging nodes after local refinement. I have to turn this off because it adds simplices to the grid, which QkDG
cannot handle.
...pretty please?
Unfortunately, I'm not confident enough with the PDELab-internals to write one myself. I hope this might be a cakewalk for a developer. Otherwise, you can guide me through the steps for building it myself (and opening an MR)