Polynomial spline spaces of non-uniform bi-degree on T-meshes: combinatorial bounds on the dimension

Polynomial splines are ubiquitous in the fields of computer-aided geometric design and computational analysis. Splines on T-meshes, especially, have the potential to be incredibly versatile since local mesh adaptivity enables efficient modeling and approximation of local features. Meaningful use of such splines for modeling and approximation requires the construction of a suitable spanning set of linearly independent splines, and a theoretical understanding of the spline space dimension can be a useful tool when assessing possible approaches for building such splines. Here, we provide such a tool. Focusing on T-meshes, we study the dimension of the space of bivariate polynomial splines, and we discuss the general setting where local mesh adaptivity is combined with local polynomial degree adaptivity. The latter allows for the flexibility of choosing non-uniform bi-degrees for the splines, i.e., different bi-degrees on different faces of the T-mesh. In particular, approaching the problem using tools from homological algebra, we generalize the framework and the discourse presented by Mourrain (Math. Comput. 83(286):847–871, 2014) for uniform bi-degree splines. We derive combinatorial lower and upper bounds on the spline space dimension and subsequently outline sufficient conditions for the bounds to coincide.

HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB

This paper presents , an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in . Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator is provided to facilitate its application to practical engineering problems.