“…In contrast, there are no B-spline based optimally convergent smooth constructions for unstructured hexahedral meshes. Currently, optimality is achieved by either reducing continuity to C 0 around extraordinary features [56][57][58], combining B-splines with meshless approximants [59] or resorting to non-standard spline definitions [60]. The first approach is not suitable for higher-order partial differential equations whereas the latter approaches lead to schemes that are usually computationally costly.…”