On p-Robust Saturation for hp-AFEM

We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK, REFINE, with ESTIMATE being implemented using the $p$-robust equilibrated flux estimator, and MARK being D\"orfler marking. As a refinement strategy we employ $p$-refinement. We investigate the question by which amount the local polynomial degree on any marked patch has to be increase in order to achieve a $p$-independent error reduction. The resulting adaptive method can be turned into an instance optimal $hp$-adaptive method by the addition of a coarsening routine.