Skip to content

Draft: Lobatto shape functions

Simon Praetorius requested to merge feature/lobatto-shape-functions into master


An (adaptive) Lobatto basis is a set of basis functions of local Lobatto shape functions that allows to have different polynomial orders on each element. It even allows to have different orders on sub-entities and thus opens the possibility for conforming bases with p-adaptivity.

Note, for a stable discretization you need to fulfill a minimum-rule that essentially means that the polynomial degree on sub-entities shared by some entities is a most the minimum of the polynomial degrees of the entities.

The implementation is based on the local Lobatto[Cube]LocalFiniteElement implementing the functions from

"Higher-Order Finite Element Methods", P. Soling, K, Segeth, I. Dolezel, 2004, Chapman & Hall/CRC

Edited by Simon Praetorius

Merge request reports