We analyse the performance of direct solvers when applied to a system of linear equations arising from an $h$-adapted, $C^0$ finite element space. Theoretical estimates are derived for typical $h$-refinement patterns arising ...