arxiv: 2602.21685 · v2 · submitted 2026-02-25 · 🧮 math.NA · cs.NA

Recognition: 2 theorem links

· Lean Theorem

Adaptive isogeometric analysis of high-order phase-field fracture based on THB-splines

H.M. Verhelst , L. Greco , A. Reali

Authors on Pith no claims yet

Pith reviewed 2026-05-15 19:44 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords phase-field fractureTHB-splinesadaptive refinementisogeometric analysisAT1 modelAT2 modelhigh-order models2D fracture

0 comments

The pith

THB-splines enable adaptive high-order phase-field fracture simulations in 2D while cutting computational cost.

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

Phase-field models simulate fracture propagation without tracking cracks explicitly, but they demand very fine meshes near the damage zone and therefore carry high computational expense. The paper combines locally refined meshes produced by Truncated Hierarchical B-splines with the higher-order AT1 and AT2 phase-field formulations to perform adaptive simulations. The method is demonstrated on two-dimensional problems and is presented as a practical route to lower the cost of these models. A sympathetic reader would care because the combination makes detailed fracture calculations feasible for engineering use without sacrificing the quality of the predicted crack paths.

Core claim

Leveraging Truncated Hierarchical B-splines (THB-splines), we introduce adaptive simulations of higher-order phase-field formulations (AT1 and AT2), focusing primarily on two-dimensional fracture problems. This approach addresses the high computational cost of phase-field models by pairing local mesh refinement with higher-order approximations in an isogeometric setting.

What carries the argument

Truncated Hierarchical B-splines (THB-splines) for adaptive local refinement in isogeometric analysis paired with high-order AT1 and AT2 phase-field damage models.

If this is right

Local refinement around the fracture zone reduces overall degrees of freedom and runtime for 2D static and dynamic problems.
Fracture path predictions remain comparable to those obtained on uniform meshes.
Both the AT1 and AT2 higher-order models can be treated within the same adaptive THB-spline framework.
The method extends the reach of phase-field analysis to problems that would otherwise be too expensive to resolve at sufficient resolution.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same local-refinement strategy could be applied to three-dimensional fracture geometries once the THB-spline implementation is extended.
Error indicators already available for phase-field models could be used to drive the hierarchical refinement automatically.
Because the underlying geometry representation remains exact, the approach may integrate cleanly with problems that involve curved boundaries or interfaces.

Load-bearing premise

THB-splines combined with high-order phase-field models deliver substantial computational savings while preserving the accuracy of fracture path predictions in the adaptive setting.

What would settle it

Run the adaptive THB-spline code on a standard 2D benchmark fracture problem and compare the resulting crack path and total dissipated energy against a reference solution computed on a uniformly refined mesh of comparable minimum element size; large deviation in either quantity would falsify the accuracy claim.

Figures

Figures reproduced from arXiv: 2602.21685 by A. Reali, H.M. Verhelst, L. Greco.

**Figure 1.** Figure 1: The concept of Truncated Hierarchical B-spline refinement in one dimension. The THB-spline basis in the top row [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗

**Figure 2.** Figure 2: Step-wise illustration of admissible refinement on a mesh corresponding to a THB-spline basis of degree [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗

**Figure 3.** Figure 3: The effect of cross-talk and the remedy of local refinement on a 1D bar. The top row represents a 1D B-spline basis [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

**Figure 4.** Figure 4: Illustration of cross-talk elimination through refinement of support extensions, provided a damage field (in red). The [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗

**Figure 5.** Figure 5: Schematic representation of sudden phase-field propagation combined with mesh adaptivity. From left to right: (a) [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗

**Figure 6.** Figure 6: Benchmark setup for the SEN (b) shear and (a) tensile tests including geometry, boundary conditions and parameters. [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗

**Figure 7.** Figure 7: Contours of the damage field at level d = 0.5 for the (a) tensile and (b) tests. The contours are provided for load steps in the interval u ∈ [0.93, 1.2] · 10−3 [mm] (tensile) and u ∈ [0.45, 0.60] · 10−3 [mm] (shear) with step size ∆u = 0.03 · 10−3 [mm] as described in figures 6a and 6b. The rows represent the results obtained on a coarse (h = ℓ0/2) and fine (h = ℓ0/4) mesh, while the columns represent dif… view at source ↗

**Figure 8.** Figure 8: Reaction force Fx, Fy (left) and dissipated energy D (right) for the (a) tensile and (b) shear tests with respect to the displacement at the top boundary, comparing different phase-field models on different meshes. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗

**Figure 9.** Figure 9: Contour plots of the damage field d at d = 0.5 at different horizontal or vertical displacements u for the tensile (a) and shear (b) tests, comparing the tensor-product B-spline basis with the adaptive implicit, explicit and hybrid approaches using the fourth-order AT1 model on a mesh with coarsest mesh size h = ℓ0 2 . 20 [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗

**Figure 10.** Figure 10: Dissipation D (top row) and total CPU-time (bottom row) comparison for the (a) tensile and (b) shear tests using the second order AT1 (left), fourth order AT1 (middle left), second order AT2 (middle right) and fourth order AT2 (right) models, all plotted against the vertical displacement of the top boundary uy. All results are provided for the fine mesh (h = ℓ0/4) except for the fourth order AT-1 model, a… view at source ↗

**Figure 11.** Figure 11: Computational cost comparison for different phase-field models using adaptive meshing for the (a) tensile and (b) [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗

**Figure 12.** Figure 12: Comparison of load stepping approaches for adaptive THB-spline phase-field models in terms of computational cost [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗

read the original abstract

In recent decades, the study of fracture propagation in solids has increasingly relied on phase-field models. Several recent contributions have highlighted the potential of this approach in both static and dynamic frameworks. However, a major limitation remains the high computational cost. Two main strategies have been identified to mitigate this issue: the use of locally refined meshes and the adoption of higher-order models. In this work, leveraging Truncated Hierarchical B-splines (THB-splines), we introduce adaptive simulations of higher-order phase-field formulations (AT1 and AT2), focusing primarily on two-dimensional fracture problems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper combines THB-splines with adaptive high-order AT1/AT2 phase-field fracture in 2D to cut computational cost, a logical extension that holds together on paper but needs concrete benchmarks to confirm the gains.

read the letter

The main takeaway is that this work puts together truncated hierarchical B-splines for local refinement with higher-order phase-field models for 2D fracture. The approach targets the usual expense of phase-field simulations by allowing adaptive meshes while keeping the smoothness needed for the AT1 and AT2 formulations. THB-splines are a known tool that preserves partition of unity and approximation order under truncation, so the construction fits the fourth-order or second-order equations without obvious contradictions. The 2D focus also sidesteps extra stabilization issues that would appear in 3D. This is a practical step rather than a new theory, building on prior isogeometric and phase-field papers in a direct way. The method description in the abstract is consistent, and the stress-test note correctly flags that the truncation operator should maintain the required orders for residual or gradient-based adaptivity. The soft spot is that the abstract gives no numerical results, error tables, or timing comparisons, so the claim of substantial savings while preserving fracture path accuracy remains unverified from the text alone. If the full manuscript includes those benchmarks against uniform meshes or lower-order runs, the efficiency argument would land more solidly. Minor gaps like missing implementation details on marking strategies or solver choices could be filled in revision. This is written for people already working in isogeometric analysis or computational fracture mechanics who care about adaptive strategies for diffuse models. It shows clear engagement with the literature and no internal inconsistencies, so it deserves a serious referee. I would send it to peer review with the expectation that the authors add direct evidence of the cost reductions.

Referee Report

0 major / 3 minor

Summary. The manuscript introduces an adaptive isogeometric analysis framework for high-order phase-field fracture models (AT1 and AT2) in two dimensions, employing Truncated Hierarchical B-splines (THB-splines) to enable local mesh refinement while preserving partition of unity and smoothness. The central claim is that this combination mitigates the high computational cost of phase-field simulations by leveraging adaptive refinement driven by residual or gradient indicators on the phase-field variable, with the truncation operator maintaining the required approximation order for the underlying fourth- or second-order equations.

Significance. If the numerical experiments confirm substantial efficiency gains (e.g., reduced degrees of freedom and CPU time) while accurately capturing fracture paths in standard 2D benchmarks, the work would represent a useful extension of isogeometric analysis to adaptive high-order phase-field fracture. The approach builds directly on established properties of THB-splines and diffuse-interface regularization without introducing internal inconsistencies in the construction.

minor comments (3)

The abstract states that the method focuses on two-dimensional fracture problems but does not indicate which specific benchmarks (e.g., single-edge notched tension or shear tests) are used to demonstrate the adaptive refinement and efficiency claims.
In the numerical results section, ensure that tables or figures explicitly report both the number of degrees of freedom and wall-clock times for adaptive versus uniform meshes at comparable accuracy levels, so that the claimed computational savings can be directly verified.
Clarify the precise form of the adaptive marking criterion (residual-based or gradient-based) and its dependence on the phase-field regularization length scale, as this choice directly affects the observed refinement patterns.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript on adaptive isogeometric analysis of high-order phase-field fracture using THB-splines. We appreciate the recognition of the approach's potential for efficiency gains in 2D benchmarks and the recommendation for minor revision. No specific major comments were raised in the report, so we will focus on minor improvements to presentation and clarity in the revised version.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper presents an adaptive IGA framework that combines established THB-spline refinement properties with standard high-order AT1/AT2 phase-field models. No derivation chain, equation, or central claim reduces by construction to a fitted parameter, self-definition, or unverified self-citation. The truncation operator, partition-of-unity preservation, and adaptive marking are invoked from prior independent literature on THB-splines and phase-field regularization; the 2-D numerical examples serve as verification rather than tautological prediction. The construction therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the approach rests on established concepts from isogeometric analysis and phase-field modeling.

pith-pipeline@v0.9.0 · 5396 in / 971 out tokens · 19382 ms · 2026-05-15T19:44:09.734893+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

leveraging Truncated Hierarchical B-splines (THB-splines), we introduce adaptive simulations of higher-order phase-field formulations (AT1 and AT2)
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

admissible refinement algorithms for THB-splines based on the T-neighborhood

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

[1]

A. A. Griffith, G. I. Taylor, Vi. the phenomena of rupture and flow in solids, Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character 221 (1921) 163–198

work page 1921
[2]

G. R. Irwin, Analysis of stresses and strains near the end of a crack traversing a plate (1957)

work page 1957
[3]

R. P. Reed, The economic effects of fracture in the United States, volume 647, US Department of Commerce, National Bureau of Standards, 1983

work page 1983
[4]

Belytschko, T

T. Belytschko, T. Black, Elastic crack growth in finite elements with minimal remeshing, International journal for numerical methods in engineering 45 (1999) 601–620

work page 1999
[5]

Bourdin, G

B. Bourdin, G. Francfort, J.-J. Marigo, Numerical experiments in revisited brittle fracture, Journal of the Mechanics and Physics of Solids 48 (2000) 797–826

work page 2000
[6]

Borden, T

M. Borden, T. Hughes, C. Landis, C. Verhoosel, A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework, Computer Methods in Applied Mechanics and Engineering 273 (2014) 100–118

work page 2014
[7]

Greco, A

L. Greco, A. Patton, M. Negri, A. Marengo, U. Perego, A. Reali, Higher order phase-field modeling of brittle fracture via isogeometric analysis, Engineering with Computers (2024)

work page 2024
[8]

Ambrosio, V

L. Ambrosio, V. M. Tortorelli, Approximation of functional depending on jumps by elliptic functional via t-convergence, Communications on Pure and Applied Mathematics 43 (1990) 999 – 1036

work page 1990
[9]

Hughes, J

T. Hughes, J. Cottrell, Y. Bazilevs, Isogeometric analysis: Cad, finite elements, nurbs, exact geometry and mesh refine- ment, Computer Methods in Applied Mechanics and Engineering 194 (2005) 4135–4195

work page 2005
[10]

Cottrell, T

J. Cottrell, T. Hughes, A. Reali, Studies of refinement and continuity in isogeometric structural analysis, Computer Methods in Applied Mechanics and Engineering 196 (2007) 4160–4183

work page 2007
[11]

Miehe, F

C. Miehe, F. Welschinger, M. Hofacker, Thermodynamically consistent phase-field models of fracture: Variational princi- ples and multi-field fe implementations, International Journal for Numerical Methods in Engineering 83 (2010) 1273–1311

work page 2010
[12]

Gerasimov, L

T. Gerasimov, L. De Lorenzis, On penalization in variational phase-field models of brittle fracture, Computer Methods in Applied Mechanics and Engineering 354 (2019) 990–1026

work page 2019
[13]

Marengo, A

A. Marengo, A. Patton, M. Negri, U. Perego, A. Reali, A rigorous and efficient explicit algorithm for irreversibility enforcement in phase-field finite element modeling of brittle crack propagation, Computer Methods in Applied Mechanics and Engineering 387 (2021) 114137

work page 2021
[14]

Proserpio, M

D. Proserpio, M. Ambati, L. De Lorenzis, J. Kiendl, A framework for efficient isogeometric computations of phase-field brittle fracture in multipatch shell structures, Computer Methods in Applied Mechanics and Engineering 372 (2020) 113363

work page 2020
[15]

Y. Li, T. Yu, S. Natarajan, T. Q. Bui, A dynamic description of material brittle failure using a hybrid phase-field model enhanced by adaptive isogeometric analysis, European Journal of Mechanics-A/Solids 97 (2023) 104783

work page 2023
[16]

Y. Li, T. Yu, C. Xing, S. Natarajan, Modeling quasi-static and dynamic thermo-elastic coupled brittle fracture using an adaptive isogeometric hybrid phase-field method, Finite Elements in Analysis and Design 224 (2023) 103993

work page 2023
[17]

Y. Li, T. Yu, S. Natarajan, An adaptive isogeometric phase-field method for brittle fracture in rock-like materials, Engineering Fracture Mechanics 263 (2022) 108298

work page 2022
[18]

Kiran, N

R. Kiran, N. Nguyen-Thanh, K. Zhou, Adaptive isogeometric analysis–based phase-field modeling of brittle electrome- chanical fracture in piezoceramics, Engineering Fracture Mechanics 274 (2022) 108738

work page 2022
[19]

Kiran, N

R. Kiran, N. Nguyen-Thanh, H. Yu, K. Zhou, Adaptive isogeometric analysis–based phase-field modeling of interfacial fracture in piezoelectric composites, Engineering Fracture Mechanics 288 (2023) 109181

work page 2023
[20]

Kiran, N

R. Kiran, N. Nguyen-Thanh, K. Zhou, On the phase-field modeling of thermo-electromechanical brittle fracture in piezoceramics using an adaptive isogeometric approach, Engineering Fracture Mechanics 304 (2024) 110068

work page 2024
[21]

B. Xu, T. Xu, M. J. Heap, A. R. Kushnir, B.-y. Su, X.-c. Lan, An adaptive phase field approach to 3d internal crack growth in rocks, Computers and Geotechnics 173 (2024) 106551

work page 2024
[22]

Nguyen-Thanh, W

N. Nguyen-Thanh, W. Li, J. Huang, K. Zhou, Adaptive higher-order phase-field modeling of anisotropic brittle fracture in 3d polycrystalline materials, Computer Methods in Applied Mechanics and Engineering 372 (2020) 113434

work page 2020
[23]

W. Li, N. Nguyen-Thanh, K. Zhou, Phase-field modeling of brittle fracture in a 3d polycrystalline material via an adaptive isogeometric-meshfree approach, International Journal for Numerical Methods in Engineering 121 (2020) 5042–5065

work page 2020
[24]

Giannelli, B

C. Giannelli, B. J¨ uttler, H. Speleers, Thb-splines: The truncated basis for hierarchical splines, Computer Aided Geometric Design 29 (2012) 485–498

work page 2012
[25]

Lorenzo, M

G. Lorenzo, M. Scott, K. Tew, T. Hughes, H. Gomez, Hierarchically refined and coarsened splines for moving interface problems, with particular application to phase-field models of prostate tumor growth, Computer Methods in Applied Mechanics and Engineering 319 (2017) 515–548. 25

work page 2017
[26]

Bracco, C

C. Bracco, C. Giannelli, A. Reali, M. Torre, R. V´ azquez, Adaptive isogeometric phase-field modeling of the cahn–hilliard equation: Suitably graded hierarchical refinement and coarsening on multi-patch geometries, Computer Methods in Applied Mechanics and Engineering 417 (2023) 116355

work page 2023
[27]

Greco, E

L. Greco, E. Maggiorelli, M. Negri, A. Patton, A. Reali, At1 fourth-order isogeometric phase-field modeling of brittle fracture, Mathematical Models and Methods in Applied Sciences 0 (0) 1–55

work page
[28]

C. Comi, U. Perego, Fracture energy based bi-dissipative damage model for concrete, International Journal of Solids and Structures 38 (2001) 6427–6454

work page 2001
[29]

Amor, J.-J

H. Amor, J.-J. Marigo, C. Maurini, Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments, Journal of the Mechanics and Physics of Solids 57 (2009) 1209–1229

work page 2009
[30]

Negri,γ-convergence for high order phase field fracture: Continuum and isogeometric formulations, Computer Methods in Applied Mechanics and Engineering 362 (2020) 112858

M. Negri,γ-convergence for high order phase field fracture: Continuum and isogeometric formulations, Computer Methods in Applied Mechanics and Engineering 362 (2020) 112858

work page 2020
[31]

Piegl, W

L. Piegl, W. Tiller, The NURBS book, Springer Science and Business Media, 1996

work page 1996
[32]

Borden, T

M. Borden, T. Hughes, C. Landis, C. Verhoosel, A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework, Comput. Methods Appl. Mech. Engrg. 273 (2014) 100 – 118

work page 2014
[33]

Miehe, M

C. Miehe, M. Hofacker, F. Welschinger, A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits, Computer Methods in Applied Mechanics and Engineering 199 (2010) 2765–2778

work page 2010
[34]

J. M. Sargado, E. Keilegavlen, I. Berre, J. M. Nordbotten, High-accuracy phase-field models for brittle fracture based on a new family of degradation functions, Journal of the Mechanics and Physics of Solids 111 (2018) 458–489

work page 2018
[35]

H. M. Verhelst, P. Weinm¨ uller, A. Mantzaflaris, T. Takacs, D. Toshniwal, A comparison of smooth basis constructions for isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 419 (2024) 116659

work page 2024
[36]

de Prenter, C

F. de Prenter, C. V. Verhoosel, E. H. van Brummelen, M. G. Larson, S. Badia, Stability and Conditioning of Immersed Finite Element Methods: Analysis and Remedies, Archives of Computational Methods in Engineering (2023) 1–61

work page 2023
[37]

Vuong, C

A.-V. Vuong, C. Giannelli, B. J¨ uttler, B. Simeon, A hierarchical approach to adaptive local refinement in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 200 (2011) 3554–3567

work page 2011
[38]

Giannelli, B

C. Giannelli, B. J¨ uttler, H. Speleers, THB-splines: The truncated basis for hierarchical splines, Computer Aided Geometric Design 29 (2012) 485–498

work page 2012
[39]

Giannelli, B

C. Giannelli, B. J¨ uttler, S. K. Kleiss, A. Mantzaflaris, B. Simeon, J. ˇSpeh, THB-splines: An effective mathematical technology for adaptive refinement in geometric design and isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 299 (2016) 337–365

work page 2016
[40]

Buffa, C

A. Buffa, C. Giannelli, Adaptive isogeometric methods with hierarchical splines: Error estimator and convergence, Mathematical Models and Methods in Applied Sciences 26 (2016) 1–25

work page 2016
[41]

Bracco, C

C. Bracco, C. Giannelli, R. V´ azquez, Refinement Algorithms for Adaptive Isogeometric Methods with Hierarchical Splines, Axioms 7 (2018) 43

work page 2018
[42]

Carraturo, C

M. Carraturo, C. Giannelli, A. Reali, R. V´ azquez, Suitably graded THB-spline refinement and coarsening: Towards an adaptive isogeometric analysis of additive manufacturing processes, Computer Methods in Applied Mechanics and Engineering 348 (2019) 660–679

work page 2019
[43]

Gantner, D

G. Gantner, D. Haberlik, D. Praetorius, Adaptive IGAFEM with optimal convergence rates: Hierarchical B-splines, Mathematical Models and Methods in Applied Sciences 27 (2017) 2631–2674

work page 2017
[44]

Buffa, G

A. Buffa, G. Gantner, C. Giannelli, D. Praetorius, R. V´ azquez, Mathematical Foundations of Adaptive Isogeometric Analysis, Archives of Computational Methods in Engineering 29 (2022) 4479–4555

work page 2022
[45]

Z. Lian, L. Leidinger, S. Hartmann, F. Bauer, M. Pabst, C. Krisadawat, R. W¨ uchner, Cross-talk effects in trimmed isogeometric shells and the control point duplication approach, Computer Methods in Applied Mechanics and Engineering 438 (2025) 117849

work page 2025
[46]

Mantzaflaris, H

A. Mantzaflaris, H. M. Verhelst, M. M¨ oller, C. Karampatzakis, S. Takacs, S. Imperatore, P. Weinm¨ uller, J. Vogl, J. Li, Y. Ji, M. Matucci, A. Limkilde, F. Scholz, D. Mokriˇ s, J. Sogn, A. Shamanskiy, E. Vollebregt, R. Tielen, A. Bressan, A. Farahat, Geometry + Simulation Modules (G+Smo), Zenodo, 2025. doi:10.5281/zenodo.15874172

work page doi:10.5281/zenodo.15874172 2025
[47]

J¨ uttler, U

B. J¨ uttler, U. Langer, A. Mantzaflaris, S. E. Moore, W. Zulehner, Geometry + Simulation Modules: Implementing Isogeometric Analysis, PAMM 14 (2014) 961–962

work page 2014
[48]

Schenk, K

O. Schenk, K. G¨ artner, Solving unsymmetric sparse systems of linear equations with pardiso, Future Generation Computer Systems 20 (2004) 475–487. Appendix A. Algorithmic overview In this section, we provide an overview of the algorithms used to solve the coupled phase-field fracture problem with adaptive refinement. In algorithms 1 and 2 we provide subr...

work page 2004