We present an easy treatment of (multi-level) hanging nodes in hp-finite elements in two dimensions. Its simplicity is due to the fact that the connectivity matrices can directly be computed in a row-wise fashion. This is achieved by treating a global degree of freedom as a local degree of freedom on each element in an appropriate union of elements. This forms a type of support which is larger than the one used in the conventional approach and generalizes the recently presented multi-level hp-method. Numerical experiments compare the new approach with respect to accuracy and the conditioning of the resulting linear system. We conclude that this new approach achieves nearly the same properties as the conventional approach in all investigated examples and provides an excellent trade-off between implementational complexity and solvability.
«