“…Indeed, multivariate B-splines and NURBS come with a rigid tensor-product structure, therefore precluding the simple modeling of local behaviors. More advanced splines, such as locally refined splines [22], hierarchical splines [33,23], T-splines and hierarchical T-splines [65,24], have emerged to offer local mesh refinement, but the integration of arbitrary local models (e.g., inclusions) within spline patches (e.g., matrix) still remains an issue [63,10,69,68,45]. From a general point of view, mainly two approaches can be followed.…”