Multi-degree splines are piecewise polynomial functions having sections of different degrees. They offer significant advantages over the classical uniform-degree framework, as they allow for modeling complex geometries with fewer degrees of freedom and, at the same time, for a more efficient engineering analysis. Moreover they possess a set of basis functions with similar properties to standard B-splines. In this paper we develop an algorithm for efficient evaluation of multi-degree B-splines, which, unlike previous approaches, is numerically stable. The proposed method consists in explicitly constructing a mapping between a known basis and the multi-degree B-spline basis of the space of interest, exploiting the fact that the two bases are related by a sequence of knot insertion and/or degree elevation steps and performing only numerically stable operations. In addition to theoretically justifying the stability of the algorithm, we will illustrate its performance through numerical experiments that will serve us to demonstrate its excellent behavior in comparison with existing methods, which, in some cases, suffer from apparent numerical problems.
Beccari C.V., Casciola G. (2022). Stable numerical evaluation of multi-degree B-splines. JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 400, 1-21 [10.1016/j.cam.2021.113743].
Stable numerical evaluation of multi-degree B-splines
Beccari C. V.
;Casciola G.
2022
Abstract
Multi-degree splines are piecewise polynomial functions having sections of different degrees. They offer significant advantages over the classical uniform-degree framework, as they allow for modeling complex geometries with fewer degrees of freedom and, at the same time, for a more efficient engineering analysis. Moreover they possess a set of basis functions with similar properties to standard B-splines. In this paper we develop an algorithm for efficient evaluation of multi-degree B-splines, which, unlike previous approaches, is numerically stable. The proposed method consists in explicitly constructing a mapping between a known basis and the multi-degree B-spline basis of the space of interest, exploiting the fact that the two bases are related by a sequence of knot insertion and/or degree elevation steps and performing only numerically stable operations. In addition to theoretically justifying the stability of the algorithm, we will illustrate its performance through numerical experiments that will serve us to demonstrate its excellent behavior in comparison with existing methods, which, in some cases, suffer from apparent numerical problems.File | Dimensione | Formato | |
---|---|---|---|
MD_splines_GrevilleR2.pdf.pdf
accesso aperto
Tipo:
Postprint
Licenza:
Licenza per Accesso Aperto. Creative Commons Attribuzione - Non commerciale - Non opere derivate (CCBYNCND)
Dimensione
389.55 kB
Formato
Adobe PDF
|
389.55 kB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.