Quasi-optimal convergence rate for an adaptive finite element method

We analyze the simplest and most standard adaptive finite element method (AFEM), with any polynomial degree, for general second order linear, symmetric elliptic operators. As is customary in practice, the AFEM marks exclusively according to the error estimator and performs a minimal element refinement without the interior node property. We prove that the AFEM is a contraction, for the sum of the energy error and the scaled error estimator, between two consecutive adaptive loops. This geometric decay is instrumental to derive the optimal cardinality of the AFEM. We show that the AFEM yields a decay rate of the energy error plus oscillation in terms of the number of degrees of freedom as dictated by the best approximation for this combined nonlinear quantity.