Search

The space of $C^1$ cubic Clough-Tocher splines is a classical finite element approximation space over triangulations for solving partial differential equations. However, for such a space there is no B-spline basis available, which is a preferred choice in computer aided geometric design and isogeometric analysis. A B-spline basis is a locally supported basis that forms a convex partition of unity. In this paper, we explore several alternative $C^1$ cubic spline spaces over triangulations equipped with a B-spline basis. They are defined over a Powell-Sabin refined triangulation and present different types of $C^2$ super-smoothness. The super-smooth B-splines are obtained through an extraction process, i.e., they are expressed in terms of less smooth basis functions. These alternative spline spaces maintain the same optimal approximation power as Clough-Tocher splines. This is illustrated with a selection of numerical examples in the context of least squares approximation and finite element approximation for second and fourth order boundary value problems.

Tchebycheffian splines are smooth piecewise functions whose pieces are drawn from (possibly different) Tchebycheff spaces, a natural generalization of algebraic polynomial spaces. They enjoy most of the properties known in the polynomial spline case. In particular, under suitable assumptions, Tchebycheffian splines admit a representation in terms of basis functions, called Tchebycheffian B-splines (TB-splines), completely analogous to polynomial B-splines. A particularly interesting subclass consists of Tchebycheffian splines with pieces belonging to null-spaces of constant-coefficient linear differential operators. They grant the freedom of combining polynomials with exponential and trigonometric functions with any number of individual shape parameters. Moreover, they have been recently equipped with efficient evaluation and manipulation procedures. In this paper, we consider the use of TB-splines with pieces belonging to null-spaces of constant-coefficient linear differential operators as an attractive substitute for standard polynomial B-splines and rational NURBS in isogeometric Galerkin methods. We discuss how to exploit the large flexibility of the geometrical and analytical features of the underlying Tchebycheff spaces according to problem-driven selection strategies. TB-splines offer a wide and robust environment for the isogeometric paradigm beyond the limits of the rational NURBS model., Comment: 35 pages, 18 figures

In this paper we present ExSpliNet, an interpretable and expressive neural network model. The model combines ideas of Kolmogorov neural networks, ensembles of probabilistic trees, and multivariate B-spline representations. We give a probabilistic interpretation of the model and show its universal approximation properties. We also discuss how it can be efficiently encoded by exploiting B-spline properties. Finally, we test the effectiveness of the proposed model on synthetic approximation problems and classical machine learning benchmark datasets.

In this paper, we address the problem of constructing $C^2$ cubic spline functions on a given arbitrary triangulation $\mathcal{T}$. To this end, we endow every triangle of $\mathcal{T}$ with a Wang-Shi macro-structure. The $C^2$ cubic space on such a refined triangulation has a stable dimension and optimal approximation power. Moreover, any spline function in such space can be locally built on each of the macro-triangles independently via Hermite interpolation. We provide a simplex spline basis for the space of $C^2$ cubics defined on a single macro-triangle which behaves like a Bernstein/B-spline basis over the triangle. The basis functions inherit recurrence relations and differentiation formulas from the simplex spline construction, they form a nonnegative partition of unity, they admit simple conditions for $C^2$ joins across the edges of neighboring triangles, and they enjoy a Marsden-like identity. Also, there is a single control net to facilitate control and early visualization of a spline function over the macro-triangle. Thanks to these properties, the complex geometry of the Wang-Shi macro-structure is transparent to the user. Stable global bases for the full space of $C^2$ cubics on the Wang-Shi refined triangulation $\mathcal{T}$ are deduced from the local simplex spline basis by extending the concept of minimal determining sets.

In this work, we focus on a fractional differential equation in Riesz form discretized by a polynomial B-spline collocation method. For an arbitrary polynomial degree $p$, we show that the resulting coefficient matrices possess a Toeplitz-like structure. We investigate their spectral properties via their symbol and we prove that, like for second order differential problems, also in this case the given matrices are ill-conditioned both in the low and high frequencies for large $p$. More precisely, in the fractional scenario the symbol has a single zero at $0$ of order $\alpha$, with $\alpha$ the fractional derivative order that ranges from $1$ to $2$, and it presents an exponential decay to zero at $\pi$ for increasing $p$ that becomes faster as $\alpha$ approaches $1$. This translates in a mitigated conditioning in the low frequencies and in a deterioration in the high frequencies when compared to second order problems. Furthermore, the derivation of the symbol reveals another similarity of our problem with a classical diffusion problem. Since the entries of the coefficient matrices are defined as evaluations of fractional derivatives of the B-spline basis at the collocation points, we are able to express the central entries of the coefficient matrix as inner products of two fractional derivatives of cardinal B-splines. Finally, we perform a numerical study of the approximation behavior of polynomial B-spline collocation. This study suggests that, in line with non-fractional diffusion problems, the approximation order for smooth solutions in the fractional case is $p+2-\alpha$ for even $p$, and $p+1-\alpha$ for odd $p$., Comment: 23 pages, 13 figures

We show that isogeometric Galerkin discretizations of eigenvalue problems related to the Laplace operator subject to any standard type of homogeneous boundary conditions have no outliers in certain optimal spline subspaces. Roughly speaking, these optimal subspaces are obtained from the full spline space defined on certain uniform knot sequences by imposing specific additional boundary conditions. The spline subspaces of interest have been introduced in the literature some years ago when proving their optimality with respect to Kolmogorov $n$-widths in $L^2$-norm for some function classes. The eigenfunctions of the Laplacian -- with any standard type of homogeneous boundary conditions -- belong to such classes. Here we complete the analysis of the approximation properties of these optimal spline subspaces. In particular, we provide explicit $L^2$ and $H^1$ error estimates with full approximation order for Ritz projectors in the univariate and in the multivariate tensor-product setting. Besides their intrinsic interest, these estimates imply that, for a fixed number of degrees of freedom, all the eigenfunctions and the corresponding eigenvalues are well approximated, without loss of accuracy in the whole spectrum when compared to the full spline space. Moreover, there are no spurious values in the approximated spectrum. In other words, the considered subspaces provide accurate outlier-free discretizations in the univariate and in the multivariate tensor-product case. This main contribution is complemented by an explicit construction of B-spline-like bases for the considered spline subspaces. The role of such spaces as accurate discretization spaces for addressing general problems with non-homogeneous boundary behavior is discussed as well., Comment: 44 pages, 19 figures. To appear in CMAME

In this paper we construct Ritz-type projectors with boundary interpolation properties in finite dimensional subspaces of the usual Sobolev space and we provide a priori error estimates for them. The abstract analysis is exemplified by considering spline spaces and we equip the corresponding error estimates with explicit constants. This complements our results recently obtained for explicit spline error estimates based on the classical Ritz projectors in [Numer. Math. 144(4):889--929, 2020]., Comment: 19 pages, 3 figures. Improved the presentation

Multi-degree Tchebycheffian splines are splines with pieces drawn from extended (complete) Tchebycheff spaces, which may differ from interval to interval, and possibly of different dimensions. These are a natural extension of multi-degree polynomial splines. Under quite mild assumptions, they can be represented in terms of a so-called MDTB-spline basis; such basis possesses all the characterizing properties of the classical polynomial B-spline basis. We present a practical framework to compute MDTB-splines, and provide an object-oriented implementation in Matlab. The implementation supports the construction, differentiation, and visualization of MDTB-splines whose pieces belong to Tchebycheff spaces that are null-spaces of constant-coefficient linear differential operators. The construction relies on an extraction operator that maps local Tchebycheffian Bernstein functions to the MDTB-spline basis of interest.

In this paper, we describe a general class of $C^1$ smooth rational splines that enables, in particular, exact descriptions of ellipses and ellipsoids - some of the most important primitives for CAD and CAE. The univariate rational splines are assembled by transforming multiple sets of NURBS basis functions via so-called design-through-analysis compatible extraction matrices; different sets of NURBS are allowed to have different polynomial degrees and weight functions. Tensor products of the univariate splines yield multivariate splines. In the bivariate setting, we describe how similar design-through-analysis compatible transformations of the tensor-product splines enable the construction of smooth surfaces containing one or two polar singularities. The material is self-contained, and is presented such that all tools can be easily implemented by CAD or CAE practitioners within existing software that support NURBS. To this end, we explicitly present the matrices (a) that describe our splines in terms of NURBS, and (b) that help refine the splines by performing (local) degree elevation and knot insertion. Finally, all $C^1$ spline constructions yield spline basis functions that are locally supported and form a convex partition of unity., Comment: 16 pages

We relate the problem of best low-rank approximation in the spectral norm for a matrix $A$ to Kolmogorov $n$-widths and corresponding optimal spaces. We characterize all the optimal spaces for the image of the Euclidean unit ball under $A$ and we show that any orthonormal basis in an $n$-dimensional optimal space generates a best rank-$n$ approximation to $A$. We also present a simple and explicit construction to obtain a sequence of optimal $n$-dimensional spaces once an initial optimal space is known. This results in a variety of solutions to the best low-rank approximation problem and provides alternatives to the truncated singular value decomposition. This variety can be exploited to obtain best low-rank approximations with problem-oriented properties., Comment: 26 pages, 1 figure. Article published in SIMAX

In this paper we describe an adaptive refinement strategy for LR B-splines. The presented strategy ensures, at each iteration, local linear independence of the obtained set of LR B-splines. This property is then exploited in two applications: the construction of efficient quasi-interpolation schemes and the numerical solution of elliptic problems using the isogeometric Galerkin method., Comment: 23 pages, 14 figures

In this paper we present an efficient and robust approach to compute a normalized B-spline-like basis for spline spaces with pieces drawn from extended Tchebycheff spaces. The extended Tchebycheff spaces and their dimensions are allowed to change from interval to interval. The approach works by constructing a matrix that maps a generalized Bernstein-like basis to the B-spline-like basis of interest. The B-spline-like basis shares many characterizing properties with classical univariate B-splines and may easily be incorporated in existing spline codes. This may contribute to the full exploitation of Tchebycheffian splines in applications, freeing them from the restricted role of an elegant theoretical extension of polynomial splines. Numerical examples are provided that illustrate the procedure described.

In this paper we provide a priori error estimates with explicit constants for both the $L^2$-projection and the Ritz projection onto spline spaces of arbitrary smoothness defined on arbitrary grids. This extends the results recently obtained for spline spaces of maximal smoothness. The presented error estimates are in agreement with the numerical evidence found in the literature that smoother spline spaces exhibit a better approximation behavior per degree of freedom, even for low smoothness of the functions to be approximated. First we introduce results for univariate spline spaces, and then we address multivariate tensor-product spline spaces and isogeometric spline spaces generated by means of a mapped geometry, both in the single-patch and in the multi-patch case., Comment: 39 pages, 4 figures. Improved the presentation. Article published in Numerische Mathematik

In this paper we provide a priori error estimates in standard Sobolev (semi-)norms for approximation in spline spaces of maximal smoothness on arbitrary grids. The error estimates are expressed in terms of a power of the maximal grid spacing, an appropriate derivative of the function to be approximated, and an explicit constant which is, in many cases, sharp. Some of these error estimates also hold in proper spline subspaces, which additionally enjoy inverse inequalities. Furthermore, we address spline approximation of eigenfunctions of a large class of differential operators, with a particular focus on the special case of periodic splines. The results of this paper can be used to theoretically explain the benefits of spline approximation under $k$-refinement by isogeometric discretization methods. They also form a theoretical foundation for the outperformance of smooth spline discretizations of eigenvalue problems that has been numerically observed in the literature, and for optimality of geometric multigrid solvers in the isogeometric analysis context., Comment: 31 pages, 2 figures. Fixed a typo. Article published in M3AS

Multi-degree splines are smooth piecewise-polynomial functions where the pieces can have different degrees. We describe a simple algorithmic construction of a set of basis functions for the space of multi-degree splines, with similar properties to standard B-splines. These basis functions are called multi-degree B-splines (or MDB-splines). The construction relies on an extraction operator that represents all MDB-splines as linear combinations of local B-splines of different degrees. This enables the use of existing efficient algorithms for B-spline evaluations and refinements in the context of multi-degree splines. A Matlab implementation is provided to illustrate the computation and use of MDB-splines.

Alfv\'en-like operators are of interest in magnetohydrodynamics, which is used in plasma physics to study the macroscopic behavior of plasma. Motivated by this important and complex application, we focus on a parameter-dependent curl-div problem that can be seen as a prototype of an Alfv\'en-like operator, and we discretize it using isogeometric analysis based on tensor-product B-splines. The involved coefficient matrices can be very ill-conditioned, so that standard numerical solution methods perform quite poorly here. In order to overcome the difficulties caused by such ill-conditioning, a two-step strategy is proposed. First, we conduct a detailed spectral study of the coefficient matrices, highlighting the critical dependence on the different physical and approximation parameters. Second, we exploit such spectral information to design fast iterative solvers for the corresponding linear systems. For the first goal we apply the theory of (multilevel block) Toeplitz and generalized locally Toeplitz sequences, while for the second we use a combination of multigrid techniques and preconditioned Krylov solvers. Several numerical tests are provided both for the study of the spectral problem and for the solution of the corresponding linear systems., Comment: 26 pages, 20 figures

Isogeometric analysis is a recently developed framework based on finite element analysis, where the simple building blocks in geometry and solution space are replaced by more complex and geometrically-oriented compounds. Box splines are an established tool to model complex geometry, and form an intermediate approach between classical tensor-product B-splines and splines over triangulations. Local refinement can be achieved by considering hierarchically nested sequences of box spline spaces. Since box splines do not offer special elements to impose boundary conditions for the numerical solution of partial differential equations (PDEs), we discuss a weak treatment of such boundary conditions. Along the domain boundary, an appropriate domain strip is introduced to enforce the boundary conditions in a weak sense. The thickness of the strip is adaptively defined in order to avoid unnecessary computations. Numerical examples show the optimal convergence rate of box splines and their hierarchical variants for the solution of PDEs.