Manifold-based surface construction techniques from geometric modelling provide an elegant and flexible framework for generating basis functions on surfaces with arbitrary topology grimm1995modeling ; navau2000modeling ; ying2004simple ; della2008construction . They combine manifold descriptions from differential geometry, see e.g. schutz1980geometrical ; doCarmo1976differential , with the flexibility of the partition of unity framework from numerical analysis Melenk1996a ; babuvska1997partition . If a manifold surface in can be mapped onto a single planar parametric domain in , it is straightforward to obtain partition of unity basis functions of any desired regularity on the parametric domain. Although it is impossible to map a surface with arbitrary topology onto a single parametric domain, it can always be represented as an atlas composed of a number of charts. The charts consist of a planar domains in that map onto the manifold surface in . The planar chart domains in are not connected and transition functions are used to navigate between the different domains. The manifold-based basis functions are obtained by simply applying the partition of unity method on the collection of chart domains majeedCirak:2016 ; zhang:2019 . The flexibility of the original partition of unity method carries over to the manifold case. The partition of unity functions, referred to as blending functions, in this paper, the local approximants on each chart domain and the transition functions can all be chosen to fit the requirements of the application at hand.
The definition of smooth functions over unstructured meshes, such as shown in Figure 1, or multi-patch geometries has always been of vital interest in the context of isogeometric analysis, cf. Hughes:2005aa ; cottrell2009isogeometric . Approaches for the definition of such functions include the subdivision surfaces Cirak:2000aa ; Cirak:2011aa ; Peters:2008aa ; zhang2018subdivision , constructions that are almost everywhere sangalli2016unstructured ; buchegger2016adaptively , constructions scott2013isogeometric ; scott2014isogeometric ; nguyen2016c ; collin2016analysis ; kapl2017isogeometric ; kapl2017analysis ; kapl2018construction and constructions with singular parameterisations toshniwal2017smooth ; toshniwal2017multi . Since most conventional CAD models are based on B-spline or NURBS surfaces, any isogeometric analysis suitable construction should be able to reproduce tensor-product B-splines and NURBS. For this reason we especially focus on manifold-based surfaces that can reproduce tensor-product B-splines in regular portions of the mesh.
In this paper, we exploit the flexibility of the partition of unity method to devise manifold-based basis functions that can reproduce or are identical to B-splines. The proposed techniques are introduced for the sake of clarity with the help of univariate B-splines. Evidently, manifold-based basis functions can reproduce B-splines only on structured regions of an unstructured mesh with extraordinary vertices, i.e. non-boundary vertices with different than four attached elements. In the vicinity of extraordinary vertices the basis functions consist of a local polynomial approximant that can smoothly blend with the surrounding B-spline reproducing basis functions. The extent of the transition region depends on the size of the chosen chart domain, which consists of an -ring of elements around each vertex. We consider several different choices for the weight functions that lead to B-spline reproducing basis functions. Especially promising are weight functions which are defined as a linear combination of B-splines defined on a grid obtained by subdividing the elements multiple times. They provide a partition of unity without normalisation and, hence, lead to polynomial manifold-based basis functions. To obtain manifold-basis functions that are identical to B-splines the local polynomial approximants have to be altered. Whereas in the original manifold constructions the local polynomial approximants are , they have to be chosen to have the same smoothness as the considered B-splines.
On structured meshes the approximation properties of manifold-based basis functions can be inferred from the theory presented in Melenk and Babuska Melenk1996a
. The summation of the local errors on the charts gives a global error estimate under some smoothness assumptions on the weight functions. The local error, for instance inor norms, is bounded by , where is the diameter of the chart domain and is the degree of the polynomials contained in the local approximant. Each chart domain consists of -rings of elements around a vertex so that there are on unstructured meshes multiple types of chart domains depending on the local connectivity of the control mesh. The local connectivity of the control mesh determines the type of transition function used in the manifold construction. When a control mesh is refined by quadrisecting its elements, all the newly introduced vertices are ordinary. That is, for points close to the extraordinary vertices the type of the transition function used depends on the refinement level of the mesh. Therefore, the theory presented in Melenk and Babuska Melenk1996a has to be extended to cover the extraordinary vertices, which we do not attempt in this paper.
The outline of this paper is as follows. In Section 2 we review the manifold basis functions as introduced in majeedCirak:2016 . Although only univariate basis functions on polygonal control meshes are considered, no specific choices for the weight functions, local approximants and the transition functions are given so that the presented theory is applicable to the multivariate case as well. Subsequently, in Section 3 several specific choices first for weight functions and then for local approximants are introduced. More specifically, in Section 3.1 five choices for the weight functions are proposed, two of which yield manifold basis functions that can reproduce B-splines. In Section 3.2 it is illustrated how to choose the local approximants so that manifold-basis functions are identical to B-splines. Finally, in Section 4 we provide a summary and comparison of the different proposed constructions.
2 Review of manifold-based basis functions
We provide a brief informal review of univariate manifold basis functions for curves with the aim to fix ideas and notation.
2.1 Basic approach
Given is a control polygon with the vertex coordinates which describes the manifold curve . To begin with, the control polygon and the curve are assumed to be closed to sidestep the discussion of boundaries. The curve is composed of a set of overlapping subdomains
Each subdomain is associated with a vertex of the control polygon in a manner yet to be described. The subdomains are obtained from corresponding planar domains with a mapping
The pair consisting of is called a chart. In the following we refer to as the chart domain or simply as the chart. If it is clear from context, we neglect the mapping . Each chart domain has its own coordinate system with the coordinates . The coordinates of points on the intersection between two subdomains and can be transformed with transition functions, that is,
Here is the pull-back of the intersection of the two subdomains.
For constructing a smooth approximant, on each chart domain we have given a blending (or, weight) function with
which have to satisfy
and have to be smooth. In addition, at the chart domain boundaries, and its derivatives up to -th order have to be zero. On each chart domain also a local approximant is defined. The approximant is usually expressed in a polynomial basis, like the power, Lagrangian or the Bézier basis. In this paper both the Lagrangian and Bézier basis of a fixed degree are used. However, the basis and the degree of the approximant may be different on every chart. Hence, having the local bases of degree on chart gives the local approximant
and, in turn, the global approximant
as well as the global basis
Note that the index from the basis may be dropped when on each chart domain the same basis is used, as in the present paper.
Next, each chart domain and its image are associated with segments/elements in the -neighbourhood of the control polygon around the vertex . That is, there are as many charts as vertices in the mesh. An -neighbourhood of a vertex is defined as the union of elements that contain the vertex. The -neighbourhood is defined recursively as the union of all 1-neighbourhoods of the -neighbourhood vertices. The number of segments associated with a chart is hence . In turn, each segment is present in charts. See Figure 2 for a construction with each chart domain consisting of the two segments in the -neighbourhood of the vertex .
In Section 3 we consider the span of in (9) as the analysis space on . That is, the basis together with the corresponding set of coefficients is used for -fitting. However, the set of coefficients lack an intuitive interpretation, similar to the control vertices of splines, so that the basis (9) is not suitable for geometric modelling. Hence, in the following we define a design space
as a suitable subspace having degrees of freedom corresponding to the vertex coordinates.
2.2 Mesh-based approach
On each chart domain the coefficients of the local approximant can be assigned to vertices in the -neighbourhood, see Figure 2. Each vertex is present on charts which leads to a coupling between the coefficients of the local approximants (7) of the involved charts. If the number of the coefficients of the local approximant is less than the number of vertices in the chart a least squares fitting has to be applied
where denotes the least-squares projection matrix, →f is an array of the scalar vertex coefficients for the entire polygon (one scalar per vertex) and a gather matrix filled with ones and zeros to pick up the control vertex coefficients for a particular chart from . Hence, the global approximant (8) can be rewritten as
Finally, the manifold curve can be obtained by replacing the array of vertex scalar coefficients with the array of given vertex coordinates , so that each map (2) reads
To summarise so far, each segment on the control mesh has a unique set of corresponding segments on several planar chart domains . The introduced manifold construction ensures that the images of the set of segments from disparate planar charts are identical on the manifold . To advance a more classical finite element interpretation, each segment on the manifold with the index represents an element and has a corresponding reference element to evaluate the element integrals111The index for the reference element is usually dropped because all of them can be assumed to have the same domain.. The mapping of the parent element onto the manifold is composed of two maps
with and being a segment on the chart domain or manifold, respectively. This implies for the field variables in a reference element ,
so that the required smoothness of depends on the collection of maps and on the respective chart domains.
The approximation of the field variables with (11) leads for the element to the following definition of finite element basis functions :
In the following we denote the basis functions with and the mapping from the reference element to the chart with . This notation is not precise when charts have different geometries and number of vertices.
It is clear that the smoothness of the basis functions depends on the smoothness of blending functions , local basis as well as the mappings . For instance, in the two-dimensional continuous construction introduced in majeedCirak:2016 , the blending functions are chosen to be (normalised) B-splines of degree , the local basis are chosen to be a polynomial basis and are conformal maps. Furthermore, each chart domain consists of the elements in the 1-neigborhood of the corresponding vertex. Figure 2 illustrates this construction for continuous basis functions in the univariate case.
3 Reproduction of B-splines
We consider again the global basis (9) for analysis, repeated here for convenience,
and discuss how to choose the blending functions and the local basis so that B-splines are a subset of the basis . To reproduce B-splines it is sufficient to choose either the blending functions or the local basis suitably. In the following we introduce several choices for the blending functions and local basis and comment on their extendability to the bivariate case. Evidently, B-splines are defined on a structured mesh so that manifold-based basis functions will only reproduce B-splines on the parts of the mesh with no extraordinary vertices.
In the univariate case, the parameter domain of the manifold curve can be assumed to be one single finite interval . Due to the choice of the single finite interval the transition functions are identity maps and are affine maps, both are omitted in the following. Without loss of generality, the parameter domain is uniformly partitioned with inner nodes with the coordinates . Moreover, we use the notation for polynomials of degree , for B-splines of degree and continuity and define the space
as a function space on .
3.1 Blending function choices
We fix the local basis to be polynomials of some prescribed degree and study how to choose the blending functions to obtain B-spline reproducing manifold-based basis functions. As the B-splines form a local, non-negative partition of unity, they are used as blending functions. Specifically, the blending functions will be chosen either as
standard B-splines of maximum smoothness,
rational B-spline functions, or
linear combinations of B-splines.
We compare the different approaches in terms of maximum number of overlapping charts at any point of the domain, expected approximation order as well as smoothness properties.
Piecewise linear continuous blending functions
In case of linear B-spline blending functions of degree each chart domain contains three knots. This means that any point on the parameter domain is present on two different chart domains. Having hat functions as blending functions and polynomials of degree as local functions on every chart, we reproduce continuous, piecewise polynomials of degree , i.e., B-splines. Hence, and . In Figure 3 the manifold-based basis obtained with linear B-spline blending functions and a cubic Bézier local basis is shown. The corresponding hat function and the cubic Bézier basis are depicted in Figure 4.
Extending the construction to surfaces, we obtain functions that are tensor-product polynomials within every quadrilateral element and continuous across every edge. Hence, we can only reproduce continuous basis functions on unstructured quadrilateral Bézier meshes.
Higher order continuous B-spline blending functions
The generalisation of the linear B-spline blending functions to the higher order B-splines with is straightforward. As shown in Figure 5 the support of each basis function on the parameter domain is defined as a chart. Therefore, for B-splines of degree and smoothness in any point of the domain charts overlap.
Taking splines of degree as blending functions, together with polynomials of degree as local basis, we reproduce B-splines of degree and smoothness . Hence, we have . Let be the subspace of without global polynomials, i.e.,
then we have
The dimension of is independent of and is equal to the number of inner knots of the spline space , i.e., , and
As the span of contains splines of degree , we can expect an approximation order of in . Note that since the space of blending functions contains polynomials, the functions in are linearly dependent for , as we then have .
The treatment of the boundary is not straightforward. In Figure 5 two different options for choosing cubic blending functions with is presented. The treatment of boundaries becomes relevant when extending the construction to surfaces. This leads to surfaces that are if the mesh is regular. The construction on the top contains all splines, but leads to smooth surfaces at the extraordinary vertices. The construction on the bottom reproduces only a subspace of all B-splines, but generates surfaces that are smooth everywhere. However, as the overlap between the charts is large, the construction becomes cumbersome. Especially, in the bivariate case there can be several extraordinary vertices within one chart domain which can render their smooth parametrisation challenging.
For this reason, in the following, we consider blending functions that have a small support, but generate a smooth basis.
Rational B-spline blending functions
To circumvent the difficulties that arise from using chart domains with a large number of overlaps, we construct blending functions that lead to only two overlapping chart domains in any point of the domain. That is, we consider blending functions with a overlaps similar to the hat functions in Figure 4 (left), but possess the smoothness of higher order B-spline blending functions as in Figure 5 (bottom). To reproduce this behaviour, we first select a linear combination of the B-splines in Figure 5 as blending functions, so that that the supports of no more than two blending functions overlap at the same time. In addition the B-splines to be used as blending functions are chosen from a suitably scaled coordinate system. See Figure 6 for a construction with cubic B-splines with defined on a coordinate axis scaled by a factor . Here, the blending functions are defined as
where are the B-spline basis functions. Note that every other B-spline is chosen as a blending function and their sum does not add up to one. To obtain a partition of unity, all functions are divided by their sum. The resulting blending functions are then piecewise rational, as depicted in Figure 6. The resulting manifold-based basis functions are also rational and their numerical integration may need more quadrature points than polynomial basis functions of similar order and smoothness.
Linear combinations of B-splines as blending functions
One can take linear combinations of consecutive B-splines as blending functions to obtain polynomial blending functions. We consider only cubic B-splines with and express the coefficients in the linear combination as masks . To begin with, the B-splines to be used as blending functions are defined on a coordinate axis scaled by a factor and the blending function mask is , see Figure 7. That is, the blending functions are obtained as
Alternatively, it is possible to use B-splines defined on a coordinate axis scaled by a factor with a mask so that
Comparison of different blending function choices
To summarise, the constructions using rational functions or sums of B-splines generate charts that have only a one ring overlap. The same is true for piecewise linear hat functions as blending functions. Moreover, the dimension of the global function space is the same in these three cases. However, for piecewise linears, the functions are only , whereas for the other three approaches the smoothness is . In the following we compare the approximation power of the respective approaches.
In Figure 10 we show log-error plots when performing -fitting onto a given function. In this example we considered a sine function over the unit interval. Here, the mesh size satisfies , where we used levels . We compare linear blending functions with local polynomials of degree (blue line) and (red line), as in Subsection 3.1. The former has a theoretical convergence rate of in , while the latter has a theoretical rate of . Both discretizations are only. We moreover compare rational blending functions as in (17) (yellow line), linear combinations of splines as blending functions as in (18) (green line) or in (19) (purple line) of degree and local polynomials of degree . In all three cases, the expected convergence rate is .
We compare all constructions with uniform cubic B-splines of mesh size (orange line). They can be interpreted as a manifold construction as in Subsection 3.1 with and . This construction yields the highest error. Similarly high errors are observed for the basis from Subsection 3.1 with (resulting in piecewise cubics) and the lowest error for the same construction with (resulting in piecewise quartics). The bases constructed in Subsection 3.1 produce rates depending on the polynomial degree of the resulting splines. As all splines of a given degree and varying smoothness converge similarly (here e.g. , ), we have omitted this case in Figure 10. Note that the weight functions with boundary correction in Figure 5 (bottom) will not converge optimally without increasing the degree of local functions close to the boundary.
When comparing the constructions from Subsections 3.1 and 3.1, it turns out that all three converge with optimal rates of order and with significantly smaller constant when compared to piecewise polynomials of degree . This means that, when fitting onto a smooth, univariate function, manifold constructions yield a better approximation than standard B-splines, even though the manifold constructions do not reproduce B-splines. Among the manifold constructions, the non-rational variant (18) from Subsection 3.1 seems to be the faster.
It is reasonable to assume that a B-spline compatible manifold construction further improves the approximation properties. Therefore we introduce constructions with modified local functions that reproduce B-splines.
3.2 Local approximants
We discuss next how to choose the local basis so that the manifold-based construction reproduces B-splines of maximum smoothness. Here, the blending functions have only to satisfy the partition of unity property. Hence, any one of the blending functions introduced in Section 3.1 can be used. To avoid the complications arising from chart domains with large number of overlaps, we consider only blending functions which lead to two overlapping charts. In addition, for the sake of concreteness we focus in the following on cubic B-splines and note that the proposed construction carries over to arbitrary degree.
The global manifold-based approximant (8) on a parameter domain consisting of a single finite interval is given by
It is required that this approximant is equal to a B-spline over all or some of the chart domains . The cubic B-spline approximant is defined as
where are the control point coefficients. This approximant can be expressed as a weighted sum of chart domain contributions by multiplication with the partition of unity function, that is,
The support of each of the terms is strictly restricted to one chart domain, see Figure 11. Note that the local basis in Figure 11 (bottom) consists of five functions, which are scaled differently due to the multiplication with the blending function. Term by term matching of the manifold-based (20) and the weighted B-spline approximants (22) requires on every chart domain
This equation yields a set of equations for determining the coefficients in dependence of the coefficients . The cubic B-spline has one knot with , , at the centre of the chart domain so that has to consist out of two pieces. Indeed, it is sufficient to consider in each half and of the chart domain a separate polynomial approximant. Choosing in each half a Bézier basis the coefficients can simply be obtained by Bézier extraction as a linear combination of the B-spline coefficients . Or more generally, the coefficients are obtained by solving a small linear system of equations obtained by collocating (23) at four distinct points (for a cubic B-spline) within the segment. Since is not known (23) has to be considered for each of the four non-zero B-spline basis functions individually. Similar to (10), this gives a relation between the two sets of coefficients expressed as
Note that the projection matrix depends on the specific local basis chosen and is here the same on all the chart domains. Introducing the obtained coefficients into (20) yields the manifold-based basis functions, which are by design the same as the B-spline basis functions.
In the bivariate case the B-spline approximant (21) is only available on parts of the mesh with a tensor-product structure. In the vicinity of extraordinary vertices there is no representation as in (23). In such regions, as in the original manifold construction introduced in Section 2.2 a continuous polynomial approximant has to be fitted to the control mesh coefficients. The manifold construction ensures that the global approximant has the desired smoothness properties.
We developed new manifold-based b-spline basis functions by using the manifold-based surface construction techniques from geometric modelling. As illustrated the manifold-based surface construction techniques can be understood as the extension of the partition of unity method to manifolds. Specific choices for the blending functions and local approximants yield B-splines on structured control meshes. Due to the flexibility of the partition of unity method several such choices are possible. We introduced in total five different choices for the blending functions two of which reproduce B-splines. In addition, we introduced one choice for the local approximant that leads to B-splines.
In Table 1 the properties of the manifold-based basis functions obtained from each of the six different choices are listed. For finite elements polynomial basis functions are to be preferred because they usually require less quadrature points to integrate. The number of breaking points within a finite element gives out of how many smoothly attached pieces a basis function consists. For efficient numerical integration the breaking points of the basis function have to be considered so that constructions with fewer breaking points are to be preferred. In the multivariate case, on unstructured meshes only constructions which require only one-ring of elements around each vertex as a chart domain are viable. If the chart domain consists out of more than one ring of elements, there can be several extraordinary vertices in a chart which makes their parametrisation challenging. In the regular setting the approximation order of the introduced constructions can be inferred from Melenk and Babuska Melenk1996a . In Table 1 the higher order convergence of the first two constructions is remarkable. The first construction yields however only basis functions and the second construction requires charts with several rings of elements. Overall, the most promising constructions for finite elements appear to be the blending functions assembled from B-splines introduced in Section 3.1 and the local approximant introduced in Section 3.2. We note that the smoothness of the two resulting basis functions is with . In closing, we note that the mathematical and numerical study of the introduced constructions on unstructured meshes provides a promising area for future research.
|rational /||B-spline||chart size||breaking||approx.|
|in Sect. 3.1||p||✓||one-ring||0|
|in Sect. 3.1||p||✓||-ring||0|
|in Sect. 3.1||r||one-ring||1|
|in Sect. 3.1||p||one-ring||2|
|in Sect. 3.1||p||one-ring||3|
|as in Sect. 3.2||p||✓||one-ring|
The authors would like to acknowledge the kind hospitality of the Erwin Schrödinger International Institute for Mathematics and Physics (ESI), where some of this research was carried out as part of the thematic programme Numerical Analysis of Complex PDE Models in the Sciences.
- (1) Babuška, I., Melenk, J.M.: The partition of unity method. International journal for numerical methods in engineering 40, 727–758 (1997)
- (2) Buchegger, F., Jüttler, B., Mantzaflaris, A.: Adaptively refined multi-patch b-splines with enhanced smoothness. Applied Mathematics and Computation 272, 159–172 (2016)
- (3) do Carmo, M.P.: Differential geometry of curves and surfaces. Prentice-Hall, Englewood Cliffs, NJ (1976)
- (4) Cirak, F., Long, Q.: Subdivision shells with exact boundary control and non-manifold geometry. International Journal for Numerical Methods in Engineering 88, 897–923 (2011)
- (5) Cirak, F., Ortiz, M., Schröder, P.: Subdivision surfaces: A new paradigm for thin-shell finite-element analysis. International Journal for Numerical Methods in Engineering 47, 2039–2072 (2000)
- (6) Collin, A., Sangalli, G., Takacs, T.: Analysis-suitable multi-patch parametrizations for isogeometric spaces. Computer Aided Geometric Design 47, 93–113 (2016)
- (7) Cottrell, J.A., Hughes, T.J., Bazilevs, Y.: Isogeometric analysis: toward integration of CAD and FEA. John Wiley & Sons (2009)
- (8) Della Vecchia, G., Jüttler, B., Kim, M.S.: A construction of rational manifold surfaces of arbitrary topology and smoothness from triangular meshes. Computer Aided Geometric Design 25, 801–815 (2008)
- (9) Grimm, C.M., Hughes, J.F.: Modeling surfaces of arbitrary topology using manifolds. In: SIGGRAPH 1995 Conference Proceedings, pp. 359–368 (1995)
- (10) Hughes, T.J.R., Cottrell, J.A., Bazilevs, Y.: Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering 194, 4135–4195 (2005)
- (11) Kapl, M., Buchegger, F., Bercovier, M., Jüttler, B.: Isogeometric analysis with geometrically continuous functions on planar multi-patch geometries. Computer Methods in Applied Mechanics and Engineering 316, 209–234 (2017)
- (12) Kapl, M., Sangalli, G., Takacs, T.: Analysis-suitable multi-patch isogeometric spaces: basis and dual basis. arXiv preprint arXiv:1711.05161 (2017)
- (13) Kapl, M., Sangalli, G., Takacs, T.: Construction of analysis-suitable planar multi-patch parameterizations. Computer-Aided Design 97, 41–55 (2018)
- (14) Majeed, M., Cirak, F.: Isogeometric analysis using manifold-based smooth basis functions. Computer Methods in Applied Mechanics and Engineering 316, 547–567 (2017)
- (15) Melenk, J.M., Babuska, I.: The partition of unity finite element method: Basic theory and applications. Computer Methods in Applied Mechanics and Engineering 139, 289–314 (1996)
- (16) Navau, J.C., Garcia, N.P.: Modeling surfaces from meshes of arbitrary topology. Computer Aided Geometric Design 17, 643–671 (2000)
- (17) Nguyen, T., Karčiauskas, K., Peters, J.: finite elements on non-tensor-product 2d and 3d manifolds. Applied Mathematics and Computation 272, 148–158 (2016)
- (18) Peters, J., Reif, U.: Subdivision Surfaces. Springer Series in Geometry and Computing. Springer (2008)
- (19) Sangalli, G., Takacs, T., Vázquez, R.: Unstructured spline spaces for isogeometric analysis based on spline manifolds. Computer Aided Geometric Design 47, 61–82 (2016)
- (20) Schutz, B.F.: Geometrical methods of mathematical physics. Cambridge University Press, Cambridge, UK (1980)
- (21) Scott, M.A., Simpson, R.N., Evans, J.A., Lipton, S., Bordas, S.P., Hughes, T.J., Sederberg, T.W.: Isogeometric boundary element analysis using unstructured T-splines. Computer Methods in Applied Mechanics and Engineering 254, 197–221 (2013)
- (22) Scott, M.A., Thomas, D.C., Evans, E.J.: Isogeometric spline forests. Computer Methods in Applied Mechanics and Engineering 269, 222–264 (2014)
- (23) Toshniwal, D., Speleers, H., Hiemstra, R.R., Hughes, T.J.: Multi-degree smooth polar splines: A framework for geometric modeling and isogeometric analysis. Computer Methods in Applied Mechanics and Engineering 316, 1005–1061 (2017)
- (24) Toshniwal, D., Speleers, H., Hughes, T.J.: Smooth cubic spline spaces on unstructured quadrilateral meshes with particular emphasis on extraordinary points: Geometric design and isogeometric analysis considerations. Computer Methods in Applied Mechanics and Engineering 327, 411–458 (2017)
- (25) Ying, L., Zorin, D.: A simple manifold-based construction of surfaces of arbitrary smoothness. In: SIGGRAPH 2004 Conference Proceedings, pp. 271–275 (2004)
- (26) Zhang, Q., Cirak, F.: Manifold-based isogeometric analysis basis functions with prescribed sharp features. Submitted for publication. (2019)
- (27) Zhang, Q., Sabin, M., Cirak, F.: Subdivision surfaces with isogeometric analysis adapted refinement weights. Computer-Aided Design 102, 104–114 (2018)