pith. sign in

arxiv: 2604.09113 · v1 · submitted 2026-04-10 · 🧮 math.NA · cs.NA

A ROM-based BDDC solver for unfitted p-FEM level-set-based lattice structures

Pith reviewed 2026-05-10 17:06 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords BDDCreduced order modelunfitted finite elementslevel setlattice structuresdomain decompositionp-FEMMDEIM
0
0 comments X

The pith

A BDDC solver accelerated by reduced-order models computes large level-set lattice structures with bounded iterations on standard hardware.

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

The paper develops a domain-decomposition approach for simulating lattice structures whose cells are defined by level-set functions and arbitrary mappings. Each cell is treated as a single high-order unfitted finite element, and the Balanced Domain Decomposition by Constraints method assigns one subdomain per cell. Assembly of the cell stiffness matrices is accelerated by a reduced-order model based on matrix discrete empirical interpolation that is trained once offline and reused across different geometries. A stabilization term is added to preserve the solver's theoretical scalability properties despite the approximation. Tests on a graded 2D lattice exceeding 17,000 cells confirm that iteration counts remain bounded when the subdomain-to-mesh-size ratio is held fixed and that the full solve completes in roughly 30 seconds on a laptop.

Core claim

The central claim is that the combination of unfitted p-FEM discretization, BDDC domain decomposition with one subdomain per cell, and an MDEIM-based reduced-order model for stiffness-matrix assembly, together with a stabilization term, produces an accurate and scalable solver for arbitrarily graded and topologically varying lattice structures described by level sets, without invoking homogenization or periodicity assumptions.

What carries the argument

The MDEIM reduced-order surrogate for cell stiffness matrices, trained offline on geometric-mapping contributions and combined with a stabilization term inside the BDDC preconditioner for unfitted elements.

If this is right

  • Large graded lattices with thousands of distinct cell geometries can be solved directly at full resolution rather than through homogenization.
  • Solver iterations stay bounded when the ratio of subdomain size to mesh size is fixed, matching classical BDDC theory.
  • Expensive quadrature on cut elements is confined to an offline training stage and does not recur during repeated solves or design iterations.
  • The method accommodates arbitrary mappings and topology changes between cells without requiring scale separation.

Where Pith is reading between the lines

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

  • The same offline-training strategy could be applied to three-dimensional lattices once the ROM basis is extended to higher-dimensional mappings.
  • Because the reduced-order model is independent of the specific physical problem, the framework might be reused across different constitutive laws inside the same lattice geometry.
  • The controllable error introduced by stabilization opens a route to trading a modest accuracy loss for further speed gains in early-stage design exploration.

Load-bearing premise

The stabilization term added to the reduced-order approximation keeps the introduced error small enough that it does not degrade solution quality or destroy BDDC scalability in lattices with varying cell topologies.

What would settle it

If the number of BDDC iterations grows without bound as the number of subdomains increases while the ratio of subdomain diameter to local mesh size is kept constant, the claimed scalability would be refuted.

Figures

Figures reproduced from arXiv: 2604.09113 by Giuliano Guarino, Gonzalo Bonilla Moreno, Pablo Antolin.

Figure 1
Figure 1. Figure 1: Construction of a single cell of the lattice structure. The untrimmed parametric domain [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: 2D TPMS corresponding to the level-set functions in Table [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: A lattice structure is constructed by juxtaposing six cells in a compatible way. The compatibility [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Position of the Gauss-Lobatto-Legendre nodes of the Lagrangian basis superimposed to the trimmed [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Classification of the DoFs associated with a single cell (a) and with a two-by-two cells structure (b). [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Classification of the coarse DoFs associated with a single cell (a) and with a two-by-two cells [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison between unfitted p-FEM element method and CutFEM for the manufactured problem [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Benchmark geometry for the tests in Sections [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (a) Accuracy of the fast assembly interpolation measured in the [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Error of the fast assembly tensor depending on the basis size [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Effect of the stabilization in the solver iterations and the solution error. (a) Iterations with and [PITH_FULL_IMAGE:figures/full_fig_p023_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Performances of three different solvers are compared: BDDC, Cholesky decomposition, and SOR. [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Iterations (a) and computational time (b) associated with the solution of the problem using a [PITH_FULL_IMAGE:figures/full_fig_p025_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Iterations (a) and computational time (b) as a function of number of cells. In (b), the computa [PITH_FULL_IMAGE:figures/full_fig_p026_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: (a) Geometry of the sandwich wing example in Section [PITH_FULL_IMAGE:figures/full_fig_p026_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: (a) Geometry of lattice structure representing a wrench as described in Section [PITH_FULL_IMAGE:figures/full_fig_p027_16.png] view at source ↗
read the original abstract

We present a domain decomposition method for the fast simulation of large lattice structures described by level set functions. The method does not rely on homogenization or multiscale techniques, and therefore avoids their underlying assumptions such as scale separation and periodicity. Individual cells are defined through level set functions and mapped into physical space using arbitrary order mappings, allowing the creation of complex graded designs with varying geometries and topologies. The discretization is based on unfitted p-FEM, where each cell is approximated by a single high order element. This choice naturally handles the implicit geometric description and provides high accuracy with a moderate number of degrees of freedom. The solver is built on the Balanced Domain Decomposition by Constraints (BDDC) method, where each cell corresponds to one subdomain. To accelerate the assembly of the cell stiffness matrices, we combine a fast assembly technique that separates the contributions of the geometric mapping from the trimmed domain with a reduced order model (ROM) based on the matrix discrete empirical interpolation method (MDEIM). The ROM surrogate is trained offline and reused for any geometric mapping, restricting the expensive quadrature on cut elements to the training stage. A stabilization term ensures the scalability of the solver when using the ROM approximation, at the cost of a small and controllable error. We validate the method through numerical experiments and demonstrate its performance on a complex 2D problem with more than 17,000 cells of varying geometry, solved in approximately 30 seconds on a standard laptop. The number of solver iterations remains bounded as the number of subdomains grows, provided the ratio between subdomain and mesh sizes is kept constant, in agreement with classical BDDC scalability properties.

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.

Referee Report

3 major / 3 minor

Summary. The paper proposes a BDDC domain decomposition solver for unfitted p-FEM discretizations of level-set-based lattice structures. Each cell is treated as a subdomain and discretized with a single high-order element; cell stiffness matrices are assembled via a combination of geometric-mapping separation and an MDEIM-based reduced-order model trained offline. A stabilization term is added to restore the stability properties required for BDDC scalability when the ROM surrogate is used. The method is demonstrated on a single large 2D graded lattice containing more than 17,000 cells of varying geometry and topology, solved in roughly 30 seconds on a laptop, with iteration counts remaining bounded when the subdomain-to-mesh-size ratio is held constant.

Significance. If the stabilization term can be shown to preserve both BDDC iteration bounds and solution accuracy across topologically varying cells, the approach would provide a practical route to direct high-order simulation of complex, non-periodic lattices without homogenization assumptions. The offline ROM training and reuse for arbitrary mappings addresses a key assembly bottleneck in unfitted p-FEM, and the reported 30-second solve time on 17k cells indicates engineering relevance. However, the current evidence rests on one numerical example and the assumption that the ROM-induced perturbation remains controllable.

major comments (3)
  1. [Numerical experiments] Numerical experiments section: the single 17,000-cell demonstration reports bounded iterations and acceptable runtime, yet contains no systematic convergence study, no comparison against full quadrature assembly, and no quantification of the L2 or energy error introduced by the ROM-plus-stabilization approximation across cells with differing cut configurations. This is load-bearing for the central claim that the error remains “small and controllable.”
  2. [Stabilization term] Stabilization term description (likely §3): the manuscript states that the added stabilization restores scalability “at the cost of a small and controllable error,” but provides neither an a-priori bound on the perturbation to the local stiffness matrices nor an analysis of how the stabilization coefficient interacts with the varying Jacobians and cut topologies present in graded lattices. Without such analysis the extension of classical BDDC theory to the ROM setting remains empirical.
  3. [ROM training] ROM training and reuse (Section on MDEIM): the offline training set is described as sufficient for “any geometric mapping,” yet no coverage study or sensitivity test is reported for the range of mappings and topology changes encountered in the target graded structures. If the training set under-samples certain cut configurations, the local matrix error can propagate through the BDDC coarse space, undermining the observed iteration bound.
minor comments (3)
  1. [Abstract] The abstract claims agreement with “classical BDDC scalability properties,” but the manuscript does not restate the precise assumptions (e.g., quasi-uniform subdomain partitioning, sufficient overlap) under which the ROM perturbation is expected to preserve those properties.
  2. [Numerical experiments] Hardware and implementation details for the 30-second timing are not provided; a brief statement of processor, memory, and whether the code is serial or parallel would allow readers to assess reproducibility.
  3. [Notation] Notation for the stabilization coefficient and the MDEIM interpolation points should be introduced once and used consistently; occasional re-definition of symbols across sections reduces readability.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive and detailed comments. We address each major comment point by point below, indicating the revisions we will make to strengthen the manuscript while remaining faithful to the existing results.

read point-by-point responses
  1. Referee: [Numerical experiments] Numerical experiments section: the single 17,000-cell demonstration reports bounded iterations and acceptable runtime, yet contains no systematic convergence study, no comparison against full quadrature assembly, and no quantification of the L2 or energy error introduced by the ROM-plus-stabilization approximation across cells with differing cut configurations. This is load-bearing for the central claim that the error remains “small and controllable.”

    Authors: We agree that additional numerical evidence would better support the claim of controllable error. In the revised manuscript we will add: (i) a systematic convergence study for representative single cells with varying cut configurations and polynomial degrees, comparing ROM-based solutions against full-quadrature references in both L2 and energy norms; (ii) a direct accuracy comparison between ROM and full assembly on a smaller multi-cell lattice; and (iii) tabulated relative matrix and solution errors across a range of cut topologies. These additions will quantify the approximation error and its propagation. revision: yes

  2. Referee: [Stabilization term] Stabilization term description (likely §3): the manuscript states that the added stabilization restores scalability “at the cost of a small and controllable error,” but provides neither an a-priori bound on the perturbation to the local stiffness matrices nor an analysis of how the stabilization coefficient interacts with the varying Jacobians and cut topologies present in graded lattices. Without such analysis the extension of classical BDDC theory to the ROM setting remains empirical.

    Authors: The stabilization term is introduced to recover the spectral equivalence properties needed for BDDC scalability when the local matrices are replaced by their ROM surrogates. While the current manuscript relies on numerical evidence that iteration counts remain bounded and comparable to the non-ROM case, we acknowledge the absence of an a-priori perturbation bound. In the revision we will expand the description of the stabilization coefficient, report its sensitivity to Jacobian variation and cut topology through additional numerical tests, and clarify that the theoretical extension remains empirical at present. revision: partial

  3. Referee: [ROM training] ROM training and reuse (Section on MDEIM): the offline training set is described as sufficient for “any geometric mapping,” yet no coverage study or sensitivity test is reported for the range of mappings and topology changes encountered in the target graded structures. If the training set under-samples certain cut configurations, the local matrix error can propagate through the BDDC coarse space, undermining the observed iteration bound.

    Authors: The MDEIM training set was constructed from a diverse collection of mappings and cut configurations representative of graded lattices. We agree that explicit coverage and sensitivity results would increase confidence in generalization. The revised manuscript will include a new subsection detailing the training-set composition together with approximation-error metrics for both in-distribution and out-of-distribution mappings and topologies, including extreme cut ratios, to assess potential impact on the coarse-space correction. revision: yes

standing simulated objections not resolved
  • A rigorous a-priori bound on the ROM-induced perturbation to the local stiffness matrices and its interaction with arbitrary Jacobians and cut topologies.

Circularity Check

0 steps flagged

No significant circularity; derivation builds on external BDDC and MDEIM theory with independent numerical validation

full rationale

The paper's core chain—unfitted p-FEM discretization of level-set lattices, BDDC domain decomposition with one cell per subdomain, MDEIM-based ROM for fast stiffness assembly, and an added stabilization term—relies on established external references for BDDC scalability and MDEIM training. The stabilization is introduced to restore properties lost under ROM approximation and is validated numerically on a 17k-cell graded lattice example, with iteration counts shown to remain bounded under standard subdomain-to-mesh ratio conditions. No step reduces a claimed prediction or uniqueness result to a self-fit, self-citation chain, or definitional renaming; the error controllability is treated as an empirical assumption supported by experiments rather than derived tautologically from the target data.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the ROM surrogate preserving sufficient accuracy for BDDC scalability when paired with stabilization; the offline training step and the constant subdomain-to-mesh-size ratio are key external requirements.

free parameters (1)
  • stabilization coefficient
    Introduced to restore scalability under ROM approximation; its specific value is chosen for the reported experiments but not quantified in the abstract.
axioms (2)
  • domain assumption Classical BDDC iteration bounds hold when the ratio of subdomain size to mesh size is held constant
    Invoked to explain the observed bounded iteration count.
  • ad hoc to paper The ROM approximation error remains small and controllable after stabilization
    Added to justify use of the surrogate in the solver.

pith-pipeline@v0.9.0 · 5601 in / 1424 out tokens · 83701 ms · 2026-05-10T17:06:42.826849+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

74 extracted references · 74 canonical work pages

  1. [1]

    T. D. Ngo, A. Kashani, G. Imbalzano, K. T. Nguyen, D. Hui, Additive Manufacturing (3D Printing): A Review of Materials, Methods, Applications and Challenges, Composites Part B: Engineering 143 (2018) 172–196

  2. [2]

    Vyatskikh, S

    A. Vyatskikh, S. Delalande, A. Kudo, X. Zhang, C. M. Portela, J. R. Greer, Additive Manufacturing of 3D Nano-Architected Metals, Nature communications 9 (1) (2018) 1–8

  3. [3]

    C. Pan, Y. Han, J. Lu, Design and Optimization of Lattice Structures: A Review, Applied Sciences 10 (18) (2020) 6374

  4. [4]

    Zheng, H

    X. Zheng, H. Lee, T. H. Weisgraber, M. Shusteff, J. DeOtte, E. B. Duoss, J. D. Kuntz, M. M. Biener, Q. Ge, J. A. Jackson, et al., Ultralight, Ultrastiff Mechanical Metamaterials, Science 344 (6190) (2014) 1373–1377

  5. [5]

    Bauer, S

    J. Bauer, S. Hengsbach, I. Tesari, R. Schwaiger, O. Kraft, High-Strength Cellular Ceramic Composites with 3D Microarchitecture, Proceedings of the National Academy of Sciences 111 (7) (2014) 2453–2458

  6. [6]

    S. C. Han, J. W. Lee, K. Kang, A New Type of Low Density Material: Shellular, Advanced Materials 27 (37) (2015) 5506–5511

  7. [7]

    L. R. Meza, A. J. Zelhofer, N. Clarke, A. J. Mateos, D. M. Kochmann, J. R. Greer, Resilient 3D Hierarchical Architected Metamaterials, Proceedings of the National Academy of Sciences 112 (37) (2015) 11502–11507

  8. [8]

    A. J. D. Shaikeea, H. Cui, M. O’Masta, X. R. Zheng, V. S. Deshpande, The Toughness of Mechanical Metamaterials, Nature Materials 21 (3) (2022) 297–304

  9. [9]

    X. Ren, R. Das, P. Tran, T. D. Ngo, Y. M. Xie, Auxetic Metamaterials and Structures: A Review, Smart materials and structures 27 (2) (2018) 023001

  10. [10]

    Tancogne-Dejean, A

    T. Tancogne-Dejean, A. B. Spierings, D. Mohr, Additively-Manufactured Metallic Micro-Lattice Materi- als for High Specific Energy Absorption under Static and Dynamic Loading, Acta Materialia 116 (2016) 14–28

  11. [11]

    S. Shan, S. H. Kang, J. R. Raney, P. Wang, L. Fang, F. Candido, J. A. Lewis, K. Bertoldi, Multistable Architected Materials for Trapping Elastic Strain Energy, Advanced Materials 27 (29) (2015) 4296–4301

  12. [12]

    J. U. Surjadi, L. Gao, H. Du, X. Li, X. Xiong, N. X. Fang, Y. Lu, Mechanical Metamaterials and Their Engineering Applications, Advanced Engineering Materials 21 (3) (2019) 1800864

  13. [13]

    D. M. Kochmann, J. B. Hopkins, L. Valdevit, Multiscale Modeling and Optimization of the Mechanics of Hierarchical Metamaterials, MRS Bulletin 44 (10) (2019) 773–781

  14. [14]

    M. Buck, O. Iliev, H. Andr¨ a, Multiscale Finite Element Coarse Spaces for the Application to Linear Elasticity, Open Mathematics 11 (4) (2013) 680–701

  15. [15]

    T. Y. Hou, X.-H. Wu, A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media, Journal of computational physics 134 (1) (1997) 169–189

  16. [16]

    Castelletto, H

    N. Castelletto, H. Hajibeygi, H. A. Tchelepi, Multiscale Finite-Element Method for Linear Elastic Ge- omechanics, Journal of Computational Physics 331 (2017) 337–356

  17. [17]

    Doˇ sk´ aˇ r, J

    M. Doˇ sk´ aˇ r, J. Zeman, P. Krysl, J. Nov´ ak, Microstructure-Informed Reduced Modes Synthesized with Wang Tiles and the Generalized Finite Element Method, Computational Mechanics 68 (2) (2021) 233– 253

  18. [18]

    Savvas, G

    D. Savvas, G. Stefanou, M. Papadrakakis, G. Deodatis, Homogenization of Random Heterogeneous Media with Inclusions of Arbitrary Shape Modeled by XFEM, Computational mechanics 54 (5) (2014) 1221–1235. 30

  19. [19]

    Strouboulis, I

    T. Strouboulis, I. Babuˇ ska, K. Copps, The Design and Analysis of the Generalized Finite Element Method, Computer methods in applied mechanics and engineering 181 (1-3) (2000) 43–69

  20. [20]

    F. Feyel, A Multilevel Finite Element Method (FE2) to Describe the Response of Highly Non-Linear Structures Using Generalized Continua, Computer Methods in applied Mechanics and engineering 192 (28-30) (2003) 3233–3244

  21. [21]

    Schr¨ oder, A Numerical Two-Scale Homogenization Scheme: the FE 2-Method, in: Plasticity and beyond, Springer, 2014, pp

    J. Schr¨ oder, A Numerical Two-Scale Homogenization Scheme: the FE 2-Method, in: Plasticity and beyond, Springer, 2014, pp. 1–64

  22. [22]

    Charalambakis, Homogenization Techniques and Micromechanics: A Survey and Perspectives, Ap- plied Mechanics Reviews 63 (3) (2010)

    N. Charalambakis, Homogenization Techniques and Micromechanics: A Survey and Perspectives, Ap- plied Mechanics Reviews 63 (3) (2010)

  23. [23]

    Moulinec, P

    H. Moulinec, P. Suquet, A FFT-Based Numerical Method for Computing the Mechanical Properties of Composites from Images of Their Microstructures, in: IUTAM symposium on microstructure-property interactions in composite materials, Springer, 1995, pp. 235–246

  24. [24]

    Zeman, J

    J. Zeman, J. Vondˇ rejc, J. Nov´ ak, I. Marek, Accelerating a FFT-Based Solver for Numerical Homoge- nization of Periodic Media by Conjugate Gradients, Journal of Computational Physics 229 (21) (2010) 8065–8071

  25. [25]

    M¨ uller, M

    V. M¨ uller, M. Kabel, H. Andr¨ a, T. B¨ ohlke, Homogenization of Linear Elastic Properties of Short-Fiber Reinforced Composites—A Comparison of Mean Field and Voxel-Based Methods, International Journal of Solids and Structures 67 (2015) 56–70

  26. [26]

    R. N. Glaesener, E. A. Tr¨ aff, B. Telgen, R. M. Canonica, D. M. Kochmann, Continuum Representation of Nonlinear Three-Dimensional Periodic Truss Networks by On-The-Fly Homogenization, International Journal of Solids and Structures 206 (2020) 101–113

  27. [27]

    P. Onck, E. Andrews, L. Gibson, Size Effects in Ductile Cellular Solids. Part I: Modeling, International Journal of Mechanical Sciences 43 (3) (2001) 681–699

  28. [28]

    Yoder, L

    M. Yoder, L. Thompson, J. Summers, Size Effects in Lattice Structures and a Comparison to Micropolar Elasticity, International Journal of Solids and Structures 143 (2018) 245–261

  29. [29]

    C. M. Portela, J. R. Greer, D. M. Kochmann, Impact of Node Geometry on the Effective Stiffness of Non-Slender Three-Dimensional Truss Lattice Architectures, Extreme Mechanics Letters 22 (2018) 138–148

  30. [30]

    Jamshidian, N

    M. Jamshidian, N. Boddeti, D. W. Rosen, O. Weeger, Multiscale Modelling of Soft Lattice Metamaterials: Micromechanical Nonlinear Buckling Analysis, Experimental Verification, and Macroscale Constitutive Behaviour, International Journal of Mechanical Sciences 188 (2020) 105956

  31. [31]

    Weeger, Isogeometric sizing and shape optimization of 3D beams and lattice structures at large deformations, Structural and Multidisciplinary Optimization 65 (2) (2022) 1–22

    O. Weeger, Isogeometric sizing and shape optimization of 3D beams and lattice structures at large deformations, Structural and Multidisciplinary Optimization 65 (2) (2022) 1–22

  32. [32]

    Bonatti, D

    C. Bonatti, D. Mohr, Large deformation response of additively-manufactured FCC metamaterials: From octet truss lattices towards continuous shell mesostructures, International Journal of Plasticity 92 (2017) 122–147

  33. [33]

    K. D. Nguyen, T. V. Duong, J. Lee, Y. Bazilevs, H. Nguyen-Xuan, An Efficient NURBS-Based Recon- struction Approach to Isogeometric Analysis of Triply Periodic Minimal Surface Structures, Engineering with Computers 41 (6) (2025) 4479–4509

  34. [34]

    Korshunova, G

    N. Korshunova, G. Alaimo, S. B. Hosseini, M. Carraturo, A. Reali, J. Niiranen, F. Auricchio, E. Rank, S. Kollmannsberger, Image-Based Numerical Characterization and Experimental Validation of Tensile Behavior of Octet-Truss Lattice Structures, Additive Manufacturing 41 (2021) 101949

  35. [35]

    Korshunova, G

    N. Korshunova, G. Alaimo, S. Hosseini, M. Carraturo, A. Reali, J. Niiranen, F. Auricchio, E. Rank, S. Kollmannsberger, Bending Behavior of Octet-Truss Lattice Structures: Modelling Options, Numerical Characterization and Experimental Validation, Materials & Design 205 (2021) 109693. 31

  36. [36]

    D¨ uster, J

    A. D¨ uster, J. Parvizian, Z. Yang, E. Rank, The Finite Cell Method for Three-Dimensional Problems of Solid Mechanics, Computer methods in applied mechanics and engineering 197 (45-48) (2008) 3768–3782

  37. [37]

    Hirschler, P

    T. Hirschler, P. Antolin, A. Buffa, Fast and Multiscale Formation of Isogeometric Matrices of Microstruc- tured Geometric Models, Computational Mechanics 69 (2) (2022) 439–466

  38. [38]

    Hirschler, R

    T. Hirschler, R. Bouclier, P. Antolin, A. Buffa, Reduced Order Modeling Based Inexact FETI-DP Solver for Lattice Structures, International Journal for Numerical Methods in Engineering 125 (8) (2024) e7419

  39. [39]

    Guillet, T

    C. Guillet, T. Hirschler, P. Jolivet, P. Antolin, R. Bouclier, Efficient Fine-Scale Simulation of Nonlinear Hyperelastic Lattice Structures, arXiv preprint arXiv:2603.10741 (2026)

  40. [40]

    Guillet, T

    C. Guillet, T. Hirschler, P. Jolivet, R. Bouclier, Multilevel Matrix-Free Method for High-Performance Isogeometric Analysis of Lattice Structures, Journal of Computational Physics 537 (2025) 114136

  41. [41]

    Kumar, S

    S. Kumar, S. Tan, L. Zheng, D. M. Kochmann, Inverse-Designed Spinodoid Metamaterials, npj Com- putational Materials 6 (1) (2020) 73

  42. [42]

    J. Feng, J. Fu, X. Yao, Y. He, Triply Periodic Minimal Surface (TPMS) Porous Structures: from Multi- Scale Design, Precise Additive Manufacturing to Multidisciplinary Applications, International Journal of Extreme Manufacturing 4 (2) (2022) 022001

  43. [43]

    Badia, F

    S. Badia, F. Verdugo, A. F. Mart´ ın, The Aggregated Unfitted Finite Element Method for Elliptic Problems, Computer Methods in Applied Mechanics and Engineering 336 (2018) 533–553

  44. [44]

    A. Main, G. Scovazzi, The Shifted Boundary Method for Embedded Domain Computations. Part I: Poisson and Stokes Problems, Journal of Computational Physics 372 (2018) 972–995

  45. [45]

    Burman, P

    E. Burman, P. Hansbo, M. G. Larson, S. Zahedi, Cut Finite Element Methods, Acta Numerica 34 (2025) 1–121

  46. [46]

    Chaturantabut, D

    S. Chaturantabut, D. C. Sorensen, Nonlinear Model Reduction via Discrete Empirical Interpolation, SIAM Journal on Scientific Computing 32 (5) (2010) 2737–2764

  47. [47]

    Negri, A

    F. Negri, A. Manzoni, D. Amsallem, Efficient Model Reduction of Parametrized Systems by Matrix Discrete Empirical Interpolation, Journal of Computational Physics 303 (2015) 431–454.doi:https: //doi.org/10.1016/j.jcp.2015.09.046

  48. [48]

    E. N. Karatzas, F. Ballarin, G. Rozza, Projection-Based Reduced Order Models for a Cut Finite Element Method in Parametrized Domains, Computers & Mathematics with Applications 79 (3) (2020) 833–851

  49. [49]

    Chasapi, P

    M. Chasapi, P. Antolin, A. Buffa, A Localized Reduced Basis Approach for Unfitted Domain Methods on Parameterized Geometries, Computer Methods in Applied Mechanics and Engineering 410 (2023) 115997

  50. [50]

    Mueller, S

    N. Mueller, S. Badia, Y. Zhao, Reduced Basis Solvers for Unfitted Methods on Parameterized Domains, Computer Methods in Applied Mechanics and Engineering 451 (2026) 118610

  51. [51]

    Parvizian, A

    J. Parvizian, A. D¨ uster, E. Rank, Finite Cell Method: h- and p-Extension for Embedded Domain Problems in Solid Mechanics, Computational Mechanics 41 (1) (2007) 121–133

  52. [52]

    Lehrenfeld, A

    C. Lehrenfeld, A. Reusken, Analysis of a High-Order Unfitted Finite Element Method for Elliptic Inter- face Problems, IMA Journal of Numerical Analysis 38 (3) (2018) 1351–1387

  53. [53]

    P. A. Martorell, S. Badia, High Order Unfitted Finite Element Discretizations for Explicit Boundary Representations, Journal of Computational Physics 511 (2024) 113127

  54. [54]

    Idesman, M

    A. Idesman, M. Mobin, J. Bishop, 10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes, Computational Mechanics 76 (2025) 1085–1115.doi:10.1007/s00466-025-02641-1. 32

  55. [55]

    Badia, F

    S. Badia, F. Verdugo, Robust and Scalable Domain Decomposition Solvers for Unfitted Finite Element Methods, Journal of Computational and Applied Mathematics 344 (2018) 740–759

  56. [56]

    Antolin, QUGaR: Quadratures for Unfitted GeometRies,https://github.com/pantolin/qugar, GitHub repository (2025)

    P. Antolin, QUGaR: Quadratures for Unfitted GeometRies,https://github.com/pantolin/qugar, GitHub repository (2025)

  57. [57]

    I. A. Baratta, J. P. Dean, J. S. Dokken, M. Habera, J. S. Hale, C. N. Richardson, M. E. Rognes, M. W. Scroggs, N. Sime, G. N. Wells, DOLFINx: The next generation FEniCS problem solving environment (Dec. 2023).doi:10.5281/zenodo.10447666

  58. [58]

    R. I. Saye, Algoim: Algorithms for Implicit Domains,https://github.com/algoim/algoim, GitHub repository (2025)

  59. [59]

    R. I. Saye, High-Order Quadrature Methods for Implicitly Defined Surfaces and Volumes in Hyperrect- angles, SIAM Journal on Scientific Computing 37 (2) (2015) A993–A1019

  60. [60]

    R. I. Saye, High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials, Journal of Computational Physics 448 (2022) 110720

  61. [61]

    Dauge, A

    M. Dauge, A. D¨ uster, E. Rank, Theoretical and Numerical Investigation of the Finite Cell Method, Journal of Scientific Computing 65 (3) (2015) 1039–1064

  62. [62]

    T. P. A. Mathew, Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations, Springer, 2008

  63. [63]

    Heinlein, A

    A. Heinlein, A. Klawonn, M. Lanser, J. Weber, A Frugal FETI-DP and BDDC Coarse Space for Hetero- geneous Problems, Tech. rep., TR series, Center for Data and Simulation Science, University of Cologne (2019)

  64. [64]

    Hirschler, P

    T. Hirschler, P. Antolin, A. Buffa, Fast and Multiscale Formation of Isogeometric Matrices of Mi- crostructured Geometric Models, Computational Mechanics 69 (2) (2022) 439–466.doi:10.1007/ s00466-021-02098-y

  65. [65]

    M. D. Mckay, R. J. Beckman, W. J. Conover, A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code, Technometrics 42 (1) (2000) 55–61. doi:10.1080/00401706.2000.10485979

  66. [66]

    Barrault, Y

    M. Barrault, Y. Maday, N. C. Nguyen, A. T. Patera, An ‘Empirical Interpolation’ Method: Application to Efficient Reduced-Basis Discretization of Partial Differential Equations, Comptes Rendus Mathema- tique 339 (9) (2004) 667–672.doi:https://doi.org/10.1016/j.crma.2004.08.006

  67. [67]

    M. J. Powell, The Theory of Radial Basis Function Approximation in 1990, Advances in numerical analysis (1992) 105–210

  68. [68]

    Bungartz, M

    H.-J. Bungartz, M. Griebel, Sparse Grids, Acta Numerica, Cambridge University Press, 2004, pp. 147– 270

  69. [69]

    Bonilla Moreno, FLAS h: Fast simulation tools for Lattice Structures,https://github.com/ Gonzalobm99/FLASh, GitHub repository (2025)

    G. Bonilla Moreno, FLAS h: Fast simulation tools for Lattice Structures,https://github.com/ Gonzalobm99/FLASh, GitHub repository (2025)

  70. [70]

    Balay, S

    S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, E. M. Constantinescu, L. Dalcin, S. Benson, A. Dener, V. Eijkhout, J. Faibussowitsch, W. D. Gropp, V. Hapla, T. Isaac, P. Jolivet, D. Karpeev, D. Kaushik, M. G. Knepley, F. Kong, S. Kruger, D. A. May, L. C. McInnes, R. T. Mills, L. Mitchell, T. Munson, J. E. Roman, K. Rupp, P. Sanan, ...

  71. [71]

    L. D. Dalcin, R. R. Paz, P. A. Kler, A. Cosimo, Parallel distributed computing using Python, Advances in Water Resources 34 (9) (2011) 1124–1139.doi:10.1016/j.advwatres.2011.04.013. 33

  72. [72]

    Antolin, G

    P. Antolin, G. Bonilla Moreno, Reduced-Order Model data for FLASh: Fast simulation tools for Lattice Structures, dataset (Mar. 2026).doi:10.5281/zenodo.19254389

  73. [73]

    Burman, S

    E. Burman, S. Claus, P. Hansbo, M. G. Larson, A. Massing, CutFEM: Discretizing Geometry and Partial Differential Equations, International Journal for Numerical Methods in Engineering 104 (7) (2015) 472– 501

  74. [74]

    D. M. Young, Iterative Solution of Large Linear Systems, Elsevier, 2014. 34