A parallel hp-adaptive finite element method

This paper describes a combination of automatic hp-adaptive finite elements and domain decomposition. The combination is based on the Bank-Hoist parallel adaptive meshing paradigm. The hp-adaptivity is based on a derivative recovery technique while the domain decomposition method is formulated using a mortar-like formulation with Dirac delta functions for the mortar element space. Numerical results show that the approach can scale to hundreds of processors and the convergence of the domain decomposition method is independent of the number of processors as well as the distribution of element sizes and element degrees among the local meshes.