“…It was concluded from these comparative studies that high-order polynomial methods in combination with static condensation are able to deliver comparable or, in some cases, even superior performance. On the other hand, high order hp-discontinuous Galerkin methods have been demonstrated as competitive discretisation schemes (see, e.g., [43,44] and the references cited therein), due to their flexibility in dealing with hp-adaptivity and general shaped elements, and their ability in achieving the exponential convergence rate of spectral techniques.…”