Pith. sign in

REVIEW 3 major objections 4 minor 41 references

Isogeometric collocation on planar multi-patch domains

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read The paper constructs a globally C2-smooth isogeometric spline space whose dimension is independent of the geometry, and uses collocation on it to solve the Poisson equation over planar bilinear multi-patch domains.

desk verdict A useful, honest C^2 multi-patch collocation paper whose main construction is plausible and well tested, but whose dimension/basis claim needs proof before I'd fully trust it. read the letter →

arxiv 1908.00813 v2 pith:WACWMRFH submitted 2019-08-02 math.NA cs.NA

classification math.NAcs.NA MSC 65N3565D1768U07
keywords isogeometricanalysiscollocationC2-smoothmulti-patchsplinespacesuperconvergentpointsPoisson'sequationGrevilleabscissaedomain
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to solve the Poisson equation on planar domains split into several bilinearly parameterized patches by collocation, using a globally C2-smooth spline space as the discretization space. The construction splits the space into patch, edge, and vertex parts and enforces extra smoothness at vertices so that the space dimension is independent of the initial multi-patch geometry. Two collocation point strategies are compared: tensor-product Greville abscissae and clustered superconvergent points. The numerical experiments indicate that the method reproduces one-patch convergence behavior, with superconvergent points improving the odd-degree case.

What carries the argument

The load-bearing object is the direct-sum decomposition W = (sum of patch subspaces) ⊕ (sum of edge subspaces) ⊕ (sum of vertex subspaces) of the C2-smooth space, with each subspace spanned by locally supported functions. C2-smoothness across interfaces is enforced through the geometric-continuity conditions (1)-(3), expressed via the functions α and β that describe how adjacent bilinear parameterizations meet. The vertex subspaces are modified by requiring C4-smoothness at the vertex, which fixes their dimension to 15 (or 18 for valency-2 boundary vertices) and makes the total dimension (19) geometry-independent. The collocation system is assembled by mapping either Greville abscissae or clustered superconvergent points through each patch geometry, leading to an overdetermined linear system solved by least squares.

What would settle it

Compute, for a small concrete configuration (for example the three-patch domain of Fig. 6(a) with p=5, r=2, k=2), the actual rank of the matrix whose columns are the evaluations of the candidate basis functions at a sufficiently large set of sample points; if the rank is strictly less than the value from (19), the basis and direct-sum claim fails. Alternatively, rerun Example 3 with a manufactured solution and check whether the L2 error for p=5 with clustered superconvergent points decreases as O($h^{5}$); an O($h^{4}$) rate would indicate that the claimed convergence behavior is not reproduced.

Watch

Extended reading notes

Core claim

The paper's central claim is that for any planar bilinearly parameterized multi-patch domain one can construct, by explicit local formulas or small linear systems, a basis of locally supported C2-smooth isogeometric functions whose dimension is given by formula (19) and is independent of the geometry. The basis comprises functions associated with patches, edges, and vertices; vertex functions are forced to be C4 at the corresponding vertex, which removes the geometry-dependent degrees of freedom present in earlier constructions. With this space as ansatz, collocation for the Poisson equation with Dirichlet boundary conditions yields the reported convergence rates under h-refinement, and the same rates as the one-patch case are observed for the tested degrees p=5 and p=6 and regularities r=2,3.

Load-bearing premise

The load-bearing premise is that the collection of locally supported functions built in Section 2.3 is genuinely a basis: linear independence of each family and trivial pairwise subspace intersections are asserted rather than proved, and the same holds for the C4-smoothness of the vertex functions; without this, the dimension count (19) and the nonsingularity of the least-squares system are not guaranteed.

Editorial extensions

If this is right

  • A C2-smooth collocation solve on multi-patch domains can be assembled without special treatment of interfaces, unlike methods that only couple patches with C0 continuity.
  • For the tested spaces (p,r) in {(5,2),(6,2),(6,3)}, Greville points give the same suboptimal odd-degree rates as in the one-patch case, while clustered superconvergent points recover O(h^p) in L2 and H1 and O(h^{p-1}) in H2.
  • The dimension formula (19) lets users predict the number of degrees of freedom before solving, regardless of the patch layout.
  • The method applies to any planar bilinear multi-patch domain and also to a bicubic three-patch geometry that is bilinear close to the patch interfaces, since the construction only needs bilinear behavior near interfaces.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The C4-at-vertices trick is a general dimension-stabilization mechanism and could likely be reused to make C1 or C3 multi-patch spline spaces geometry-independent in the same fashion.
  • The polynomial roots for superconvergent points (for example 15x^4 - 12x^2 + 1 for p=5, r=2) may serve as building blocks for other reduced-regularity spline spaces, not only the three cases tested.
  • Since the collocation system is overdetermined, a natural next step not pursued in the paper is to search for a square subset of points that keeps the convergence rates, which would remove the least-squares solve.
  • The observed loss of one L2 order for p=5 with superconvergent points on multi-patch domains suggests that the vertex or interface coupling, rather than the point selection, is the limiting factor; testing a single-patch analog with the same reduced regularity would isolate the cause.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper presents an isogeometric collocation framework for the Poisson equation on planar bilinearly parameterized multi-patch domains. The discretization space is a globally C2-smooth spline space W defined as the direct sum of patch, edge, and vertex subspaces; the vertex construction is modified by enforcing C4-smoothness at vertices so that the space dimension is independent of the initial geometry (Eq. 19). Two choices of collocation points are investigated: tensor-product Greville abscissae and clustered superconvergent points. Numerical experiments on several multi-patch domains (Section 4) report convergence rates of O(h^p), O(h^p), and O(h^{p-1}) in the L2, H1, and H2 norms for the superconvergent points, with one observed reduction to O(h^5) for p=5 in the L2 norm.

Significance. If the construction is sound, the paper makes a useful contribution to isogeometric collocation for multi-patch domains: it provides a uniform, locally supported basis for a C2-smooth space whose dimension is independent of the geometry, and it extends superconvergent-point collocation from one-patch to multi-patch settings. The numerical results support the claimed convergence rates for the tested configurations, and the authors are explicit about the one suboptimal case. The main weakness is that the theoretical foundation of the space construction—the linear independence and direct-sum decomposition in Section 2.4—is asserted rather than proved, and this is exactly the part that is new relative to the authors' prior work [27].

major comments (3)
  1. [§2.4] The claim that the collection of functions (5)–(18) forms a basis of W, and the resulting dimension formula (19), rests on the assertion that the decomposition (4) is a direct sum and that the individual functions are linearly independent. This is not proved for the new vertex subspaces constructed in §2.3.2. In particular, the cited [27, Theorem 1] does not cover the newly added C4 vertex conditions or the interpolation systems (13)–(14). The authors need to show that the systems (14) are nonsingular for all admissible valencies and bilinear parameterizations, that the 15 (or 18) vertex functions are linearly independent, and that no vertex function is a nontrivial linear combination of edge functions with overlapping support near the vertex. Without this, the dimension formula (19) and the well-posedness of the collocation system (21) are not established.
  2. [§2.3.2] The sentence after Eq. (14) states that the first set of interpolation conditions 'uniquely determine all coefficients a^{Γ(i_l)}_{j1,j2}' and that the second set determines all a^{(i_l)}_{j1,j2}. No argument is given for the count of independent conditions versus unknowns, nor for the invertibility of the resulting linear systems for the geometry-dependent coefficients α and β. Since this is the mechanism that defines the vertex basis functions, a proof or a reference to a published proof is needed; the assertion that the construction 'works' is precisely the point that the stress-test identifies as load-bearing.
  3. [§3.2.3, Table 1] The superconvergent point polynomials (15x^4−12x^2+1 and 99x^5−130x^3+31x) are stated without derivation or reference for the regularity cases r=2 and r=3. These points are used for all multi-patch numerical experiments, so their derivation from the 1D model problem (23) should be documented, or a reference provided. This is not a flaw in the method itself, but it makes the convergence claims in Section 4 harder to verify independently.
minor comments (4)
  1. [§2.3.1] In the definition of the interface edge subspace, the index range 'j2 = 5−j1,...,n_{j1}+j1−6' appears to contain a typo; the subscript on n is unexplained and the same expression is garbled in §2.4 as 'n j1 +j1−6'.
  2. [§2.3.2] For boundary vertices of valency ν_i ≥ 3, the text says that g^{(i0)}_{Γ(i0);j1,j2}(ξ2,ξ1) is 'just defined as standard B-splines' N^{p,r}_{j1,j2}(ξ2,ξ1); this notation is confusing because the same symbol is used for interface edge functions in Eq. (8), and the argument ordering should be stated more carefully.
  3. [Example 1] The count of Greville collocation points, |J| = ν n^2 − ν n + 1, is stated without derivation; it would be helpful to explain how repetitions along interfaces and at the central vertex are counted.
  4. [§3.2.3] For (p,r)=(6,2), the text first says 'we have even to add two points' and then states that 'the set of all superconvergent points coincides with the set of clustered superconvergent points'; this phrasing is confusing and should be clarified.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the C2-space construction relies on published prior theorems, and the convergence claims are empirical or derived from an independent 1D model problem.

full rationale

The paper's central derivation is the construction of the C2-smooth space W and its use in isogeometric collocation. The patch/edge subspaces and the C2 interface conditions are taken from the authors' prior works [26,27]; those are published, parameter-free mathematical results, cited as such, and do not incorporate the present paper's numerical outputs, so they constitute independent support rather than circular input. The genuinely new part, the vertex subspaces in Section 2.3.2, is asserted to produce 15 linearly independent C4-at-vertex functions leading to the dimension formula (19), but this is an unproved linear-independence/direct-sum assertion and a potential correctness gap, not a reduction of the conclusion to the assumptions; it is not an instance of a fitted parameter being renamed as a prediction. The superconvergent points in Section 3.2.3 are obtained from an independent 1D model Poisson problem (23) and an asymptotic analysis, not fitted to the multi-patch results, so the convergence rates reported in Section 4 are not forced by construction. No equation in the paper is equivalent by definition to its claimed output, and no cited uniqueness theorem is used to forbid alternatives in a way that would make the choice circular. Therefore the appropriate finding is no significant circularity.

Assumptions & free parameters 1 free parameters · 7 assumptions · 0 invented entities

The central claim rests on several unproved or externally sourced assumptions: the local reparameterization convention, the G2/C2 equivalence, the C2 interface conditions from [26], the h-refinability restrictions, the well-definedness of the vertex functions, the linear independence of the assembled basis, and the transfer of 1D superconvergent points to 2D multi-patch collocation. None of these are demonstrated in depth in this paper; the numerical experiments are the main evidence that the method works.

free parameters (1)
  • Clustered superconvergent point selection rule S_k = For (5,2): depends on k mod 4; for (6,3): omit all inner knots except first and last; for (6,2): add second and…
    Ad hoc choice introduced in Section 3.2.3 to make the number of collocation points equal to the dimension of the underlying spline space; not derived from theory and it affects the claimed convergence behavior.
assumptions (7)
  • domain assumption Local reparameterizations of the geometry mappings as in Fig. 1 are always possible without loss of generality.
    Section 2.1 states this assumption to fix the parameterization of neighboring patches; the C2-smoothness conditions (1)-(3) and the vertex construction (10)-(14) depend on it.
  • standard math An isogeometric function is C2-smooth on a multi-patch domain if and only if its associated graph surface patches are G2-smooth (cf. [16,28]).
    Section 2.2 invokes this equivalence to justify the interface conditions (1)-(3).
  • standard math The conditions (1)-(3) from [26, Lemma 1] characterize C2-smoothness across a single interface.
    Section 2.2 uses these conditions as the definition of the space V; the edge basis functions (7)-(8) are built to satisfy them.
  • domain assumption The assumptions p>=5, 2<=r<=p-3, and k>=(9-p)/(p-r-2) ensure that the constructed C2-smooth spaces are h-refinable and well-defined (see [26,27]).
    Section 2.1 states these restrictions; all numerical examples are limited to (p,r) in {(5,2),(6,2),(6,3)}.
  • ad hoc to paper The vertex functions exist, are uniquely determined by the interpolation conditions (13), and are C2-smooth on the domain and C4-smooth at the vertex.
    Section 2.3.2 asserts this construction is well-defined based on four reasons, but the fourth reason and the analogous boundary-vertex case are not fully proved.
  • ad hoc to paper The functions listed in Section 2.4 are linearly independent and the decomposition (4) is a direct sum, giving the dimension formula (19).
    Section 2.4 states this by construction and by trivial intersection of subspaces, without a detailed proof.
  • domain assumption Superconvergent points computed from the 1D model problem (23) retain their superconvergence for the 2D multi-patch collocation problem.
    Section 3.2.3 generalizes one-patch superconvergent points to multi-patch domains; this is validated only numerically in Section 4, not proven.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Isogeometric collocation on planar multi-patch domains." pith.science (2026). https://pith.science/paper/WACWMRFH

@misc{pith2026190800813,
  author       = {Pith},
  title        = {Pith review of: Isogeometric collocation on planar multi-patch domains},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WACWMRFH}},
  note         = {Machine review of arXiv:1908.00813}
}
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 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.

Figures

Figures reproduced from arXiv: 1908.00813 by the authors.

Figure 1
Figure 1. Left: Considering two neighboring patches Ω [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. As already mentioned in the previous section, the number of (global) collocation [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 2
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. points (without repetitions of same points) is slightly larger than the number of basis func￾tions of the C 2 -smooth space W. But the quotient of the nu… view at source ↗
Figures from the paper (7 more)
Figure 3
Figure 3. Figure 3: All superconvergent points (black and white points) and clustered superconvergent points (only [PITH_FULL_IMAGE:figures/full_fig_p018_3.png]
Figure 4
Figure 4. Figure 4: First row: All superconvergent points (black and white points) and clustered superconvergent [PITH_FULL_IMAGE:figures/full_fig_p018_4.png]
Figure 5
Figure 5. Figure 5: Isogeometric collocation on the bilinearly parameterized unit-square using different sets of col [PITH_FULL_IMAGE:figures/full_fig_p020_5.png]
Figure 6
Figure 6. Figure 6: Computational domains and exact solutions for performing isogeometric collocation in Example 3 [PITH_FULL_IMAGE:figures/full_fig_p022_6.png]
Figure 7
Figure 7. Figure 7: Error plots w.r.t. the number of degrees of freedom (NDOF) of performing isogeometric collo [PITH_FULL_IMAGE:figures/full_fig_p023_7.png]
Figure 8
Figure 8. Figure 8: Error plots w.r.t. the number of degrees of freedom (NDOF) of performing isogeometric collo [PITH_FULL_IMAGE:figures/full_fig_p024_8.png]
Figure 9
Figure 9. Figure 9: Isogeometric collocation on a bicubic three-patch spline geometry using different sets of colloca [PITH_FULL_IMAGE:figures/full_fig_p025_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 35 canonical work pages

  1. [27]

    Kapl and V

    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. 29

  2. [1]

    Anitescu, Y

    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

  3. [2]

    Auricchio, L

    F. Auricchio, L. Beir˜ ao 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

  4. [3]

    Auricchio, L

    F. Auricchio, L. Beir˜ ao 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

  5. [4]

    Beir˜ ao da Veiga, A

    L. Beir˜ ao da Veiga, A. Buffa, G. Sangalli, and R. V´ azquez. Mathematical analysis of variational isogeometric methods. Acta Numerica, 23:157–287, 5 2014

  6. [5]

    Beir˜ ao da Veiga, C

    L. Beir˜ ao 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

  7. [6]

    Bercovier and T

    M. Bercovier and T. Matskewich. Smooth B´ ezier Surfaces over Unstructured Quadri- lateral Meshes. Lecture Notes of the Unione Matematica Italiana, Springer, 2017

  8. [7]

    Blidia, B

    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

Show all 41 references
  1. [8]

    C.L. Chan, C. Anitescu, and T. Rabczuk. Isogeometric analysis with strong multipatch C1-coupling. Comput. Aided Geom. Design , 62:294–310, 2018

  2. [9]

    Collin, G

    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

  3. [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

  4. [11]

    De Lorenzis, J

    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

  5. [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

  6. [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. 28

  7. [14]

    Fahrendorf, L

    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

  8. [15]

    Gomez and L

    H. Gomez and L. De Lorenzis. The variational collocation method. Comput. Methods Appl. Mech. Engrg., 309:152–181, 2016

  9. [16]

    Groisser and J

    D. Groisser and J. Peters. Matched G k-constructions always yield C k-continuous isogeometric elements. Comput. Aided Geom. Des. , 34:67 – 72, 2015

  10. [17]

    Hoschek and D

    J. Hoschek and D. Lasser. Fundamentals of computer aided geometric design . A K Peters Ltd., Wellesley, MA, 1993

  11. [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

  12. [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

  13. [20]

    M. Kapl, F. Buchegger, M. Bercovier, and B. J¨ uttler. Isogeometric analysis with geo- metrically continuous functions on planar multi-patch geometries. Comput. Methods Appl. Mech. Engrg., 316:209 – 234, 2017

  14. [21]

    M. Kapl, G. Sangalli, and T. Takacs. Dimension and basis construction for analysis- suitable G1 two-patch parameterizations. Comput. Aided Geom. Des. , 52–53:75 – 89, 2017

  15. [22]

    M. Kapl, G. Sangalli, and T. Takacs. Isogeometric analysis with C1 functions on unstructured quadrilateral meshes. Technical Report 1812.09088, arXiv.org, 2018

  16. [23]

    M. Kapl, G. Sangalli, and T. Takacs. An isogeometric C1 subspace on unstructured multi-patch planar domains. Comput. Aided Geom. Des. , 69:55–75, 2019

  17. [24]

    Kapl and V

    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

  18. [25]

    Kapl and V

    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

  19. [26]

    Kapl and V

    M. Kapl and V. Vitrih. Dimension and basis construction for C2-smooth isogeometric spline spaces over bilinear-like G2 two-patch parameterizations. J. Comput. Appl. Math., 335:289–311, 2018

  20. [28]

    M. Kapl, V. Vitrih, B. J¨ uttler, and K. Birner. Isogeometric analysis with geometrically continuous functions on two-patch geometries. Comput. Math. Appl. , 70(7):1518 – 1538, 2015

  21. [29]

    Karˇ ciauskas, T

    K. Karˇ ciauskas, T. Nguyen, and J. Peters. Generalizing bicubic splines for modeling and IGA with irregular layout. Comput.-Aided Des., 70:23 – 35, 2016

  22. [30]

    Kiendl, E

    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

  23. [31]

    Marino, J

    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

  24. [32]

    Montardini, G

    M. Montardini, G. Sangalli, and L. Tamellini. Optimal-order isogeometric collocation at Galerkin superconvergent points. Comput. Methods Appl. Mech. Engrg. , 316:741– 757, 2017

  25. [33]

    Mourrain, R

    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

  26. [34]

    Nguyen and J

    T. Nguyen and J. Peters. Refinable C1 spline elements for irregular quad layout. Comput. Aided Geom. Des. , 43:123 – 130, 2016

  27. [35]

    J. Peters. Geometric continuity. In Handbook of computer aided geometric design , pages 193–227. North-Holland, Amsterdam, 2002

  28. [36]

    Reali and T

    A. Reali and T. J. R. Hughes. An introduction to isogeometric collocation methods. In Isogeometric Methods for Numerical Simulation , pages 173–204. Springer, 2015

  29. [37]

    Schillinger, M

    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

  30. [38]

    Schillinger, J

    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

  31. [39]

    Toshniwal, H

    D. Toshniwal, H. Speleers, R. Hiemstra, and T. J. R. Hughes. Multi-degree smooth polar splines: A framework for geometric modeling and isogeometric analysis.Comput. Methods Appl. Mech. Engrg. , 2016

  32. [40]

    Toshniwal, H

    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. 30

  33. [41]

    Toshniwal, H

    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. 31

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.