## Isogeometric collocation on planar multi-patch domains

Mario Kapl a, ∗ , Vito Vitrih b a Johann Radon Institute for Computational and Applied Mathematics, Austrian Academy of Sciences, Linz, Austria b IAM and FAMNIT, University of Primorska, Koper, Slovenia

## Abstract

We present an isogeometric framework based on collocation to construct a C 2 -smooth approximation of the solution of the Poisson's equation over planar bilinearly parameterized multi-patch domains. The construction of the used globally C 2 -smooth discretization space for the partial differential equation is simple and works uniformly for all possible multipatch configurations. The basis of the C 2 -smooth space can be described as the span of three different types of locally supported functions corresponding to the single patches, edges and vertices of the multi-patch domain. For the selection of the collocation points, which is important for the stability and convergence of the collocation problem, two different choices are numerically investigated. The first approach employs the tensor-product Greville abscissae as collocation points, and shows for the multi-patch case the same convergence behavior as for the one-patch case [2], which is suboptimal in particular for odd spline degree. The second approach generalizes the concept of superconvergent points from the one-patch case (cf. [1, 15, 32]) to the multi-patch case. Again, these points possess better convergence properties than Greville abscissae in case of odd spline degree.

Keywords: isogeometric analysis; collocation; superconvergent points; second order continuity; multi-patch domain; Poisson's equation 2010 MSC: 65N35, 65D17, 68U07

## 1. Introduction

Isogeometric Analysis (IgA) is an approach for numerically solving a partial differential equation (PDE) by using the same spline or NURBS space for describing the geometry and for representing the solution of the considered PDE, cf. [4, 10, 18]. Usually, the developed PDE solvers are based on finding a solution of the variational (weak) form of the PDE, since using the weak form requires less smooth functions compared to directly employing the strong form. However, this technique needs the evaluation of integrals by using some particular quadrature rules, and the accuracy of the solution depends on the quality of the numerical integration. In contrast, solving the strong form of the PDE via collocation eliminates integration but requires spaces of higher regularity.

∗ Corresponding author

Email addresses: mario.kapl@ricam.oeaw.ac.at (Mario Kapl), vito.vitrih@upr.si (Vito Vitrih)

While in case of a one-patch domain, the higher continuity of the functions can be easily guaranteed by using spline functions with the desired smoothness within the patch, the construction of C s -smooth ( s ≥ 1) isogeometric spline spaces over multi-patch domains is challenging, and is the task of current research, see e.g. [6-9, 20-23, 28, 29, 33, 34, 40, 41] for s = 1 and [24-27, 39] for s = 2. The design of C s -smooth multi-patch spline spaces is linked to the concept of geometric continuity of multi-patch surfaces, cf. [17, 35], due to the fact that an isogeometric function is C s -smooth on a multi-patch domain if and only if its associated multi-patch graph surface is G s -smooth , cf. [16, 28].

So far, the problem of isogeometric collocation has been mostly explored on one-patch domains, see e.g. [1, 2, 5, 13-15, 30-32, 36, 38]. For analyzing the convergence behavior of a particular approach, the errors are generally computed with respect to L ∞ , W 1 , ∞ and W 2 , ∞ norm, or equivalently with respect to L 2 , H 1 and H 2 norm, respectively. The study of isogeometric collocation methods has started in [2] by using the Greville and Demko abscissae as collocation points. In comparison to the Galerkin approach, both choices show a suboptimal convergence behavior with respect to L 2 ( L ∞ ) norm, namely of orders O ( h p ) and O ( h p -1 ) under h -refinement for even and odd spline degree p , respectively, and additionally a suboptimal convergence order O ( h p -1 ) with respect to H 1 ( W 1 , ∞ ) norm but just in case of odd spline degree p .

In [1], an isogeometric collocation method has been presented, which is based on the computation and on the use of specific collocation points called superconvergent points . The proposed technique possesses only for even spline degree p in case of the L 2 ( L ∞ ) norm a suboptimal convergence behavior of order O ( h p ), and is optimal for all other cases. This is achieved at the expense of solving an overdetermined linear system due to the fact that the number of superconvergent points is larger than the number of degrees of freedom. The two methods [15] and [32] select different specific subsets of the superconvergent points to achieve that the number of collocation points coincides with the number of degrees of freedom. While the points 1 in [15] are chosen in an alternating way, which reduces for the case of odd spline degree p the convergence order also to O ( h p ) in the L 2 ( L ∞ ) norm, the points in [32] are selected in a clustered way to maintain the convergence behavior of [1].

In [1, 2, 14, 32, 38], the concept of isogeometric collocation has been mostly used for solving the Poisson's equation. Further problems and PDEs of interest in the framework of isogeometric collocation are amongst others linear and nonlinear elasticity, elastostatistics and elastodynamics [1, 3, 13-15, 19, 36, 38], plate/shell problems [15, 30] and beam problems [5, 31].

As already mentioned before, most of the existing isogeometric collocation methods have in common that they are mainly restricted to the case of one-patch domains. Two techniques which also deal with the case of multi-patch domains are [3, 19]. While in [3] standard NURBS functions are employed, the work in [19] is based on the use of PHTsplines [12]. Since for both methods the used spline spaces are in general just C 0 -smooth across the patch interfaces, special techniques for the collocation on the interfaces are

1 In [15], the superconvergent points are also denoted as Cauchy Galerkin points , since collocation at these points produces the Galerkin solution exactly.

## required.

The goal of this paper is to compute a globally C 2 -smooth approximation of the solution of the Poisson's equation on planar bilinearly parameterized multi-patch domains by means of isogeometric collocation. For this purpose, a globally C 2 -smooth discretization space is generated and used. The construction of this C 2 -smooth multi-patch space is based on and extends the work [27]. There, a subspace of the entire C 2 -smooth spline space [24] has been generated, which possesses a simpler structure as the entire C 2 -smooth space but still maintains the full approximation properties. The construction of the locally supported basis functions is simple via explicit formulae or by computing the null space of a small system of linear equations. Whereas the basis construction in [27] works uniformly for all possible multi-patch configurations, the number of basis functions depends on the initial geometry. This work overcomes the latter limitation by constructing a C 2 -smooth spline space whose dimension is independent of the initial geometry by e.g. additionally enforcing that the functions have to be C 4 -smooth at the inner vertices and at the boundary vertices of patch valencies greater or equal to three of the multi-patch domain.

The generated C 2 -smooth spline space can be directly employed to perform isogeometric collocation on the multi-patch domain, and no special treatment for the patch interfaces as in [3, 19] is needed. For assembling the system matrices of the collocation problem two different choices of collocation points are studied. On the one hand, we employ the tensorproduct Greville abscissae for the single patches, and on the other hand we generalize the so-called superconvergent points for the one-patch case (cf. [1, 15, 32]) to the multi-patch case. Like in the one-patch case, the superconvergent points possess a better convergence behavior as the Greville points for odd spline degree, namely of orders O ( h p ), O ( h p ) and O ( h p -1 ) in comparison to O ( h p -1 ), O ( h p -1 ) and O ( h p -1 ) with respect to the L 2 ( L ∞ ), H 1 ( W 1 , ∞ ) and H 2 ( W 2 , ∞ ) norm.

The remainder of the paper is organized as follows. Section 2 presents the construction of a particular C 2 -smooth isogeometric spline space over bilinear multi-patch domains, which will be used as a discretization space for solving the Poisson's equation by means of isogeometric collocation. The multi-patch isogeometric collocation method as well as two different strategies for choosing the collocation points are described in Section 3. On the basis of several numerical examples, the convergence behavior under h -refinement of the collocation method is studied in Section 4, and demonstrates the potential of the method for the use in IgA. Finally, we conclude the paper.

## 2. A C 2 -smooth isogeometric spline space

The construction of a particular C 2 -smooth multi-patch isogeometric spline space will be presented, which will be based on and will extend the work [27]. The C 2 -smooth space will be used in Section 3 to build an isogeometric collocation method for solving the Poisson's equation over bilinear multi-patch domains. We will start with the introduction of the multi-patch setting, which will be used throughout the paper.

## 2.1. The multi-patch setting

Let I Ω be an index set, and let Ω and Ω ( ) i , i ∈ I Ω , be open domains in R 2 , such that Ω = ∪ i ∈I Ω Ω ( ) i . Furthermore, let Ω ( ) i , i ∈ I Ω , be quadrangular patches, which are mutually disjoint, and the closures of any two of them have either an empty intersection, possess exactly one common vertex or share the whole common edge. Additionally, the deletion of any vertex of the multi-patch domain Ω does not split Ω into subdomains, whose union would be unconnected. We will further assume that each patch Ω ( ) i is parameterized by a bilinear, bijective and regular geometry mapping F ( ) i ,

↦

$$F ^ { ( i ) } \, \colon [ 0, 1 ] ^ { 2 } \to \mathbb { R } ^ { 2 }, \quad \xi = ( \xi _ { 1 }, \xi _ { 2 } ) \mapsto F ^ { ( i ) } ( \xi ) = F ^ { ( i ) } ( \xi _ { 1 }, \xi _ { 2 } ), \quad i \in \mathcal { I } _ { \Omega },$$

such that Ω ( ) i = F ( ) i ([0 , 1] 2 ). We will also use the splitting of the multi-patch domain Ω into the single patches Ω ( ) i , i ∈ I Ω , edges Γ ( ) i , i ∈ I Γ and vertices v ( ) i , i ∈ I Ξ , i.e.

$$\overline { \Omega } = \bigcup _ { i \in \mathcal { I } _ { \Omega } } \Omega ^ { ( i ) } \, \dot { \cup } \, \bigcup _ { i \in \mathcal { I } _ { \Gamma } } \Gamma ^ { ( i ) } \, \dot { \cup } \, \bigcup _ { i \in \mathcal { I } _ { \Xi } } v ^ { ( i ) },$$

where I Γ and I Ξ are the index sets of the indices of the edges Γ ( ) i and vertices v ( ) i , respectively. We further have I Γ = I Γ I ˙ ∪I Γ B and I Ξ = I Ξ I ˙ ∪I Ξ B with I Ξ B = I Ξ 1 ˙ ∪I Ξ 2 ˙ ∪I Ξ 3 , where I Γ I and I Γ B contain all indices of interfaces and boundary edges, respectively, I Ξ I and I Ξ B collect all indices of inner and boundary vertices, respectively, and I Ξ 1 , I Ξ 2 and I Ξ 3 are further the index sets of the indices of boundary vertices of patch valencies one, two and greater or equal to three, respectively. In addition, the patch valency of a vertex v ( ) i will be denoted by ν i .

Let S p,r h ([0 , 1]) be the univariate spline space of degree p , regularity r and mesh size h = 1 k +1 on the unit interval [0 , 1], constructed from the uniform open knot vector

$$( t _ { 0 } ^ { p, r }, \dots, t _ { 2 p + k ( p - r ) + 1 } ^ { p, r } ) = ( \underbrace { 0, \dots, 0 } _ { ( p + 1 ) - \text{times} }, \underbrace { \frac { 1 } { k + 1 }, \dots, \frac { 1 } { k + 1 } } _ { ( p - r ) - \text{times} }, \dots, \underbrace { \frac { k } { k + 1 }, \dots, \frac { k } { k + 1 } } _ { ( p - r ) - \text{times} }, \underbrace { 1, \dots, 1 } _ { ( p + 1 ) - \text{times} } ),$$

where k is the number of different inner knots. Furthermore, let S p r , h ([0 , 1] 2 ) be the tensorproduct spline space S p,r h ([0 , 1]) ⊗S p,r h ([0 , 1]) on the unit-square [0 , 1] 2 . We denote the Bsplines of the spaces S p,r h ([0 , 1]) and S p r , h ([0 , 1] 2 ) by N p,r j and N p r , j 1 ,j 2 = N p,r j 1 N p,r j 2 , respectively, with j, j 1 , j 2 = 0 1 , , . . . , n -1, and n = p +1+ ( k p -r ). We assume that p ≥ 5, 2 ≤ r ≤ p -3, and that the number of inner knots satisfies k ≥ 9 -p p - -r 2 . These assumptions are necessary to ensure that the constructed C 2 -smooth spline spaces in Section 2.3 will be h -refineable and well-defined, see also [26, 27]. Since the geometry mappings F ( ) i , i ∈ I Ω , are bilinearly parameterized, we trivially have that

$$F ^ { ( i ) } \in \mathcal { S } _ { h } ^ { p, r } ( [ 0, 1 ] ^ { 2 } ) \times \mathcal { S } _ { h } ^ { p, r } ( [ 0, 1 ] ^ { 2 } ).$$

Moreover, we assume for two neighboring patches Ω ( i 0 ) and Ω ( i 1 ) with the common interface Γ ( ) i ⊂ Ω ( i 0 ) ∩ Ω ( i 1 ) that the associated geometry mappings F ( i 0 ) and F ( i 1 ) are parameterized as given in Fig. 1 (left), and assume for the neighboring patches enclosing

an inner or boundary vertex v ( ) i of patch valency ν i ≥ 3 that the associated geometry mappings are parameterized as shown in Fig. 1 (right). For the latter case, we additionally assume that the patches enclosing the vertex v ( ) i are relabeled by Ω ( i 0 ) , . . . , Ω ( i ν i -1 ) , and that the common interfaces 2 of the two-patch subdomains Ω ( i glyph[lscript] ) ∪ Ω ( i glyph[lscript] +1 ) , glyph[lscript] = 0 1 , , . . . , ν i -1, are relabeled by Γ ( i glyph[lscript] +1 ) . In case of an inner vertex v ( ) i we further consider the lower index glyph[lscript] of the indices i glyph[lscript] modulo ν i , which just means that the lower index glyph[lscript] is replaced by the remainder of the division of glyph[lscript] by ν i . Note that the local reparameterizations as shown in Fig. 1 are always possible (if necessary), and hence can be assumed to be satisfied without loss of generality.

I

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-ac569430-ec13-4cb4-9db4-553fb71db6e0.jpg)

i

M

Figure 1: Left: Considering two neighboring patches Ω ( i 0 ) and Ω ( i 1 ) , we can always assume that the two corresponding geometry mappings F ( i 0 ) and F ( i 1 ) are parameterized as shown. Right: For the patches Ω ( i glyph[lscript] ) , glyph[lscript] = 0 1 , , . . . , ν i -1, which enclose the vertex v ( ) i , the corresponding geometry mappings F ( i glyph[lscript] ) can be always reparameterized to be given as visualized.

## 2.2. C 2 -smoothness across patch interfaces

The space of C 2 -smooth isogeometric functions on Ω is given as

$$\mathcal { V } = \left \{ \phi \in C ^ { 2 } ( \overline { \Omega } ) \, \colon \, \phi | _ { \overline { \Omega } ^ { ( i ) } } \in \mathcal { S } _ { h } ^ { p, r } ( [ 0, 1 ] ^ { 2 } ) \circ ( F ^ { ( i ) } ) ^ { - 1 }, \ i \in \mathcal { I } _ { \Omega } \right \}.$$

This space can be characterized by means of the concept of geometric continuity of multipatch surfaces, cf. [17, 35]. An isogeometric function φ belongs to the space V if and only if for any two neighboring patches Ω ( i 0 ) and Ω ( i 1 ) with the common interface Γ ( ) i = Ω ( i 0 ) ∩ Ω ( i 1 ) , i ∈ I Γ I , the associated graph surface patches [ F ( i 0 ) φ ◦ F ( i 0 ) ] and [ F ( i 1 ) φ ◦ F ( i 1 ) ]

2 In case of a boundary vertex v ( ) i , we denote the boundary edges by Γ ( i 0 ) and Γ ( i ν i ) .

are G 2 -smooth across their common interface , see e.g. [16, 28]. An equivalent condition to the G 2 -smoothness of the two neighboring graph surface patches is that the two associated spline functions g ( i 0 ) = φ ◦ F ( i 0 ) and g ( i 1 ) = φ ◦ F ( i 1 ) satisfy

$$\underbrace { g ^ { ( i _ { 0 } ) } ( 0, \xi ) = g ^ { ( i _ { 1 } ) } ( 0, \xi ) } _ { \coloneqq g _ { 0 } ( \xi ) },$$

$$\underbrace { \underbrace { D _ { \xi _ { 2 } } g ^ { ( i _ { 0 } ) } ( 0, \xi ) - \beta ^ { ( i _ { 0 } ) } _ { \Gamma ^ { ( i ) } ( \xi ) g ^ { \prime } _ { 0 } ( \xi ) } } { \alpha ^ { ( i _ { 0 } ) } _ { \Gamma ^ { ( i ) } ( \xi ) } } = \frac { D _ { \xi _ { 2 } } g ^ { ( i _ { 1 } ) } ( 0, \xi ) - \beta ^ { ( i _ { 1 } ) } _ { \Gamma ^ { ( i ) } ( \xi ) g ^ { \prime } _ { 0 } ( \xi ) } } { \alpha ^ { ( i _ { 1 } ) } _ { \Gamma ^ { ( i ) } ( \xi ) } }, } _ { \coloneqq g _ { 1 } ( \xi ) }, \\.$$

and

$$\frac { D _ { \xi _ { 1 } \xi _ { 1 } } g ^ { ( i _ { 0 } ) } ( 0, \xi ) - ( \beta _ { \Gamma ^ { ( i ) } ( i ) } ^ { ( i _ { 0 } ) } ( \xi ) ) ^ { 2 } g _ { 0 } ^ { \prime \prime } ( \xi ) - 2 \alpha _ { \Gamma ^ { ( i ) } ( i ) } \beta _ { \Gamma ^ { ( i ) } ( i ) } ^ { ( i _ { 0 } ) } ( \xi ) g _ { 1 } ^ { \prime } ( \xi ) } { ( \alpha ^ { ( i _ { 0 } ) } ( \xi ) ) ^ { 2 } } } = \\ \frac { D _ { \xi _ { 1 } \xi _ { 1 } } g ^ { ( i _ { 1 } ) } ( 0, \xi ) - ( \beta _ { \Gamma ^ { ( i _ { 1 } ) } ( i ) } ^ { ( i _ { 1 } ) } ( \xi ) ) ^ { 2 } g _ { 0 } ^ { \prime \prime } ( \xi ) - 2 \alpha _ { \Gamma ^ { ( i ) } ( i ) } ^ { ( i _ { 1 } ) } \beta _ { \Gamma ^ { ( i _ { 1 } ) } ( i ) } ( \xi ) g _ { 1 } ^ { \prime } ( \xi ) } { ( \alpha ^ { ( i _ { 1 } ) } ( \xi ) ) ^ { 2 } },$$

for all ξ ∈ [0 , 1], cf. [26, Lemma 1], where α ( τ ) Γ ( ) i and β ( τ ) Γ ( ) i , τ ∈ { i 0 , i 1 } , are linear polynomials given by

$$\alpha _ { \Gamma ^ { ( i ) } } ^ { ( \tau ) } ( \xi ) = c _ { 1 } \det [ D _ { \xi _ { 1 } } \mathcal { F } ^ { ( \tau ) } ( 0, \xi ), D _ { \xi _ { 2 } } \mathcal { F } ^ { ( \tau ) } ( 0, \xi ) ],$$

and

$$\beta _ { \Gamma ^ { ( i ) } } ^ { ( \tau ) } ( \xi ) = \frac { D _ { \xi _ { 1 } } \mathcal { F } ^ { ( \tau ) } ( 0, \xi ) \cdot D _ { \xi _ { 2 } } \mathcal { F } ^ { ( \tau ) } ( 0, \xi ) } { | | D _ { \xi _ { 2 } } \mathcal { F } ^ { ( \tau ) } ( 0, \xi ) | | ^ { 2 } },$$

respectively, with c 1 ∈ R such that

$$| | \alpha _ { \Gamma ^ { ( i ) } } ^ { ( i _ { 0 } ) } + 1 | | _ { L ^ { 2 } } ^ { 2 } + | | \alpha _ { \Gamma ^ { ( i ) } } ^ { ( i _ { 1 } ) } - 1 | | _ { L ^ { 2 } } ^ { 2 }$$

is minimized, cf. [27]. Note that α ( i 0 ) Γ ( ) i &lt; 0 and α ( i 1 ) Γ ( ) i &gt; 0, since the geometry mappings F ( i 0 ) and F ( i 1 ) are regular.

The space V has been studied in different configurations in [24-26]. There, it has been observed that the space V possesses a complex structure and its dimension depends on the initial geometry. In the next subsection, we will consider instead a subspace of the space V , whose structure is simpler and whose dimension is independent of the initial geometry.

## 2.3. A C 2 -smooth isogeometric space

We will describe the construction of a particular C 2 -smooth space W ⊆ V , which will be similar to the one in [27]. There will be two main differences between the newly constructed space W and the space from [27]. Firstly, the space W will not be anymore restricted to containing just C 2 -smooth functions fulfilling the homogeneous boundary conditions of order 2 given by

$$u ( x ) = \frac { \partial u } { \partial n } ( x ) = \triangle u ( x ) = 0, \quad x \in \partial \Omega,$$

which have been used in [27] as boundary conditions for solving the triharmonic equation via its weak form and standard Galerkin discretization. Secondly, the dimension of the space W will now be independent of the initial geometry.

Analogous to [27], the space W is defined as the direct sum

$$\mathcal { W } = \left ( \bigoplus _ { i \in \mathcal { I } _ { \Omega } } \mathcal { W } _ { \Omega ^ { ( i ) } } \right ) \oplus \left ( \bigoplus _ { i \in \mathcal { I } _ { \Gamma } } \mathcal { W } _ { \Gamma ^ { ( i ) } } \right ) \oplus \left ( \bigoplus _ { i \in \mathcal { I } _ { \Xi } } \mathcal { W } _ { v ^ { ( i ) } } \right ),$$

where W Ω ( ) i , W Γ ( ) i and W v ( ) i are subspaces corresponding to the single patches Ω ( ) i , edges Γ ( ) i and vertices v ( ) i , respectively. While the construction of the subspaces W Ω ( ) i and W Γ ( ) i will work analogously to or similar as in [27], and hence will be kept short below, the construction of the subspaces W v ( ) i will differ.

2.3.1. The patch subspace W Ω ( ) i and the edge subspace W Γ ( ) i

As in [27], the patch subspaces W Ω ( ) i , i ∈ I Ω , are given as

$$\mathcal { W } _ { \Omega ^ { ( i ) } } = \text{span} \{ \phi _ { \Omega ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 1 }, j _ { 2 } = 3, 4, \dots, n - 4 \},$$

with the functions

$$\phi _ { \Omega ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \begin{cases} ( N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } \circ ( F ^ { ( i ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \Omega ^ { ( i ) }, \\ 0 & \text{otherwise}. \end{cases}$$

For the edge subspaces W Γ ( ) i , i ∈ I Γ , we have to distinguish between the case of a boundary edge Γ ( ) i , i.e. i ∈ I Γ B , and the case of an interface Γ ( ) i , i ∈ I Γ I . Let us start with the case of a boundary edge Γ ( ) i , i ∈ I Γ B . Assume without loss of generality, that the boundary edge Γ ( ) i is contained in the geometry mapping F ( i 0 ) and is given by F ( i 0 ) ( { } × 0 (0 , 1)). Then the subspace W Γ ( ) i possesses the form

$$\mathcal { W } _ { \Gamma ^ { ( i ) } } = \text{span} \{ \phi _ { \Gamma ^ { ( i ) } _ { B } ; j _ { 1 }, j _ { 2 } } | \ j _ { 2 } = 5 - j _ { 1 }, \dots, n + j _ { 1 } - 6, \ j _ { 1 } = 0, 1, 2 \},$$

where φ Γ ( ) i B ; j 1 ,j 2 is defined similar to (5) as

$$\phi _ { \Gamma _ { B } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \begin{cases} ( N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } \circ ( F ^ { ( i _ { 0 } ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \overline { \Omega ^ { ( i _ { 0 } ) } }, \\ 0 & \text{otherwise}. \end{cases}$$

In case of an interface Γ ( ) i ⊂ Ω ( i 0 ) ∩ Ω ( i 1 ) , i ∈ I Γ I , and assuming without loss of generality that the associated geometry mappings F ( i 0 ) and F ( i 1 ) are parameterized as in Fig. 1 (left), the edge subspace W Γ ( ) i is defined as

$$\mathcal { W } _ { \Gamma ^ { ( i ) } } = \text{span} \{ \phi _ { \Gamma _ { I } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 2 } = 5 - j _ { 1 }, \dots, n _ { j _ { 1 } } + j _ { 1 } - 6, \ j _ { 1 } = 0, 1, 2 \},$$

with the functions

$$\phi _ { \Gamma ^ { ( i ) } _ { I } ; j _ { 1 }, j _ { 2 } } ( x ) = \begin{cases} ( g ^ { ( i _ { 0 } ) } _ { \Gamma ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } \circ ( F ^ { ( i _ { 0 } ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \overline { \Omega ^ { ( i _ { 0 } ) } }, \\ ( g ^ { ( i _ { 1 } ) } _ { \Gamma ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } \circ ( F ^ { ( i _ { 1 } ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \overline { \Omega ^ { ( i _ { 1 } ) } }, \end{cases}$$

where

$$\text{where} \\ g _ { \Gamma ^ { ( \tau ) } ; 0, j _ { 2 } } ^ { ( \tau ) } ( \xi _ { 1 }, \xi _ { 2 } ) & = N _ { j _ { 2 } } ^ { p, r + 2 } ( \xi _ { 2 } ) M _ { 0 } ( \xi _ { 1 } ) + \beta _ { \Gamma ^ { ( \tau ) } } ^ { ( \tau ) } ( \xi _ { 2 } ) ( N _ { j _ { 2 } } ^ { p, r + 2 } ) ^ { \prime } ( \xi _ { 2 } ) M _ { 1 } ( \xi _ { 1 } ) \\ & + \left ( \beta _ { \Gamma ^ { ( \tau ) } ( \xi _ { 2 } ) } ^ { ( \tau ) } ( \xi _ { 2 } ) \right ) ^ { 2 } ( N _ { j _ { 2 } } ^ { p, r + 2 } ) ^ { \prime \prime } ( \xi _ { 2 } ) M _ { 2 } ( \xi _ { 1 } ), \\ g _ { \Gamma ^ { ( \tau ) } ( i ) ; 1, j _ { 2 } } ^ { ( \tau ) } ( \xi _ { 1 }, \xi _ { 2 } ) & = \frac { p } { h } \left ( \alpha _ { \Gamma ^ { ( \tau ) } ( i ) } ( \xi _ { 2 } ) N _ { j _ { 2 } } ^ { p - 1, r + 1 } ( \xi _ { 2 } ) M _ { 1 } ( \xi _ { 1 } ) \\ & + 2 \alpha _ { \Gamma ^ { ( \tau ) } ( i ) } ( \xi _ { 2 } ) \beta _ { \Gamma ^ { ( \tau ) } ( i ) } ( \xi _ { 2 } ) ( N _ { j _ { 2 } } ^ { p - 1, r + 1 } ) ^ { \prime } ( \xi _ { 2 } ) M _ { 2 } ( \xi _ { 1 } ) \right ), \\ g _ { \Gamma ^ { ( \tau ) } ( i ) ; 2, j _ { 2 } } ^ { ( \tau ) } ( \xi _ { 1 }, \xi _ { 2 } ) & = \frac { p ( p - 1 ) } { h ^ { 2 } } \left ( \alpha _ { \Gamma ^ { ( \tau ) } ( i ) } ( \xi _ { 2 } ) \right ) ^ { 2 } N _ { j _ { 2 } } ^ { p - 2, r } ( \xi _ { 2 } ) M _ { 2 } ( \xi _ { 1 } ), \\ \text{for} \, \tau \in \{ i _ { 0 }, i _ { 1 } \}, \text{with}$$

for τ ∈ { i 0 , i 1 } , with

$$M _ { 0 } ( \xi ) = \sum _ { j = 0 } ^ { 2 } N _ { j } ^ { p, r } ( \xi ), \ M _ { 1 } ( \xi ) = \frac { p } { h } \left ( N _ { 1 } ^ { p, r } ( \xi ) + 2 N _ { 2 } ^ { p, r } ( \xi ) \right ), \ M _ { 2 } ( \xi ) = \frac { p ( p - 1 ) } { h ^ { 2 } } N _ { 2 } ^ { p, r } ( \xi ),$$

and

$$n _ { 0 } = \dim \left ( \mathcal { S } _ { h } ^ { p, r + 2 } ( [ 0, 1 ] ) \right ), \ n _ { 1 } = \dim \left ( \mathcal { S } _ { h } ^ { p - 1, r + 1 } ( [ 0, 1 ] ) \right ), \ n _ { 2 } = \dim \left ( \mathcal { S } _ { h } ^ { p - 2, r } ( [ 0, 1 ] ) \right ).$$

All functions of the patch subspaces W Ω ( ) i , i ∈ I Ω , and of the edge subspaces W Γ ( ) i , i ∈ I Γ , are C 2 -smooth on Ω, since the Eqs. (1)-(3) are satisfied for all interfaces Γ ( j ) , j ∈ I Γ I . In case of the patch subspaces W Ω ( ) i and of the edge subspaces W Γ ( ) i for i ∈ I Γ B , the equations are even trivially satisfied, since the functions have vanishing values, gradients and Hessians on all interfaces Γ ( j ) , j ∈ I Γ I . In case of the edge subspaces W Γ ( ) i for i ∈ I Γ I , whereas the equations are also trivially fulfilled for the interfaces Γ ( j ) , j ∈ I Γ I \ { } i , the particular construction of the functions φ Γ ( ) i I ; j 1 ,j 2 ensures that the Eqs. (1)-(3) are satisfied for the non-trivial case of the corresponding interface Γ ( ) i , too, cf. [26, 27].

## 2.3.2. The vertex subspace W v ( ) i

The vertex subspaces W v ( ) i , i ∈ I Ξ , constructed in this work, will differ from the ones in [27]. The new construction of the subspaces will ensure that the dimension of the vertex subspaces will now be independent of the initial geometry. Furthermore, the design of the vertex subspaces W v ( ) i , i ∈ I Ξ B , will not be anymore restricted to any boundary conditions. We will distinguish between four cases, namely if the vertex v ( ) i is an inner vertex, i.e. i ∈ I Ξ I , a boundary vertex of patch valency one, i.e. i ∈ I Ξ 1 , a boundary vertex of patch valency two, i.e. i ∈ I Ξ 2 , or a boundary vertex of patch valency greater or equal to three, i.e. i ∈ I Ξ 3 .

Let us start with the case of an inner vertex v ( ) i , that is i ∈ I Ξ I . Recall that we assume without loss of generality that the geometry mappings F ( i glyph[lscript] ) containing the vertex v ( ) i are relabeled and parameterized as described and shown in Section 2.1 and Fig. 1 (right), respectively. Again, the lower index glyph[lscript] of the indices i glyph[lscript] will be considered modulo ν i . First

of all, we consider analogous to [27], the isogeometric function φ v ( ) i : Ω → R , which is defined as

$$\phi _ { v ^ { ( i ) } } ( x ) = \begin{cases} ( g _ { i _ { \ell } } \circ ( F ^ { ( i _ { \ell } ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \overline { \Omega ^ { ( i _ { \ell } ) } }, \ \ell = 0, 1, \dots, \nu _ { i } - 1, \\ 0 & \text{otherwise}, \end{cases}$$

where the functions g i glyph[lscript] are given as

$$& \text{where the functions } g _ { i _ { \ell } } \text{ are given as} \\ & g _ { i _ { \ell } } ( \xi _ { 1 }, \xi _ { 2 } ) = \underbrace { \sum _ { j = 0 } ^ { 2 } \sum _ { j = 0 } ^ { 4 - j _ { 1 } } a _ { j _ { 1 }, j _ { 2 } } ^ { ( i _ { \ell } ) } g _ { \Gamma ( i _ { \ell } ) ; j _ { 1 }, j _ { 2 } } ( \xi _ { 2 }, \xi _ { 1 } ) } _ { \colon = g _ { i _ { \ell } } ^ { \Gamma ( i _ { \ell } ) } ( \xi _ { 1 }, \xi _ { 2 } ) } } _ { \colon = g _ { i _ { \ell } } ^ { \Gamma ( i _ { \ell } + 1 ) } ( \xi _ { 1 }, \xi _ { 2 } ) } ( \xi _ { 1 }, \xi _ { 2 } ) } \\ & \quad \underbrace { \sum _ { j = 0 } ^ { 2 } \sum _ { j = 0 } ^ { 2 } a _ { j _ { 1 }, j _ { 2 } } ^ { ( i _ { \ell } ) } N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } ( \xi _ { 1 }, \xi _ { 2 } ) } _ { \colon = g _ { i _ { \ell } } ^ { \Gamma ( i _ { \ell } + 1 ) } ( \xi _ { 1 }, \xi _ { 2 } ) } \\ & \quad \underbrace { \sum _ { j = 0 } ^ { 2 } \sum _ { j = 1 } ^ { 2 } a _ { j _ { 1 }, j _ { 2 } } ^ { ( i _ { \ell } ) } N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } ( \xi _ { 1 }, \xi _ { 2 } ) } _ { \colon = g _ { i _ { \ell } } ^ { \Gamma ( i _ { \ell } ) } ( \xi _ { 1 }, \xi _ { 2 } ) } \\ & \text{with coefficients } a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \left }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gam { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamr ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \left } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a } _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i_ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, { \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a ^ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell } } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gam { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { \ell }, \frac { \Gamma ( i _ { \ell } ) } { a _ { i _ { \ell }, \frac { \Gamma ($$

with coefficients a Γ ( i glyph[lscript] ) j 1 ,j 2 , a ( i glyph[lscript] ) j 1 ,j 2 ∈ R , and with the functions g ( i glyph[lscript] ) Γ ( i glyph[lscript] ); j 1 ,j 2 and g ( i glyph[lscript] ) Γ ( i glyph[lscript] +1); 1 j ,j 2 given by (8).

In [27, Lemma 3], it was shown, that the function φ v ( ) i is C 2 -smooth on Ω if the coefficients a Γ ( i glyph[lscript] ) j 1 ,j 2 and a ( i glyph[lscript] ) j 1 ,j 2 are selected in such a way that the following homogeneous linear system

$$\partial _ { \xi _ { 1 } } ^ { i } \partial _ { \xi _ { 2 } } ^ { j } \left ( g _ { i _ { \ell } } ^ { \Gamma ^ { ( i _ { \ell + 1 } ) } } - g _ { i _ { \ell } } ^ { \Gamma ^ { ( i _ { \ell } ) } } \right ) ( 0 ) = 0 \quad \text{and} \quad \partial _ { \xi _ { 1 } } ^ { i } \partial _ { \xi _ { 2 } } ^ { j } \left ( g _ { i _ { \ell } } ^ { \Gamma ^ { ( i _ { \ell + 1 } ) } } - g _ { i _ { \ell } } ^ { \Omega ^ { ( i _ { \ell } ) } } \right ) ( 0 ) = 0, \quad \text{(11)}$$

with 0 ≤ i, j ≤ 2 and glyph[lscript] = 0 1 , , . . . , ν i -1, is satisfied. The construction of a C 2 -smooth isogeometric function φ v ( ) i is based on the idea of generating the spline function g i glyph[lscript] = φ v ( ) i ◦ F ( i glyph[lscript] ) for each patch Ω ( i glyph[lscript] ) , glyph[lscript] = 0 , . . . , ν i -1, by using appropriate linear combinations of functions g ( i glyph[lscript] ) Γ ( i glyph[lscript] ) ; j 1 ,j 2 and g ( i glyph[lscript] ) Γ ( i glyph[lscript] +1 ) ; j 1 ,j 2 , j 1 = 0 1 2, , , j 2 = 0 1 , , . . . , 4 -j 1 , to ensure C 2 -smoothness across the interfaces Γ ( i glyph[lscript] ) and Γ ( i glyph[lscript] +1 ) , respectively, and by subtracting a linear combination of those standard B-splines N p r , j 1 ,j 2 , j , j 1 2 = 0 1 2, , , which have been added twice, to get a well defined C 2 -smooth function φ v ( ) i , cf. [27, Section 4.4].

In [27], the vertex subspace W v ( ) i has been defined as the span of all functions (9) whose spline functions g i glyph[lscript] fulfill the homogeneous linear system (11). A drawback of the resulting space is its geometry-dependent dimension. Below, we will construct a vertex subspace W v ( ) i , whose dimension will be independent of the initial geometry. This will be achieved by additionally enforcing that the functions of the generated space W v ( ) i have to be C 4 -smooth at the vertex v ( ) i . This strategy can be seen as an extension of approach [23], where C 1 -smooth functions have been generated in the vicinity of a vertex v ( ) i by additionally enforcing C 2 -smoothness of the functions at the vertex.

Let ψ j 1 ,j 2 : Ω → R , j , j 1 2 = 0 , . . . , 4 with j 1 + j 2 ≤ 4 be functions which are C 2 -smooth on Ω and even C 4 -smooth at the vertex v ( ) i such that

$$\partial _ { x _ { 1 } } ^ { m _ { 1 } } \partial _ { x _ { 2 } } ^ { m _ { 2 } } \psi _ { j _ { 1 }, j _ { 2 } } ( \mathfrak { v } ^ { ( i ) } ) = \sigma ^ { j _ { 1 } + j _ { 2 } } \delta _ { j _ { 1 } } ^ { m _ { 1 } } \delta _ { j _ { 2 } } ^ { m _ { 2 } },$$

where δ m j is the Kronecker delta, and σ is a scaling factor given by

$$\sigma = \left ( \frac { h } { p \, \nu _ { i } } \sum _ { \ell = 0 } ^ { \nu _ { i } - 1 } | | \nabla F ^ { ( i _ { \ell } ) } ( 0 ) | | \right ) ^ { - 1 },$$

cf. [23, Definition 14]. We define the isogeometric function φ v ( ) i I ; j 1 ,j 2 , j , j 1 2 = 0 1 , , . . . , 4, j 1 + j 2 ≤ 4, as

$$\phi _ { v _ { I } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \phi _ { v ^ { ( i ) } } ( x ),$$

where the function φ v ( ) i is specified by fulfilling the interpolation conditions

$$\partial _ { x _ { 1 } } ^ { m _ { 1 } } \partial _ { x _ { 2 } } ^ { m _ { 2 } } \phi _ { v ^ { ( i ) } } ( v ^ { ( i ) } ) = \partial _ { x _ { 1 } } ^ { m _ { 1 } } \partial _ { x _ { 2 } } ^ { m _ { 2 } } \psi _ { j _ { 1 }, j _ { 2 } } ( v ^ { ( i ) } ), \ m _ { 1 }, m _ { 2 } = 0, 1, \dots, 4, \ m _ { 1 } + m _ { 2 } \leq 4. \quad \ ( 1 3 )$$

The coefficients a Γ ( i glyph[lscript] ) j 1 ,j 2 and a ( i glyph[lscript] ) j 1 ,j 2 of the spline functions g i glyph[lscript] of the isogeometric function φ v ( ) i I ; j 1 ,j 2 are uniquely determined by the conditions (13), and can be computed via the following equivalent interpolation problem

$$\partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \, g _ { i _ { \ell } } ^ { ( i _ { \ell } ) } ( 0 ) & = \partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \left ( \psi _ { j _ { 1 }, j _ { 2 } } \circ \, F ^ { ( i _ { \ell } ) } \right ) ( 0 ), \quad 0 \leq m _ { 1 } \leq 4, \, 0 \leq m _ { 2 } \leq 2, \, m _ { 1 } + m _ { 2 } \leq 4, \\ \partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \, g _ { i _ { \ell } } ^ { 0 ^ { ( i _ { \ell } ) } } ( 0 ) & = \partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \left ( \psi _ { j _ { 1 }, j _ { 2 } } \circ \, F ^ { ( i _ { \ell } ) } \right ) ( 0 ), \quad 0 \leq m _ { 1 }, m _ { 2 } \leq 2,$$

for glyph[lscript] = 0 1 , , . . . , ν i -1. The presented construction of the isogeometric function φ v ( ) i I ; j 1 ,j 2 works and leads to a well-defined C 2 -smooth function on Ω, which is additionally C 4 -smooth at the vertex v ( ) i . This is true because of the following four reasons: First, since the functions g ( i glyph[lscript] ) Γ ( i glyph[lscript] ) ;0 ,j 2 , j 2 = 0 1 , , . . . , 4, g ( i glyph[lscript] ) Γ ( i glyph[lscript] ) ;1 ,j 2 , j 2 = 0 1 2 3, and , , , g ( i glyph[lscript] ) Γ ( i glyph[lscript] ) ;2 ,j 2 , j 2 = 0 1 2, of , , g Γ ( i glyph[lscript] ) i glyph[lscript] interpolate the traces of the form N p,r +2 j 2 , the first derivatives of the form N p -1 ,r +1 j 2 and the second derivatives of the form N p -2 , j 2 , respectively, across the interface Γ ( i glyph[lscript] ) , cf. [26, Section 7], the first set of interpolation conditions in (14) uniquely determine all coefficients a Γ ( i glyph[lscript] ) j 1 ,j 2 . Moreover, as the functions g Γ ( i glyph[lscript] ) i glyph[lscript] ensure C 2 -smoothness across the interfaces Γ ( i glyph[lscript] ) , we also obtain

$$\partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \, g _ { i _ { \ell } } ^ { \Gamma ^ { ( i _ { \ell + 1 } ) } } ( 0 ) = \partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \left ( \psi _ { j _ { 1 }, j _ { 2 } } \circ \mathcal { F } ^ { ( i _ { \ell } ) } \right ) ( 0 ),$$

for 0 ≤ m 1 ≤ 2, 0 ≤ m 2 ≤ 4, m 1 + m 2 ≤ 4 and glyph[lscript] = 0 1 , , . . . , ν i -1. Second, all coefficients a ( i glyph[lscript] ) j 1 ,j 2 are specified by the second set of interpolation conditions in (14). Third, thanks to the fulfillment of conditions (14) and (15) and the fact that the function g i glyph[lscript] is given by (10), we get

$$\partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \, g _ { i _ { \ell } } ( 0 ) = \partial _ { \xi _ { 1 } } ^ { m _ { 1 } } \partial _ { \xi _ { 2 } } ^ { m _ { 2 } } \left ( \psi _ { j _ { 1 }, j _ { 2 } } \circ \mathbf F ^ { ( i _ { \ell } ) } \right ) ( 0 ),$$

for 0 ≤ m 1 ≤ 4, 0 ≤ m 2 ≤ 4, m 1 + m 2 ≤ 4 and glyph[lscript] = 0 1 , , . . . , ν i -1, which is further equivalent to satisfying conditions (13). Moreover, we obtain that the function φ v ( ) i I ; j 1 ,j 2 is well defined and C 2 -smooth on Ω, where the latter property is the direct consequence of satisfying the homogeneous linear system (11). Fourth, the C 4 -smoothness of the function φ v ( ) i I ; j 1 ,j 2 at vertex v ( ) i follows directly from the fulfillment of the interpolation conditions (13).

Finally, we define for each inner vertex v ( ) i , i.e. i ∈ I Ξ I , the vertex subspace W v ( ) i as

$$\mathcal { W } _ { v ^ { ( i ) } } = \text{span} \{ \phi _ { v _ { l } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 1 }, j _ { 2 } = 0, 1, \dots, 4, \ j _ { 1 } + j _ { 2 } \leq 4 \}.$$

Let us continue with the case of a boundary vertex v ( ) i of patch valency ν i ≥ 3, that is i ∈ I Ξ 3 . This can be treated similarly to the case of an inner vertex as above. The vertex subspaces W v ( ) i , i ∈ I Ξ 3 , are defined as

$$\mathcal { W } _ { v ^ { ( i ) } } = \text{span} \{ \phi _ { v _ { 3 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 1 }, j _ { 2 } = 0, 1, \dots, 4, \ j _ { 1 } + j _ { 2 } \leq 4 \},$$

where the functions φ v ( ) i 3 ; j 1 ,j 2 are again determined just via

$$\phi _ { v _ { 3 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \phi _ { v ^ { ( i ) } } ( x ),$$

by interpolating the conditions (13). The only two differences in the construction of the functions are that the lower index glyph[lscript] of the indices i glyph[lscript] is not taken modulo ν i , and that for the patches Ω ( i 0 ) and Ω ( i glyph[lscript] -1 ) the functions g ( i 0 ) Γ ( i 0 ) ; j 1 ,j 2 and g ( i ν i -1 ) Γ ( i ν i ) ; j 1 ,j 2 in (10), which are given on the boundary edges Γ ( i 0 ) or Γ ( i ν i ) , are just defined as standard B-splines, i.e.

$$g ^ { ( i _ { 0 } ) } _ { \Gamma ^ { ( i _ { 0 } ) } ; j _ { 1 }, j _ { 2 } } ( \xi _ { 2 }, \xi _ { 1 } ) = N ^ { p, r } _ { j _ { 1 }, j _ { 2 } } ( \xi _ { 2 }, \xi _ { 1 } ) \quad \text{and} \quad g ^ { ( i _ { \nu _ { i } - 1 } ) } _ { \Gamma ^ { ( i _ { \nu _ { i } - 1 } ) } ; j _ { 1 }, j _ { 2 } } ( \xi _ { 1 }, \xi _ { 2 } ) = N ^ { p, r } _ { j _ { 1 }, j _ { 2 } } ( \xi _ { 1 }, \xi _ { 2 } ).$$

Analogously to the case of an inner vertex, one can argue that the functions φ v ( ) i 3 ; j 1 ,j 2 are C 2 -smooth on Ω and even C 4 -smooth at the vertex v ( ) i .

In case of a boundary vertex v ( ) i of patch valency ν i = 2, i.e. i ∈ I Ξ 2 , we assume for the two neighboring patches Ω ( i 0 ) and Ω ( i 1 ) enclosing the vertex v ( ) i , that the corresponding geometry mappings F ( i 0 ) and F ( i 1 ) are parameterized as shown in Fig. 1 (left), that we have v ( ) i = F ( i 0 ) ( 0 ) = F ( i 1 ) ( 0 ), and that the common interface is denoted by Γ ( j 0 ) , i.e. Γ ( j 0 ) = Ω ( i 0 ) ∩ Ω ( i 1 ) . Then, we construct the vertex subspace W v ( ) i as

$$\mathcal { W } _ { v ^ { ( i ) } } = \text{span} \{ \phi _ { v _ { 2 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 1 } = 0, 1, \dots, 4, \ j _ { 2 } = 0, 1, \dots, 4 - \min ( j _ { 1 }, 2 ) \},$$

where φ v ( ) i 2 ; j 1 ,j 2 is given as

$$\phi _ { v _ { 2 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \begin{cases} \phi _ { \Gamma _ { I } ^ { ( j _ { 0 } ) } ; j _ { 1 }, j _ { 2 } } ( x ) \quad \text{if $j_{1}=0,1,2,$} \\ \phi _ { \Omega ^ { ( i _ { 0 } ) } ; j _ { 1 } + \lfloor \frac { j _ { 2 } } { 2 } \rfloor, j _ { 2 } \bmod 2 } ( x ) \quad \text{if $j_{1}=3,$} \\ \phi _ { \Omega ^ { ( i _ { 1 } ) } ; j _ { 1 } - 1 + \lfloor \frac { j _ { 2 } } { 2 } \rfloor, j _ { 2 } \bmod 2 } ( x ) \quad \text{if $j_{1}=4,$} \end{cases}$$

with the functions φ Ω ( i glyph[lscript] ) ; j 1 ,j 2 and φ Γ ( j 0 ) I ; j 1 ,j 2 as defined in (5) and (7), respectively. The functions φ v ( ) i 2 ; j 1 ,j 2 are C 2 -smooth on Ω, since the functions φ Γ ( j 0 ) I ; j 1 ,j 2 , j 1 = 0 1 2, , , j 2 = 0 , . . . , 4 -j 1 , are C 2 -smooth on Ω by construction, and the functions φ Ω ( i glyph[lscript] ) ; j 1 ,j 2 , j 1 = 3 4, , j 2 = 0 1 , , . . . 1 -j 1 +3, glyph[lscript] = 0 1, possess vanishing values, gradients and Hessians along all , interfaces Γ ( j ) , j ∈ I Γ I

Finally, for the case of a boundary vertex v ( ) i of patch valency ν i = 1, i.e. i ∈ I Ξ 1 , and assuming without loss of generality that the vertex v ( ) i is given by v ( ) i = F ( i 0 ) ( 0 ), the vertex subspace W v ( ) i is of the form

$$\mathcal { W } _ { v ^ { ( i ) } } = \text{span} \{ \phi _ { v _ { 1 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } | \ j _ { 1 }, j _ { 2 } = 0, 1, \dots, 4, \ j _ { 1 } + j _ { 2 } \leq 4 \},$$

where φ v ( ) i 1 ; j 1 ,j 2 is defined as

$$\phi _ { v _ { 1 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } } ( x ) = \begin{cases} ( N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } \circ ( F ^ { ( i _ { 0 } ) } ) ^ { - 1 } ) ( x ) & \text{if } x \in \overline { \Omega ^ { ( i _ { 0 } ) } }, \\ 0 & \text{otherwise}. \end{cases}$$

The functions φ v ( ) i 1 ; j 1 ,j 2 are trivially C 2 -smooth on Ω, since the functions have vanishing values, gradients and Hessians along all interfaces Γ ( j ) , j ∈ I Γ I .

One could generate for the case of a boundary vertex v ( ) i of patch valency ν i = 2 or ν i = 1 a different vertex subspace W v ( ) i as described above by just following also for these cases the strategy explained for a boundary vertex of patch valency ν i ≥ 3. However, we will use in our numerical examples, see Section 4, the construction presented above, since it is simple and does not require the solving of a system of linear equations.

## 2.4. Basis and dimension of the space W

Since the space W is the direct sum (4), a basis of the space W can be obtained by just constructing bases of the individual subspaces W Ω ( ) i , i ∈ I Ω , W Γ i , i ∈ I Γ , and W v ( ) i , i ∈ I Ξ . As for each collection of functions (5), (6), (7), (12), (16), (17) or (18), spanning one of the subspaces W Ω ( ) i , W Γ ( ) i or W v ( ) i , the single functions are linearly independent by construction, the functions form a basis of the corresponding subspace. Moreover, the intersection of any two of the individual subspaces does not have any common function except the zero function. Therefore, the entire collection of all functions, i.e.

$$\text{opt the zero function. } I \text{ herefore, the entire collection of all functions, i.e.} \\ \phi _ { \Omega ^ { ( i ) } _ { B } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 }, j _ { 2 } = 3, 4, \dots, n - 4, \, i \in \mathcal { I } _ { \Omega } ; \\ \phi _ { \Gamma ^ { ( i ) } _ { B } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 } = 0, 1, 2, \, j _ { 2 } = 5 - j _ { 1 }, \dots, n + j _ { 1 } - 6, \, i \in \mathcal { I } _ { \Gamma _ { B } } ; \\ \phi _ { \Gamma ^ { ( i ) } _ { 1 } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 } = 0, 1, 2, \, j _ { 2 } = 5 - j _ { 1 }, \dots, n _ { j _ { 1 } } + j _ { 1 } - 6, \, i \in \mathcal { I } _ { \Gamma _ { I } } ; \\ \phi _ { v ^ { ( i ) } _ { 1 } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 }, j _ { 2 } = 0, 1, \dots, 4, \, j _ { 1 } + j _ { 2 } \leq 4, \, i \in \mathcal { I } _ { \Xi _ { I } } ; \\ \phi _ { v ^ { ( i ) } _ { 2 } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 } = 0, 1, \dots, 4, \, j _ { 2 } = 0, 1, \dots, 4 - \min ( j _ { 1 }, 2 ), \, i \in \mathcal { I } _ { \Xi _ { 2 } } ; \\ \dots \\ \Gamma _ { i } \quad \colon \, \colon \quad \omega \quad \alpha \quad \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \end{array}$$

and

$$\phi _ { v _ { 1 } ^ { ( i ) } ; j _ { 1 }, j _ { 2 } }, \, j _ { 1 }, j _ { 2 } = 0, 1, \dots, 4, \, j _ { 1 } + j _ { 2 } \leq 4, \, i \in \mathcal { I } _ { \Xi _ { 1 } } ;$$

builds a basis of the space W , cf. [27, Theorem 1].

In addition, the decomposition (4) of the space W implies that

$$\dim \mathcal { W } = \sum _ { i \in \mathcal { I } _ { \Omega } } \dim \mathcal { W } _ { \Omega ^ { ( i ) } } + \sum _ { i \in \mathcal { I } _ { \Gamma } } \dim \mathcal { W } _ { \Gamma ^ { ( i ) } } + \sum _ { i \in \mathcal { I } _ { \Xi } } \dim \mathcal { W } _ { v ^ { ( i ) } }.$$

Since the dimensions of the single subspaces W Ω ( ) i , i ∈ I Ω , W Γ i , i ∈ I Γ , and W v ( ) i , i ∈ I Ξ , are given by

$$\dim \mathcal { W } _ { \Omega ^ { ( i ) } } = ( n - 6 ) ^ { 2 },$$

$$\dim \mathcal { W } _ { \Gamma ^ { ( i ) } } = \begin{cases} 3 ( n - 8 ) \text{ if } i \in \mathcal { I } _ { \Gamma _ { B } }, \\ 3 ( n - 2 k - 9 ) \text{ if } i \in \mathcal { I } _ { \Gamma _ { I } }, \end{cases}$$

and

$$\dim \mathcal { W } _ { v ^ { ( i ) } } = \begin{cases} 1 5 \text{ if } i \in \mathcal { I } _ { \Xi _ { I } } \cup \mathcal { I } _ { \Xi _ { 1 } } \cup \mathcal { I } _ { \Xi _ { 3 } }, \\ 1 8 \text{ if } i \in \mathcal { I } _ { \Xi _ { 2 } }, \end{cases}$$

respectively, the dimension of W is equal to

$$\dim \mathcal { W } = | \mathcal { I } _ { \Omega } | ( n - 6 ) ^ { 2 } + 3 | \mathcal { I } _ { \Gamma ^ { I } } | ( n - 2 k - 9 ) + 3 | \mathcal { I } _ { \Gamma ^ { B } } | ( n - 8 ) + 1 5 \left ( | \mathcal { I } _ { \Xi _ { I } } | + | \mathcal { I } _ { \Xi _ { 1 } } | + | \mathcal { I } _ { \Xi _ { 3 } } | ) + 1 8 | \mathcal { I } _ { \Xi _ { 2 } } |. \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

Therefore, in contrast to [27], the dimension of the C 2 -smooth space W is independent of the initial geometry.

## 3. Isogeometric collocation

We will present the idea of using the concept of isogeometric collocation for solving the Poisson's equation on planar multi-patch domains.

## 3.1. Problem statement and isogeometric formulation

Let f : Ω → R , f 1 : ∂ Ω → R . We are interested in finding u : Ω → R , u ∈ C 2 (Ω), which solves the Poisson's equation

$$\triangle u ( x ) & = f ( x ), \quad x \in \Omega, \\ u ( x ) & = f _ { 1 } ( x ), \quad x \in \partial \Omega.$$

The idea is to employ isogeometric collocation to compute a C 2 -smooth approximation u h ∈ W of the solution u . This requires the use of global collocation points y j , j ∈ J , which are separated into two distinct sets of collocation points, namely of inner collocation points y I j , j ∈ J I , and of boundary collocation points y B j , j ∈ J B . Inserting these points into (20), we obtain

$$\triangle u ( \boldsymbol y _ { j } ^ { I } ) & = f ( \boldsymbol y _ { j } ^ { I } ), \quad j \in \mathcal { J } _ { I }, \\ u ( \boldsymbol y _ { j } ^ { B } ) & = f _ { 1 } ( \boldsymbol y _ { j } ^ { B } ), \quad j \in \mathcal { J } _ { B }.$$

To use the isogeometric approach for solving problem (21), we have to express each global collocation point y j with respect to local coordinates via

$$\zeta _ { j } ^ { ( \delta ( j ) ) ; I } = \left ( \mathcal { F } ^ { ( \delta ( j ) ) } \right ) ^ { - 1 } ( y _ { j } ^ { I } ) \,, \, j \in \mathcal { J } _ { I }, \quad \text{and} \quad \zeta _ { j } ^ { ( \delta ( j ) ) ; B } = \left ( \mathcal { F } ^ { ( \delta ( j ) ) } \right ) ^ { - 1 } ( y _ { j } ^ { B } ) \,, \, j \in \mathcal { J } _ { B },$$

where

$$\delta \colon \mathcal { J } \rightarrow \mathcal { I } _ { \Omega }, \quad \delta ( j ) = \min \{ i \in \mathcal { I } _ { \Omega }, \, \mathbf y _ { j } \in \overline { \Omega ^ { ( i ) } } \}.$$

Having the local collocation points ζ ( δ j ( )); I j and ζ ( δ j ( )); B j , then equations (21) can be transformed into

$$\left | \det J ^ { ( \delta ( j ) ) } \left ( \zeta _ { j } ^ { ( \delta ( j ) ) ; I } \right ) \right | & \left ( \nabla \circ \left ( N ^ { ( \delta ( j ) ) } ( \xi ) \nabla ( u \circ \mathcal { F } ^ { ( \delta ( j ) ) } ) ( \xi ) \right ) \right ) \right | _ { \xi = \zeta _ { j } ^ { ( \delta ( j ) ) ; I } } = f \left ( \mathcal { F } ^ { ( \delta ( j ) ) } \left ( \zeta _ { j } ^ { ( \delta ( j ) ) ; I } \right ) \right ), \ j \in \mathcal { J } _ { I }, \\ u \left ( \mathcal { F } ^ { ( \delta ( j ) ) } \left ( \zeta _ { j } ^ { ( \delta ( j ) ) ; B } \right ) \right ) = f _ { 1 } \left ( \mathcal { F } ^ { ( \delta ( j ) ) } \left ( \zeta _ { j } ^ { ( \delta ( j ) ) ; B } \right ) \right ), \ j \in \mathcal { J } _ { B }, \\ \quad \ \. \quad \.$$

where J ( glyph[lscript] ) is the Jacobian of F ( glyph[lscript] ) and

$$N ^ { ( \ell ) } ( \xi ) = \left ( J ^ { ( \ell ) } ( \xi ) \right ) ^ { - T } \left ( J ^ { ( \ell ) } ( \xi ) \right ) ^ { - 1 } \left | \det J ^ { ( \ell ) } ( \xi ) \right |.$$

This gives rise to a linear system for the unknown coefficients c i of the approximation u h = ∑ i ∈I c φ i i , where I = { 0 1 , , . . . , dim W- } 1 , and { φ i } i ∈I is the basis of W presented in Section 2.4.

## 3.2. Selection of collocation points

The choice of the collocation points plays an important role in the stability and convergence behavior of the numerical solution. Below, we will first present the general concept of choosing these points, and will then describe in more detail two particular selections of the collocation points.

## 3.2.1. The general concept

The idea is to choose for each patch Ω ( ) i , i ∈ I Ω , global collocation points F ( ) i ( ζ j ), which are determined by the local collocation points

$$\zeta _ { j } = ( \zeta _ { j _ { 1 } }, \zeta _ { j _ { 2 } } ) \in [ 0, 1 ] ^ { 2 },$$

which are again composed of univariate collocation points ζ j ∈ [0 , 1]. In case that we obtain the same global collocation point for more than one patch, we keep the corresponding point just for the patch Ω ( ) i with the smallest index i , to avoid repetitions of points. Moreover, we split the resulting points into inner and boundary collocation points, cf. Section 3.1.

In the one-patch case, mainly two different types of (univariate) collocation points have been considered. First, standard collocation points such as Greville abscissae or Demko abscissae are used, see e.g. [2, 30, 36], and possess a convergence behavior with respect 3

3 Instead of studying the convergence behavior with respect to the L 2 , H 1 and H 2 norm, often equivalently, the L ∞ , W 1 , ∞ and W 2 , ∞ -norm is considered in the literature.

to the L 2 , H 1 and H 2 norm of orders O ( h p -1 ), O ( h p -1 ) and O ( h p -1 ) and of orders O ( h p ), O ( h p ) and O ( h p -1 ) for odd and even spline degree p , respectively. These convergence rates are suboptimal for the L 2 norm and in the case of odd spline degree p also for the H 1 norm in comparison to the Galerkin approach, where the rates are of orders O ( h p +1 ), O ( h p ) and O ( h p -1 ).

The second type of collocation points are the so-called superconvergent points [1, 15, 32], which provide in case of odd spline degree p better convergence rates than the Greville abscissae. In more detail, in case of odd spline degree, the convergence behavior for the approach [15] is of orders O ( h p ), O ( h p ) and O ( h p -1 ) with respect to L 2 , H 1 and H 2 norm, and the convergence rates for the methods [1, 32] are even of optimal orders, i.e. O ( h p +1 ), O ( h p ) and O ( h p -1 ) with respect to L 2 , H 1 and H 2 norm. The basic concept of the superconvergent points is to choose as collocation points those points (or more precisely approximations of those points), which are the roots of the Galerkin residual D u 2 ( -u h ) of the considered problem. Since the resulting set of superconvergent points is approximately twice as large as the number of degrees of freedom, different strategies have been proposed in [1, 15, 32] to use the superconvergent points as collocation points for solving the collocation problem (21). In [1], the entire set of points is used, which leads to an overdetermined linear system that is solved by means of the least-squares method. In contrast, in [15, 32], particular subsets of superconvergent points are selected to obtain the same number of collocation points as the number of degrees of freedom. While the superconvergent points are selected in [15] in an alternating way, the points are chosen in [32] in a clustered way.

Below, we will extend the use of Greville abscissae and of superconvergent points as collocation points to the case of planar multi-patch domains. One main difference will be that splines of maximal smoothness, i.e. r = p -1, as usually employed in the onepatch case, cannot be used for the multi-patch case. In the latter case, the underlying spline spaces S p,r h have to fulfill 2 ≤ r ≤ p -3 with p ≥ 5 to allow the construction of h -refineable C 2 -smooth spline spaces W , cf. Section 2.1 and Section 2.3. For the sake of simplicity, we will restrict ourselves to the cases ( p, r ) ∈ { (5 , 2) , (6 , 2) , (6 , 3) } . In addition, the number of collocation points will be slightly larger than the number of degrees of freedom, cf. Section 3.2.2 and Section 3.2.3. Therefore, the resulting linear system will be overdetermined and will be solved as in [1] by means of the least-squares method. As already observed in [1], we first have to solve those equations in (21) which have been determined by the boundary collocation points, i.e. the second set of equations, and then use the result to solve the remaining system of linear equations.

## 3.2.2. Greville abscissae

For each patch Ω ( ) i , i ∈ I Ω , the local collocation points ζ j , j ∈ { 0 1 , , . . . , n -1 } 2 , from (22) are chosen as the tensor-product Greville points

$$\zeta _ { j } ^ { p, r } = ( \zeta _ { j _ { 1 } } ^ { p, r }, \zeta _ { j _ { 2 } } ^ { p, r } )$$

with ζ p,r j = t p,r j +1 + ...t p,r j + p p , j ∈ { 0 1 , , . . . , n - } 1 . The resulting global collocation points F ( ) i ( ζ p r , j ) for a particular example of a bilinearly parameterized three-patch domain are shown in Fig. 2. As already mentioned in the previous section, the number of (global) collocation

Figure 2: Set of all mapped Greville points F ( ) i ( ζ (5 2) , j ) (without repetitions) for a bilinearly parameterized three-patch domain and k = 4, where the boundary collocation points and the collocation points on the interfaces are specified in yellow and red, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-4a542f01-d9ea-4227-abbc-c9dc410a7c45.jpg)

points (without repetitions of same points) is slightly larger than the number of basis functions of the C 2 -smooth space W . But the quotient of the number of collocation points and of the number of basis functions converges to 1 when k grows as the following example will demonstrate.

Example 1. Let Ω be a multi-patch domain with exactly one inner vertex of edge and patch valency ν , see e.g. the three-, five- and six-patch domain in Fig 6 (a)-(c). In addition, let p = 5 and r = 2, which implies n = 6 + 3 k . Since dim W = ν (9 k 2 + 21 k + 12) + 15 (compare formula (19)) and the number of collocation points (without repetitions), i.e. |J | , is equal to νn 2 -νn +1 = ν (9 k 2 +33 k +30) + 1, we obtain

$$\lim _ { k \to \infty } \frac { | \mathcal { J } |$$

Note that the same conclusion is not only true for ( p, r ) = (6 , 2) and ( p, r ) = (6 , 3) but also for any multi-patch domain, since the number of collocation points is just slightly larger than the number of basis functions along the edges and in the vicinity of inner vertices and boundary vertices of patch valencies greater or equal to two.

The numerical results in Section 4 will show that we will obtain for the multi-patch case the same convergence behavior as for the one-patch case, namely with respect to the L 2 , H 1 and H 2 norm, convergence rates of orders O ( h p -1 ), O ( h p -1 ) and O ( h p -1 ) and of orders O ( h p ), O ( h p ) and O ( h p -1 ) for odd and even spline degree p , respectively.

## 3.2.3. Superconvergent points

The superconvergent points (cf. [1, 15, 32]) are defined as the roots of the Galerkin residual D u 2 ( -u h ) of the considered problem. Estimates of the univariate superconvergent points can be computed by solving a simple 1D Poisson's equation such as

$$u ^ { \prime \prime } ( x ) = f ( x ), \ \ x \in$$

with some particular function f . The superconvergent points have been mainly considered so far just for splines of maximal smoothness, i.e., r = p -1. Since the C 2 -smooth spline spaces W require underlying spline spaces S p,r h with 2 ≤ r ≤ p -3 and p ≥ 5, cf. Section 2.1 and Section 2.3, we will generate these points for some of these spaces, namely for the cases ( p, r ) ∈ { (5 , 2) , (6 , 2) , (6 , 3) } . To first estimate the superconvergent points, we just follow the approach presented in [15] based on a particular 1D Poisson's problem (23). Further studying the asymptotic behavior of the resulting estimates of the superconvergent points for the different levels of refinement, we obtain that the superconvergent points on each knot span with respect to the reference interval [ -1 1] are given as the roots of the , polynomial 15 x 4 -12 x 2 + 1 for p = 5 and r = 2 and as the roots of the polynomial 99 x 5 -130 x 3 +31 x for p = 6 and r = 2 3. The roots of these polynomials, and hence the , locations of the superconvergent points on each knot span with respect to the reference interval [ -1 1] , are presented in Table 1. All superconvergent points for ( p, r ) = (5 2) , and ( p, r ) ∈ { (6 , 2) , (6 , 3) } are shown in Fig. 3 and 4, respectively. In case of ( p, r ) ∈

$$\begin{array} {$$

Table 1: Locations of the superconvergent points on each knot span w.r.t. the reference interval [ -1 1] , for the cases ( p, r ) ∈ { (5 2) , , (6 2) , , (6 3) , } .

{ (6 , 2) , (6 , 3) } , the knots of the underlying spline spaces S p,r h are superconvergent points, and therefore, we count these points in case of inner knots just once for two neighboring segments. In case of ( p, r ) = (5 , 2), since the boundary points of the whole domain are not contained in the set of all superconvergent points, they have to be added to the set, see Fig. 3, to be able to impose Dirichlet boundary conditions.

The differences ∆ ( p,r ) between the number of all superconvergent points and the dimension of the underlying spline space S p,r h (i.e. dim S p,r h = p + 1 + k p ( -r )) are given by

$$& \Delta _ { ( 5, 2 ) } = 4 ( k + 1 ) + 2 - ($$

Figure 3: All superconvergent points (black and white points) and clustered superconvergent points (only black points) for the case ( p, r ) = (5 , 2). The center part (in red) is given by the second row depending on the number k of different inner knots.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-f3742877-7ccf-4d94-8413-afdac2605390.jpg)

Figure 4: First row: All superconvergent points (black and white points) and clustered superconvergent points (only black points) for ( p, r ) = (6 3). , Second row: All superconvergent points (black and yellow points), which coincide with the clustered superconvergent points, for the case ( p, r ) = (6 2), where the , yellow points specify the two added Greville abscissae.

The idea is to select from the superconvergent points in a clustered way a subset of collocation points such that the number of collocation points is equal to the number of basis functions of S p,r h , and such that the selection of points does not ruin the convergence properties as in [32] for the case of splines of maximal smoothness. While in the case of ( p, r ) = (6 , 2), we have even to add two points to the set of superconvergent points, in the cases ( p, r ) = (5 , 2) and ( p, r ) = (6 , 3) we have to omit k and k -2 points, respectively. We will call the resulting sets of collocation points clustered superconvergent points similarly to [32].

Let us first consider the case ( p, r ) = (6 , 2). We choose as additional points for the set of superconvergent points just the second and the second last Greville abscissa. This strategy has been already used for selecting additional points close to the boundary for Dirichlet boundary problems in case of splines of maximal smoothness, cf. [15, 32]. Note that for ( p, r ) = (6 2) the set of all superconvergent points coincides with the set of clustered , superconvergent points.

Recall that for the cases ( p, r ) = (5 , 2) and ( p, r ) = (6 , 3) we have to omit k and k -2 points, respectively. We denote by S k the indices of the points which will be omitted. Furthermore, we can assume k ≥ 2. For ( p, r ) = (5 2), , we choose depending on the number k of different inner knots the set S k as

$$S _ { k } = \widetilde { S } _ { k } \cup \begin{$$

with

$$\widetilde { S } _ { k } = \{ 4, 4 k + 1 \} \cup \bigcup _ { i = 0 } ^ { \lfloor \frac { k - 6 } { 4 } \rfloor } \left ( \{ 9 + 8 i, 1 2 + 8 i \} \cup \{ 4 k + 5 - ( 9 + 8 i ), 4 k + 5 - ( 1 2 + 8 i ) \} \right ).$$

In case of ( p, r ) = (6 3), we omit all inner knots except the first and the last one, which , leads to

$$S _ { k } = \{ 9 + 4 i, \ i = 0, 1, \dots, k - 3 \}.$$

All superconvergent points as well as the clustered superconvergent points are shown in Fig. 3 and 4 for the cases ( p, r ) = (5 2) and ( , p, r ) ∈ { (6 , 2) , (6 , 3) } , respectively. The following example will demonstrate that these points possess for the one-patch case the same convergence behavior with respect to the L 2 , H 1 and H 2 norm as the sets of all superconvergent points and of the clustered superconvergent points in the case of splines of maximal smoothness. Moreover, it will be shown that also for the case of Greville abscissae we get the same convergence rates for the splines of maximal smoothness and for the splines with a reduced regularity.

Example 2. We perform isogeometric collocation on the one-patch domain in Fig. 5 (first row), which is just the bilinearly parameterized unit square [0 , 1] 2 , and compare the convergence behavior under mesh-refinement for different spline spaces and different sets of collocation points. We use underlying spline spaces S p,r h ([0 , 1] 2 ) with p = 5 6 and 2 , ≤ r ≤ p -3 with mesh sizes h = 1 5 , 1 10 , 1 20 , 1 40 for ( p, r ) = (5 , 2) and h = 1 4 , 1 8 , 1 16 , 1 32 for p = 6 and r = 2 3. , Fig. 5 shows the resulting relative L 2 , H 1 and H 2 errors for the exact solution

$$u ( x _ { 1 }, x _ { 2 } ) = \cos ( 4 x _ { 1 } - 2 ) \sin \left ( 4 x _ { 2 } - \frac { 2 } { 3 } \right ),$$

(see first row), by using the Greville points (second row), all superconvergent points (third row) and the clustered superconvergent points (fourth row) as collocation points. We get as in the case of splines of maximal regularity, see [1, 32], the same rates of convergence depending on the used spline degree p . That is for p = 5 in case of Greville points the same rates of order O ( h 4 ) in the L 2 , H 1 and H 2 norm, and in case of all or clustered superconvergent points the optimal rates of orders O ( h 6 ), O ( h 5 ) and O ( h 4 ) in the L 2 , H 1 and H 2 -norm, respectively. For p = 6 we have the same rates independent of one of the three sets of collocation points, that is O ( h 6 ), O ( h 6 ) and O ( h 5 ) in the L 2 , H 1 and H 2 norm, respectively.

In the case of multi-patch domains, see Section 4, the convergence rates for the sets of superconvergent points will be the same as in the one-patch case, except for one particular case, namely for ( p, r ) = (5 , 2) with respect to the L 2 norm, where the rate will be just of order O ( h p ) instead of O ( h p +1 ). In the numerical examples in Section 4, for the sake of brevity, we will restrict ourselves just to the case of clustered superconvergent points, since experiments using all superconvergent points have led to the same convergence rates, and the number of clustered superconvergent points is dramatically smaller than the number

Relative L 2 , H 1 and H 2 errors for collocation at clustered superconvergent points

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-999623cb-a992-4100-a3a5-8429430592c8.jpg)

Figure 5: Isogeometric collocation on the bilinearly parameterized unit-square using different sets of collocation points, cf. Example. 2.

of all superconvergent points for the cases ( p, r ) = (5 2) and ( , p, r ) = (6 3). , However, in the case of multi-patch domains, also for the clustered superconvergent points the number of points is larger than the number of degrees of freedom. More precisely, the number of clustered superconvergent points is equal to the number of Greville abscissae, and therefore again, the quotient of the number of collocation points and of the number of degrees of freedom converges to one when the number k of different inner knots grows. For solving the resulting overdetermined linear system, we use as for the Greville abscissae the leastsquares method.

## 4. Numerical examples

We will perform isogeometric collocation on several multi-patch domains by using the Greville points and the clustered superconvergent points as collocation points (cf. Section 3.2.2 and Section 3.2.3), and will study the convergence under mesh refinement. Since the considered sets of collocation points will lead to an overdetermined system of linear equations, we will use the least-squares approach to solve the system, see Section 3.2. Below, we will denote the C 2 -smooth space W for a specific mesh size h by W h .

Example 3. We consider the four bilinearly parameterized multi-patch domains (a)-(d) shown in Fig. 6 (first row). For all domains we generate a sequence of C 2 -smooth spaces W h with mesh sizes h = 1 5 , 1 10 , 1 20 , 1 40 for p = 5 and r = 2 and with mesh sizes h = 1 4 , 1 8 , 1 16 , 1 32 for p = 6 and r = 2 3. The generated , C 2 -smooth spaces are then used to perform isogeometric collocation on the different multi-patch domains for a right side function f obtained from the exact solutions u visualized in Fig. 6 (a)-(d), which are given for the single domains (a)(d) by

$$u _ { a } ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { 2 x _ { 1 } } { 3 } \right ) \sin \left ( \frac { 2 x _ { 2 } } { 3 } \right ),$$

and

$$u _ { a } ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { 2 x _ { 1 } } { 3 } \right ) \sin \left ( \frac { 2 x _ { 2 } } { 3 } \right ), \\ u _ { b } ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { x _ { 1 } } { 2 } \right ) \sin \left ( \frac { x _ { 2 } - 2 } { 2 } \right ), \\ u _ { c } ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { x _ { 1 } + 3 } { 2 } \right ) \sin \left ( \frac { x _ { 2 } - 1 } { 2 } \right ), \\ u _ { d } ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { x _ { 1 } } { 3 } \right ) \sin \left ( \frac { x _ { 2 } } { 3 } \right ), \\ \text{The resulting relative } L ^ { 2 }, \, H ^ { 1 } \text{ and } H ^ { 2 } \text{ errors by using Grev} \\ \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{$$

respectively. The resulting relative L 2 , H 1 and H 2 errors by using Greville points and clustered superconvergent points are shown in Fig. 7 and 8, respectively. In case of Greville points the estimated convergence rates are the same as for the one-patch instance in Example 2, that is for p = 5 the same rates of order O ( h 4 ) in the L 2 , H 1 and H 2 norm, and for p = 6 the rates of orders O ( h 6 ), O ( h 6 ) and O ( h 5 ) in the L 2 , H 1 and H 2 norm, respectively. In case of clustered superconvergent points, we observe for both spline degrees p = 5 6 convergence rates of orders , O ( h p ), O ( h p ) and O ( h p -1 ) in the L 2 , H 1 and H 2 norm, respectively, which means for p = 5 a reduction of the order by 1 with respect

to the L 2 norm compared to the one-patch instance in Example 2. Note that using all superconvergent points would lead to the same reduced rate for p = 5 as employing clustered superconvergent points, and therefore will not be shown for the sake of brevity.

Figure 6: Computational domains and exact solutions for performing isogeometric collocation in Example 3 (see also Fig. 7 and 8).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-725e2d3a-3f30-4b67-a216-5792ca78f6a7.jpg)

Example 4. We consider the bicubic three-patch spline geometry taken from [27, Example 4], which is also visualized in Fig. 9 (first row). Since the geometry is parameterized in such a way that it is bilinear close to the patch interfaces (cf. Appendix A for the spline control points of the geometry), our approach can be also applied to this more general multi-patch domain. We perform isogeometric collocation on this domain for the right side function f obtained from the exact solution

$$u ( x _ { 1 }, x _ { 2 } ) = - 4 \cos \left ( \frac { 2 x _ { 1 } } { 3 } \right ) \sin \left ( \frac { 2 x _ { 2 } } { 3 } \right ),$$

see Fig. 6 (first row), by employing C 2 -smooth spaces W h with mesh sizes h = 1 4 , 1 8 , 1 16 , 1 32 for p = 5 6 and , r = 2. Again, we compare the resulting relative L 2 , H 1 and H 2 errors by using Greville points, see Fig. 9 (second row), and clustered superconvergent points, see Fig. 9 (third row), for collocation, and observe the same rates of convergence as in Example 3.

## 5. Conclusion

We have presented a method for computing a globally C 2 -smooth approximation of the solution of the Poisson's equation over planar bilinearly parameterized multi-patch domains. Our technique is based on the concept of isogeometric collocation and on the use

Figure 7: Error plots w.r.t. the number of degrees of freedom (NDOF) of performing isogeometric collocation on the different multi-patch domains given in Fig. 6 (first row) for right side functions obtained by the exact solutions shown in Fig 6 (second row) using the Greville points as collocation points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-28a976d1-3515-4c00-8f2a-f64855809890.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-5da1af3f-835a-4c41-b8b5-5b7c762a3e2e.jpg)

Five-patch domain (d)

Figure 8: Error plots w.r.t. the number of degrees of freedom (NDOF) of performing isogeometric collocation on the different multi-patch domains given in Fig. 6 (first row) for right side functions obtained by the exact solutions shown in Fig 6 (second row) using the clustered superconvergent points as collocation points.

Relative L 2 , H 1 and H 2 errors for collocation at clustered superconvergent points

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00813v2-figure-8b1dece2-827e-4422-99cb-7c5e25ace9f7.jpg)

Figure 9: Isogeometric collocation on a bicubic three-patch spline geometry using different sets of collocation points, cf. Example. 4.

of a globally C 2 -smooth isogeometric spline space as discretization space. The constructed C 2 -smooth space is a particular subspace of the space of all globally C 2 -smooth isogeometric spline functions on the considered multi-patch domain, and can be generated as the direct sum of simpler spaces corresponding to the individual patches, edges and vertices of the multi-patch domain. Moreover, the C 2 -smooth space possesses a dimension independent of the initial geometry, and the construction of a basis with locally supported functions is simple and works uniformly for all possible multi-patch configurations.

Two different approaches for the choice of the collocation points have been proposed, where both strategies can be seen as the extension of successfully tested techniques in the one-patch case. On the one hand, we employ the tensor-product Greville points, and on the other hand, we generalize the concept of superconvergent points (cf. [1, 15, 32]) to the case of planar multi-patch domains. While in the one-patch case splines of maximal smoothness, i.e. r = p -1, are normally used, the C 2 -smooth multi-patch spline spaces require underlying spline spaces S p,r h of regularity 2 ≤ r ≤ p -3. For the sake of simplicity, we have restricted ourselves to the cases ( p, r ) ∈ { (5 , 2) , (6 , 2) , (6 , 3) } .

For both choices of the collocation points, the numerical results have shown the same convergence behavior with respect to L 2 , H 1 and H 2 norm as in the one-patch case except in one particular case for the superconvergent points. Namely, in case of ( p, r ) = (5 , 2), the convergence rates with respect to L 2 norm seems to be always reduced by one order, that is O ( h 5 ) instead of the possible order of O ( h 6 ) in the one-patch case. However, already this slightly reduced convergence behavior in the L 2 norm is better as the one for the case of Greville points, where the obtained convergence rates are just of order O ( h 4 ) with respect to the L 2 norm.

So far, our isogeometric collocation method is directly applicable only to elliptic PDEs of second order with Dirichlet boundary conditions. To impose also Neumann boundary conditions, further investigations amongst others in connection with the solving of the overdetermined linear system with respect to the fulfillment of the boundary conditions are needed, cf. [1, Remark 2.3.2]. Starting point of these studies could be e.g. the publications [11, 19, 37], where different strategies are described to deal also with Neumann boundary conditions.

The paper leaves several further open issues, which are worth to study. Since the number of collocation points is slightly larger than the dimension of the C 2 -smooth discretization space, the finding of a set of collocation points with the same cardinality as the dimension of the space is of interest to avoid the necessity of the least-squares method for solving the resulting linear system. Clearly, collocation points which lead to optimal convergence rates are of interest, too. However, this is also still an open problem for even spline degree in the one-patch case. A further possible topic could be the extension of our approach to more general multi-patch domains such as bilinear-like parameterizations, cf. [26], or to multi-patch structured shells and multi-patch volumes.

Acknowledgment. The authors wish to thank the anonymous reviewers for their comments that helped to improve the paper. V. Vitrih was partially supported by the Slovenian Research Agency (research program P1-0404 and research projects J1-9186, J1-1715). This support is gratefully acknowledged.

## Appendix A. Spline control points of the geometry from Example 4

The spline geometry from Fig. 9 (first row), which was taken from [27, Example 4] and which is also used in Example 4 in this paper, consists of three geometry mappings F ( ) i ∈ S p r , h ([0 , 1] 2 ) ×S p r , h ([0 , 1] 2 ), i = 1 2 3, with , , p = (3 3), , r = (2 2) and , h = 1 4. / Therefore, each geometry mapping F ( ) i , i = 1 2 3, possesses a spline parameterization of the form , ,

$$F ^ { ( i ) } ( \xi _ { 1 }, \xi _ { 2 } ) = \sum _ { j _ { 1 } = 0 } ^ { 6 } \sum _ { j _ { 2 } = 0 } ^ { 6 } c _ { j _ { 1 }, j _ { 2 } } ^ { ( i ) } N _ { j _ { 1 }, j _ { 2 } } ^ { p, r } ( \xi _ { 1 }, \xi _ { 2 } ),$$

where the single control points c ( ) i j 1 ,j 2 are given in Table A.2.

Table A.2: Control points c ( ) i j 1 ,j 2 , i = 1 2 3, of the bicubic three-patch spline geometry shown in Fig. 9 (first , , row).

| c (1) j ,j                                                                                                                      | c (1) j ,j                                                                                                                                                    | c (1) j ,j                                                                                                                                       | c (1) j ,j                                                                                                                                   | c (1) j ,j                                                                                                                                        | c (1) j ,j                                                                                                                                          | c (1) j ,j                                                                                                                      |
|---------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------|
| ( 17 3 , 2) ( 101 18 , 11 6 ) ( 11 2 , 3 2 ) ( 16 3 , 1) ( 31 6 , 1 2 ) ( 91 18 , 1 6 ) (5 , 0)                                 | ( 853 144 , 169 84 ) ( 10163 1728 , 1859 1008 ) ( 371 64 , 169 112 ) ( 1633 288 , 169 168 ) ( 3193 576 , 169 336 ) ( 9433 1728 , 169 1008 ) ( 65 12 , 0)      | ( 103 16 , 57 28 ) ( 411 64 , 209 112 ) ( 409 64 , 171 112 ) ( 203 32 , 57 56 ) ( 403 64 , 57 112 ) ( 401 64 , 19 112 ) ( 25 4 , 0)              | 1 2 ( 173 24 , 29 14 ) ( 2083 288 , 319 168 ) ( 233 32 , 87 56 ) ( 37 5 , 99 100 ) ( 749 100 , 8 25 ) ( 15 2 , - 21 100 ) ( 15 2 , - 12 25 ) | ( 383 48 , 59 28 ) ( 4633 576 , 649 336 ) ( 523 64 , 177 112 ) ( 849 100 , 103 100 ) ( 873 100 , 29 100 ) ( 437 50 , - 12 25 ) ( 35 4 , - 21 25 ) | ( 1223 144 , 179 84 ) ( 14833 1728 , 1969 1008 ( 561 64 , 179 112 ) ( 923 100 , 57 50 ) ( 961 100 , 27 50 ) ( 961 100 , - 6 25 ) ( 115 12 , - 2 3 ) | ( 35 4 , 15 7 ) ) ( 425 48 , 55 28 ) ( 145 16 , 45 28 ) ( 48 5 , 121 100 ) ( 251 25 , 13 20 ) ( 201 20 , - 1 14 ) (10 , - 1 2 ) |
| c (2) j 1 ,j 2                                                                                                                  | c (2) j 1 ,j 2                                                                                                                                                | c (2) j 1 ,j 2                                                                                                                                   | c (2) j 1 ,j 2                                                                                                                               | c (2) j 1 ,j 2                                                                                                                                    | c (2) j 1 ,j 2                                                                                                                                      | c (2) j 1 ,j 2                                                                                                                  |
| ( 17 3 , 2) ( 853 144 , 169 84 ) ( 103 16 , 57 28 ) ( 173 24 , 29 14 ) ( 383 48 , 59 28 ) ( 1223 144 , 179 84 ) ( 35 4 , 15 7 ) | ( 50 9 , 13 6 ) ( 10033 1728 , 2209 1008 ) ( 1211 192 , 251 112 ) ( 2033 288 , 389 168 ) ( 4499 576 , 803 336 ) ( 14363 1728 , 2459 1008 ) ( 137 16 , 69 28 ) | ( 16 3 , 5 2 ) ( 3209 576 , 857 336 ) ( 387 64 , 297 112 ) ( 649 96 , 157 56 ) ( 1435 192 , 331 112 ) ( 4579 576 , 1027 336 ) ( 131 16 , 87 28 ) | (5 , 3) ( 167 32 , 173 56 ) ( 181 32 , 183 56 ) ( 127 20 , 359 100 ) ( 178 25 , 391 100 ) ( 769 100 , 409 100 ) ( 799 100 , 209 50 )         | ( 14 3 , 7 2 ) ( 2803 576 , 1219 336 ) ( 337 64 , 435 112 ) ( 587 100 , 441 100 ) ( 661 100 , 49 10 ) ( 73 10 , 257 50 ) ( 763 100 , 263 50 )     | ( 40 9 , 23 6 ) ( 8003 1728 , 4019 1008 ) ( 961 192 , 481 112 ) ( 109 20 , 99 20 ) ( 597 100 , 279 50 ) ( 33 5 , 583 100 ) ( 111 16 , 83 14 )       | ( 13 3 , 4) ( 325 72 , 25 6 ) ( 39 8 , 9 2 ) ( 523 100 , 523 100 ) ( 142 25 , 591 100 ) ( 31 5 , 37 6 ) ( 13 2 , 25 4 )         |
| c (3) j 1 ,j 2                                                                                                                  | c (3) j 1 ,j 2                                                                                                                                                | c (3) j 1 ,j 2                                                                                                                                   | c (3) j 1 ,j 2                                                                                                                               | c (3) j 1 ,j 2                                                                                                                                    | c (3) j 1 ,j 2                                                                                                                                      | c (3) j 1 ,j 2                                                                                                                  |
| ( 17 3 , 2) ( 50 9 , 13 6 ) ( 16 3 , 5 2 ) (5 , 3) ( 14 3 , 7 2 ) ( 40 9 , 23 6 ) ( 13 3 , 4)                                   | ( 101 18 , 11 6 ) ( 2365 432 , 143 72 ) ( 749 144 , 55 24 ) ( 115 24 , 11 4 ) ( 631 144 , 77 24 ) ( 1775 432 , 253 72 ) ( 143 36 , 11 3 )                     | ( 11 2 , 3 2 ) ( 85 16 , 13 8 ) ( 79 16 , 15 8 ) ( 35 8 , 9 4 ) ( 61 16 , 21 8 ) ( 55 16 , 23 8 ) ( 13 4 , 3)                                    | ( 16 3 , 1) ( 365 72 , 13 12 ) ( 109 24 , 5 4 ) ( 19 5 , 43 25 ) ( 61 20 , 227 100 ) ( 251 100 , 263 100 ) ( 56 25 , 141 50 )                | ( 31 6 , 1 2 ) ( 695 144 , 13 24 ) ( 199 48 , 5 8 ) ( 83 25 , 117 100 ) ( 123 50 , 181 100 ) ( 7 4 , 113 50 ) ( 141 100 , 247 100 )               | ( 91 18 , 1 6 ) ( 2015 432 , 13 72 ) ( 559 144 , 5 24 ) ( 153 50 , 7 10 ) ( 113 50 , 6 5 ) ( 151 100 , 149 100 ) ( 10 9 , 19 12 )                   | (5 , 0) ( 55 12 , 0) ( 15 4 , 0) ( 59 20 , 9 20 ) ( 43 20 , 9 10 ) ( 17 12 , 1) (1 , 1)                                         |

## References

| [1]   | C. Anitescu, Y. Jia, Y. J. Zhang, and T. Rabczuk. An isogeometric collocation method using superconvergent points. Comput. Methods Appl. Mech. Engrg. , 284:1073-1097, 2015.                                 |
|-------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [2]   | F. Auricchio, L. Beir˜o a da Veiga, T. J. R. Hughes, A. Reali, and G. Sangalli. Iso- geometric collocation methods. Math. Models Methods Appl. Sci. , 20(11):2075-2107, 2010.                                |
| [3]   | F. Auricchio, L. Beir˜o a da Veiga, T.J.R. Hughes, A. Reali, and G. Sangalli. Isoge- ometric collocation for elastostatics and explicit dynamics. Comput. Methods Appl. Mech. Engrg. , 249-252:2 - 14, 2012. |
| [4]   | L. Beir˜o a da Veiga, A. Buffa, G. Sangalli, and R. V´zquez. a Mathematical analysis of variational isogeometric methods. Acta Numerica , 23:157-287, 5 2014.                                                |
| [5]   | L. Beir˜o a da Veiga, C. Lovadina, and A. Reali. Avoiding shear locking for the Timo- shenko beam problem via isogeometric collocation methods. Comput. Methods Appl. Mech. Engrg. , 241-244:38 - 51, 2012.  |
| [6]   | M. Bercovier and T. Matskewich. Smooth B´zier e Surfaces over Unstructured Quadri- lateral Meshes . Lecture Notes of the Unione Matematica Italiana, Springer, 2017.                                         |
| [7]   | A. Blidia, B. Mourrain, and N. Villamizar. G 1 -smooth splines on quad meshes with 4-split macro-patch elements. Comput. Aided Geom. Des. , 52-53:106 - 125, 2017.                                           |
| [8]   | C.L. Chan, C. Anitescu, and T. Rabczuk. Isogeometric analysis with strong multipatch C 1 -coupling. Comput. Aided Geom. Design , 62:294-310, 2018.                                                           |
| [9]   | A. Collin, G. Sangalli, and T. Takacs. Analysis-suitable G 1 multi-patch parametriza- tions for C 1 isogeometric spaces. Comput. Aided Geom. Des. , 47:93 - 113, 2016.                                       |
| [10]  | J. A. Cottrell, T. J. R. Hughes, and Y. Bazilevs. Isogeometric Analysis: Toward Integration of CAD and FEA . John Wiley & Sons, Chichester, England, 2009.                                                   |
| [11]  | L. De Lorenzis, J. A. Evans, T. J. R Hughes, and A. Reali. Isogeometric collocation: Neumann boundary conditions and contact. Comput. Methods Appl. Mech. Engrg. , 284:21 - 54, 2015.                        |
| [12]  | J. Deng, F. Chen, X. Xin Li, C. Hu, W. Tong, Z. Yang, and Y. Feng. Polynomial splines over hierarchical T-meshes. Graphical Models , 70(4):76-86, 2008.                                                      |
| [13]  | J. A. Evans, R. R. Hiemstra, T. J. R. Hughes, and A. Reali. Explicit higher-order accurate isogeometric collocation methods for structural dynamics. Comput. Methods Appl. Mech. Engrg. , 338:208-240, 2018. |

| [14]   | F. Fahrendorf, L. De Lorenzis, and H. Gomez. Reduced integration at superconvergent points in isogeometric analysis. Comput. Methods Appl. Mech. Engrg. , 328:390-410, 2018.                                       |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [15]   | H. Gomez and L. De Lorenzis. The variational collocation method. Comput. Methods Appl. Mech. Engrg. , 309:152-181, 2016.                                                                                           |
| [16]   | D. Groisser and J. Peters. Matched G k -constructions always yield C k -continuous isogeometric elements. Comput. Aided Geom. Des. , 34:67 - 72, 2015.                                                             |
| [17]   | J. Hoschek and D. Lasser. Fundamentals of computer aided geometric design . A K Peters Ltd., Wellesley, MA, 1993.                                                                                                  |
| [18]   | T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput. Methods Appl. Mech. Engrg. , 194(39-41):4135-4195, 2005.         |
| [19]   | Y. Jia, C. Anitescu, Y. J. Zhang, and T. Rabczuk. An adaptive isogeometric analysis collocation method with a recovery-based error estimator. Comput. Methods Appl. Mech. Engrg. , 345:52-74, 2019.                |
| [20]   | M. Kapl, F. Buchegger, M. Bercovier, and B. J¨ttler. u Isogeometric analysis with geo- metrically continuous functions on planar multi-patch geometries. Comput. Methods Appl. Mech. Engrg. , 316:209 - 234, 2017. |
| [21]   | M. Kapl, G. Sangalli, and T. Takacs. Dimension and basis construction for analysis- suitable G 1 two-patch parameterizations. Comput. Aided Geom. Des. , 52-53:75 - 89, 2017.                                      |
| [22]   | M. Kapl, G. Sangalli, and T. Takacs. Isogeometric analysis with C 1 functions on unstructured quadrilateral meshes. Technical Report 1812.09088, arXiv.org, 2018.                                                  |
| [23]   | M. Kapl, G. Sangalli, and T. Takacs. An isogeometric C 1 subspace on unstructured multi-patch planar domains. Comput. Aided Geom. Des. , 69:55-75, 2019.                                                           |
| [24]   | M. Kapl and V. Vitrih. Space of C 2 -smooth geometrically continuous isogeometric functions on planar multi-patch geometries: Dimension and numerical experiments. Comput. Math. Appl. , 73(10):2319-2338, 2017.   |
| [25]   | M. Kapl and V. Vitrih. Space of C 2 -smooth geometrically continuous isogeometric functions on two-patch geometries. Comput. Math. Appl. , 73(1):37 - 59, 2017.                                                    |
| [26]   | M. Kapl and V. Vitrih. Dimension and basis construction for C 2 -smooth isogeometric spline spaces over bilinear-like G 2 two-patch parameterizations. J. Comput. Appl. Math. , 335:289-311, 2018.                 |
| [27]   | M. Kapl and V. Vitrih. Solving the triharmonic equation over multi-patch planar domains using isogeometric analysis. J. Comput. Appl. Math. , 358:385-404, 2019.                                                   |

| [28]   | M. Kapl, V. Vitrih, B. J¨ttler, u and K. Birner. Isogeometric analysis with geometrically continuous functions on two-patch geometries. Comput. Math. Appl. , 70(7):1518 - 1538, 2015.                                                                            |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [29]   | K. Karˇiauskas, c T. Nguyen, and J. Peters. Generalizing bicubic splines for modeling and IGA with irregular layout. Comput.-Aided Des. , 70:23 - 35, 2016.                                                                                                       |
| [30]   | J. Kiendl, E. Marino, and L. De Lorenzis. Isogeometric collocation for the Reissner- Mindlin shell problem. Comput. Methods Appl. Mech. Engrg. , 325:645-665, 2017.                                                                                               |
| [31]   | E. Marino, J. Kiendl, and L. De Lorenzis. Explicit isogeometric collocation for the dynamics of three-dimensional beams undergoing finite motions. Comput. Methods Appl. Mech. Engrg. , 343:530-549, 2019.                                                        |
| [32]   | M. Montardini, G. Sangalli, and L. Tamellini. Optimal-order isogeometric collocation at Galerkin superconvergent points. Comput. Methods Appl. Mech. Engrg. , 316:741- 757, 2017.                                                                                 |
| [33]   | B. Mourrain, R. Vidunas, and N. Villamizar. Dimension and bases for geometrically continuous splines on surfaces of arbitrary topology. Comput. Aided Geom. Des. , 45:108 - 133, 2016.                                                                            |
| [34]   | T. Nguyen and J. Peters. Refinable C 1 spline elements for irregular quad layout. Comput. Aided Geom. Des. , 43:123 - 130, 2016.                                                                                                                                  |
| [35]   | J. Peters. Geometric continuity. In Handbook of computer aided geometric design , pages 193-227. North-Holland, Amsterdam, 2002.                                                                                                                                  |
| [36]   | A. Reali and T. J. R. Hughes. An introduction to isogeometric collocation methods. In Isogeometric Methods for Numerical Simulation , pages 173-204. Springer, 2015.                                                                                              |
| [37]   | D. Schillinger, M. J. Borden, and H. K. Stolarski. Isogeometric collocation for phase- field fracture models. Comput. Methods Appl. Mech. Engrg. , 284:583 - 610, 2015.                                                                                           |
| [38]   | D. Schillinger, J. A. Evans, A. Reali, M. A. Scott, and T. J.R. Hughes. Isogeometric collocation: Cost comparison with Galerkin methods and extension to adaptive hi- erarchical NURBS discretizations. Comput. Methods Appl. Mech. Engrg. , 267:170 - 232, 2013. |
| [39]   | D. Toshniwal, H. Speleers, R. Hiemstra, and T. J. R. Hughes. Multi-degree smooth polar splines: Aframework for geometric modeling and isogeometric analysis. Comput. Methods Appl. Mech. Engrg. , 2016.                                                           |
| [40]   | D. Toshniwal, H. Speleers, and T. J. R. Hughes. Analysis-suitable spline spaces of arbitrary degree on unstructured quadrilateral meshes. Technical Report 16, Institute for Computational Engineering and Sciences (ICES), 2017.                                 |

- [41] D. Toshniwal, H. Speleers, and T. J. R. Hughes. Smooth cubic spline spaces on unstructured quadrilateral meshes with particular emphasis on extraordinary points: Geometric design and isogeometric analysis considerations. Comput. Methods Appl. Mech. Engrg. , 327:411-458, 2017.
Mario Kapl Vito Vitrih
Isogeometric collocation on planar multi-patch domains
We present an isogeometric framework based on collocation to construct a
$C^2$-smooth approximation of the solution of the Poisson's equation over
planar bilinearly parameterized multi-patch domains. The construction of the
used globally $C^2$-smooth discretization space for the partial differential
equation is simple and works uniformly for all possible multi-patch
configurations. The basis of the $C^2$-smooth space can be described as the
span of three different types of locally supported functions corresponding to
the single patches, edges and vertices of the multi-patch domain. For the
selection of the collocation points, which is important for the stability and
convergence of the collocation problem, two different choices are numerically
investigated. The first approach employs the tensor-product Greville abscissae
as collocation points, and shows for the multi-patch case the same convergence
behavior as for the one-patch case [2], which is suboptimal in particular for
odd spline degree. The second approach generalizes the concept of
superconvergent points from the one-patch case (cf. [1, 15, 32]) to the
multi-patch case. Again, these points possess better convergence properties
than Greville abscissae in case of odd spline degree.