pith. sign in

arxiv: 2605.31532 · v1 · pith:T4YRYZXDnew · submitted 2026-05-29 · ❄️ cond-mat.soft · cs.LG

Discovering Thermodynamically Admissible Dissipation Potentials via Grammar-Based Symbolic Regression

Pith reviewed 2026-06-28 19:52 UTC · model grok-4.3

classification ❄️ cond-mat.soft cs.LG
keywords symbolic regressiondissipation potentialsthermodynamic admissibilitygeneralized standard materialsviscoelasticityviscoplasticityelastomeroscillatory shear
0
0 comments X

The pith

A grammar-based symbolic regression discovers thermodynamically admissible dissipation potentials for inelastic materials.

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

The paper develops a symbolic regression method to recover dissipation potentials that drive the evolution of internal variables in the generalized standard materials framework. It enforces the Clausius-Duhem inequality by requiring the dual dissipation potential to be convex and non-negative, conditions that guarantee non-negative mechanical dissipation. These constraints are built directly into a composition-extended grammar so that every generated candidate automatically satisfies thermodynamic admissibility. The approach recovers known ground-truth behaviors from noisy synthetic data and, on real oscillatory shear tests of a synthetic elastomer, produces models that capture amplitude-dependent softening of the dynamic moduli while outperforming a calibrated linear Zener model.

Core claim

Candidate dual dissipation potentials are generated by a composition-extended convexity-preserving grammar that guarantees convexity and non-negativity by construction, thereby satisfying the subdifferential form of the Clausius-Duhem inequality and enabling unified discovery of rate-dependent viscoelastic and viscoplastic mechanisms, including those with genuine elastic domains.

What carries the argument

The composition-extended convexity-preserving grammar, which produces only convex and non-negative expressions for the dual dissipation potential.

If this is right

  • The method unifies modeling of viscoelastic and viscoplastic dissipative mechanisms under one thermodynamically consistent framework.
  • Discovered potentials reproduce amplitude-dependent softening of dynamic moduli observed in experimental oscillatory shear data.
  • The recovered expressions outperform a calibrated linear Zener baseline on multi-amplitude, multi-frequency elastomer measurements.
  • Performance holds under added process and measurement noise on synthetic Newtonian, power-law, and Bingham viscoplastic cases.

Where Pith is reading between the lines

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

  • Extending the grammar rules could reach a broader class of material behaviors without losing the built-in admissibility guarantee.
  • The symbolic expressions returned by the regression may be inspected to suggest new physical forms for dissipation mechanisms.
  • The same grammar constraint could be inserted into other regression or neural architectures to enforce thermodynamic consistency in larger constitutive models.

Load-bearing premise

The true dissipation potential must be expressible by expressions allowed inside the composition-extended convexity-preserving grammar.

What would settle it

Apply the method to a material whose dissipation potential requires a functional form outside the grammar and check whether recovery fails or produces inadmissible models despite abundant data.

Figures

Figures reproduced from arXiv: 2605.31532 by Federico Califano, Jacopo Ciambella.

Figure 1
Figure 1. Figure 1: Expression tree representation of a candidate dissipation potential [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the GP evolution loop: candidate potentials from the composition [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Newtonian viscosity benchmark (E1): heatmaps of mean DMA relative error across [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Newtonian viscosity benchmark (E1): strain–stress hysteresis loops for representative [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: DMA moduli comparison for the Newtonian viscosity benchmark (E1), noise-free [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Abs-power benchmark (E2): heatmaps of mean DMA relative error across the noise [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Abs-power benchmark (E2): strain–stress hysteresis loops under noise-free (top row) [PITH_FULL_IMAGE:figures/full_fig_p024_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: DMA moduli comparison for the abs-power benchmark (E2), noise-free case, trial 1. [PITH_FULL_IMAGE:figures/full_fig_p025_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: E3 Bingham viscoplastic: strain–stress hysteresis loops under the LP2 triangular [PITH_FULL_IMAGE:figures/full_fig_p026_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Representative experimental acquisition (nominal strain amplitude [PITH_FULL_IMAGE:figures/full_fig_p027_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Schematic of the identification strategy used for experimental data. Unlike the [PITH_FULL_IMAGE:figures/full_fig_p028_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Shear hysteresis loops on the synthetic elastomer specimen. Rows: nominal strain [PITH_FULL_IMAGE:figures/full_fig_p030_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Dynamic moduli of the synthetic elastomer specimen extracted from the first har [PITH_FULL_IMAGE:figures/full_fig_p031_13.png] view at source ↗
read the original abstract

Constitutive laws for inelastic materials must satisfy strict thermodynamic admissibility requirements, yet current data-driven approaches sacrifice interpretability, even when formal guarantees are provided by physics-encoded architectures. We propose a symbolic regression framework for the data-driven discovery of dissipation potentials governing the evolution of internal variables within the Generalized Standard Materials (GSM) formalism. Starting from the Clausius--Duhem inequality, we enforce the thermodynamic requirements, convexity and non-negativity, that the dual dissipation potential must satisfy to guarantee non-negative mechanical dissipation. These requirements are formulated in the general subdifferential setting, encompassing rate-dependent (viscoelastic) and viscoplastic dissipative mechanisms, including potentials with genuine elastic domains, within a unified framework. Candidate potentials are generated by a composition-extended convexity-preserving grammar that guarantees thermodynamic admissibility \emph{by construction}. The framework is validated on synthetic datasets spanning Newtonian, power-law, and Bingham viscoplastic ground truths under process and measurement noise, and on experimental oscillatory shear measurements of a synthetic elastomer across multiple strain amplitudes and frequencies, where the discovered potentials reproduce the amplitude-dependent softening of the dynamic moduli and outperform a calibrated linear Zener baseline.

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 / 2 minor

Summary. The paper proposes a symbolic regression framework that uses a composition-extended convexity-preserving grammar to discover dual dissipation potentials for Generalized Standard Materials models. Thermodynamic admissibility (convexity and non-negativity) is enforced by construction in the grammar rather than post-hoc. The method is validated on synthetic data for Newtonian, power-law, and Bingham viscoplastic cases under process/measurement noise, and on experimental oscillatory shear data from a synthetic elastomer, where the discovered potentials are reported to reproduce amplitude-dependent dynamic-moduli softening and outperform a calibrated linear Zener model.

Significance. If the central empirical claims hold, the approach would provide an interpretable, thermodynamically guaranteed route to data-driven constitutive modeling that unifies rate-dependent and viscoplastic mechanisms. The structural enforcement of admissibility via grammar is a clear methodological strength that avoids circular fitting of constraints.

major comments (3)
  1. [§3.2] §3.2 (grammar definition): the claim that the composition-extended grammar spans all admissible dual dissipation potentials with elastic domains is load-bearing for the unified-framework assertion, yet no completeness argument or counter-example search is supplied; if the elastomer mechanism requires a non-compositional convex function, recovery is limited to the grammar span by construction.
  2. [§5.3] §5.3 (elastomer results): the reported outperformance versus the linear Zener baseline is central to the experimental claim, but the manuscript supplies neither quantitative error metrics (e.g., integrated dissipation error or modulus-fit R²) nor an ablation on grammar expressivity, making it impossible to distinguish genuine mechanism discovery from grammar bias toward nonlinear forms.
  3. [§4.2] §4.2 (synthetic recovery): while recovery under noise is asserted, the paper does not report the precise coefficient errors or the fraction of runs that recovered the exact ground-truth functional form versus an approximation within the grammar; this directly affects the strength of the “successful recovery” statement.
minor comments (2)
  1. [Abstract] Abstract: the phrase “better performance than Zener” should be accompanied by at least one numerical metric for immediate context.
  2. [Results] Notation: the subdifferential setting is introduced but the precise mapping from discovered potential to the evolution equation (Eq. (X)) is not restated in the results sections, which would aid readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which help clarify the strengths and limitations of our approach. We address each major comment below and will incorporate revisions to improve the manuscript.

read point-by-point responses
  1. Referee: [§3.2] §3.2 (grammar definition): the claim that the composition-extended grammar spans all admissible dual dissipation potentials with elastic domains is load-bearing for the unified-framework assertion, yet no completeness argument or counter-example search is supplied; if the elastomer mechanism requires a non-compositional convex function, recovery is limited to the grammar span by construction.

    Authors: We agree that the manuscript does not supply a formal completeness argument or counter-example search. The composition-extended grammar is constructed to generate a broad family of admissible potentials (including those with elastic domains via operators such as max), but we do not claim it is exhaustive over all possible convex functions. This is an inherent limitation of any fixed grammar. We will revise §3.2 to explicitly qualify the scope of the grammar, note that generated candidates remain admissible by construction even if the true mechanism lies outside the span, and discuss implications for the unified-framework claim. revision: partial

  2. Referee: [§5.3] §5.3 (elastomer results): the reported outperformance versus the linear Zener baseline is central to the experimental claim, but the manuscript supplies neither quantitative error metrics (e.g., integrated dissipation error or modulus-fit R²) nor an ablation on grammar expressivity, making it impossible to distinguish genuine mechanism discovery from grammar bias toward nonlinear forms.

    Authors: The referee correctly identifies the lack of quantitative metrics and ablation. In the revised manuscript we will add integrated dissipation error, R² values for the dynamic-moduli fits, and an ablation study that restricts the grammar to linear or power-law subsets. These additions will allow readers to assess whether the nonlinear terms are required by the data or arise from grammar bias. revision: yes

  3. Referee: [§4.2] §4.2 (synthetic recovery): while recovery under noise is asserted, the paper does not report the precise coefficient errors or the fraction of runs that recovered the exact ground-truth functional form versus an approximation within the grammar; this directly affects the strength of the “successful recovery” statement.

    Authors: We acknowledge that precise coefficient errors and recovery fractions across runs are not reported. The revised §4.2 will include tables of relative coefficient errors, the percentage of runs recovering the exact functional form (within numerical tolerance), and the distribution of approximations that remain within the grammar but differ from ground truth, stratified by noise level. revision: yes

Circularity Check

0 steps flagged

No significant circularity; method enforces constraints structurally and fits data independently

full rationale

The paper constructs candidate dissipation potentials via a composition-extended convexity-preserving grammar that enforces thermodynamic admissibility (convexity, non-negativity) by design from the Clausius-Duhem inequality, then applies symbolic regression to identify coefficients from data. This is a generative search plus fitting procedure, not a derivation in which any claimed result (e.g., reproduction of amplitude-dependent softening or outperformance versus Zener) reduces by construction to the inputs or to a self-citation chain. Validation on both synthetic ground-truth cases and experimental elastomer data is standard empirical testing of the fitted model; the grammar's span limitation is an expressivity bound, not a circular reduction. No load-bearing step equates a prediction to its own fitted parameters or imports uniqueness via overlapping-author citations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The framework rests on the standard Clausius-Duhem inequality and the mathematical requirement that dual dissipation potentials be convex and non-negative; the grammar itself is an invented construct to enforce these properties during candidate generation. No explicit free parameters are named in the abstract beyond those implicitly fitted during regression.

axioms (2)
  • domain assumption Clausius-Duhem inequality must hold for all admissible processes
    Invoked at the start to derive the requirements on the dissipation potential
  • domain assumption Dual dissipation potential must be convex and non-negative in the subdifferential sense
    Central constraint used to define thermodynamic admissibility for both rate-dependent and viscoplastic cases
invented entities (1)
  • composition-extended convexity-preserving grammar no independent evidence
    purpose: Generate only thermodynamically admissible candidate potentials by construction
    New structure introduced to restrict the search space of symbolic regression

pith-pipeline@v0.9.1-grok · 5729 in / 1331 out tokens · 20762 ms · 2026-06-28T19:52:50.971425+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

62 extracted references · 4 canonical work pages · 3 internal anchors

  1. [1]

    Ciambella, M

    J. Ciambella, M. Destrade, and R. W. Ogden. On the ABAQUS FEA model of finite viscoelasticity.Rubber Chemistry and Technology, 82(2):184–193, 2009

  2. [2]

    Ciambella and P

    J. Ciambella and P. Nardinocchi. A structurally frame-indifferent model for anisotropic visco-hyperelastic materials.Journal of the Mechanics and Physics of Solids, 147:104247, 2021

  3. [3]

    A. R. Payne. The dynamic properties of carbon black-loaded natural rubber vulcanizates. Part I.Journal of Applied Polymer Science, 6(19):57–63, 1962

  4. [4]

    A. R. Payne. The dynamic properties of carbon black-loaded natural rubber vulcanizates. Part II.Journal of Applied Polymer Science, 7(3):873–885, 1963

  5. [5]

    Kabliman, A

    E. Kabliman, A. H. Kolody, J. Kronsteiner, M. Kommenda, and G. Kronberger. Applica- tion of symbolic regression for constitutive modeling of plastic deformation.Applications in Engineering Science, 6:100052, 2021

  6. [6]

    R. C. Batra and C. H. Kim. Effect of viscoplastic flow rules on the initiation and growth of shear bands at high strain rates.Journal of the Mechanics and Physics of Solids, 38(6):859– 874, 1990

  7. [7]

    Holzapfel

    Gerhard A. Holzapfel. Biomechanics of soft tissue. InHandbook of Materials Behavior Models, pages 1057–1071. Academic Press, 2001

  8. [8]

    Fundamental problems in viscoplasticity.Advances in Applied Mechanics, 9:243–377, 1966

    Piotr Perzyna. Fundamental problems in viscoplasticity.Advances in Applied Mechanics, 9:243–377, 1966

  9. [9]

    J. N. Fuhg, G. A. Padmanabha, N. Bouklas, B. Bahmani, W. Sun, N. N. Vlassis, M. Flaschel, and L. De Lorenzis. A review on data-driven constitutive laws for solids. Archives of Computational Methods in Engineering, 2024. 35

  10. [10]

    Makke and S

    N. Makke and S. Chawla. Interpretable scientific discovery with symbolic regression: a review.Artificial Intelligence Review, 57(1):2, 2024

  11. [11]

    J. R. Koza.Genetic Programming: On the Programming of Computers by Means of Natural Selection. MIT Press, 1992

  12. [12]

    Schmidt and H

    M. Schmidt and H. Lipson. Distilling free-form natural laws from experimental data. Science, 324(5923):81–85, 2009

  13. [13]

    Interpretable Machine Learning for Science with PySR and SymbolicRegression.jl

    Miles Cranmer. Interpretable machine learning for science with PySR and SymbolicRe- gression.jl.arXiv preprint arXiv:2305.01582, 2023

  14. [14]

    Abdusalamov, M

    R. Abdusalamov, M. Hillg¨ artner, and M. Itskov. Automatic generation of interpretable hyperelastic material models by symbolic regression.International Journal for Numerical Methods in Engineering, 124(9):2093–2104, 2023

  15. [15]

    J. Hou, X. Chen, T. Wu, E. Kuhl, and X. Wang. Automated data-driven discovery of material models based on symbolic regression: A case study on the human brain cortex. Acta Biomaterialia, 188:276–296, 2024

  16. [16]

    Kissas, S

    G. Kissas, S. Mishra, E. Chatzi, and L. De Lorenzis. The language of hyperelastic materials. Computer Methods in Applied Mechanics and Engineering, 428:117053, 2024

  17. [17]

    G. F. Bomarito, J. D. Hochhalter, T. J. Rubin, and K. A. Townsend. Development of inter- pretable, data-driven plasticity models with symbolic regression.Computers & Structures, 252:106557, 2021

  18. [18]

    Garbrecht, B

    K. Garbrecht, B. Birky, B. T. Lester, J. M. Emery, and J. D. Hochhalter. Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression. Journal of the Mechanics and Physics of Solids, 173:105198, 2023

  19. [19]

    Evgeniya Kabliman and Gabriel Kronberger. Identification of empirical constitutive models for age-hardenable aluminium alloy and high-chromium martensitic steel using symbolic regression.Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 384(2317), 2026

  20. [20]

    Y. Zhu, C. Su, and X. Li. A physics-guided symbolic regression framework for efficient and interpretable sand constitutive modeling.Canadian Geotechnical Journal, 62:1–14, 2025

  21. [21]

    Bahmani, H

    B. Bahmani, H. S. Suh, and W. Sun. Discovering interpretable elastoplasticity models via the neural polynomial method enabled symbolic regressions.Computer Methods in Applied Mechanics and Engineering, 422:116827, 2024

  22. [22]

    Bahmani and W

    B. Bahmani and W. Sun. Physics-constrained symbolic model discovery for polyconvex incompressible hyperelastic materials.International Journal for Numerical Methods in Engineering, 125:e7473, 2024

  23. [23]

    Halphen and Q

    B. Halphen and Q. S. Nguyen. Sur les mat´ eriaux standards g´ en´ eralis´ es.Journal de M´ ecanique, 14(1):39–63, 1975

  24. [24]

    Lemaitre and J.-L

    J. Lemaitre and J.-L. Chaboche.Mechanics of Solid Materials. Cambridge University Press, 2000

  25. [25]

    Kalina, J¨ org Brummund, WaiChing Sun, and Markus K¨ astner

    Max Rosenkranz, Karl A. Kalina, J¨ org Brummund, WaiChing Sun, and Markus K¨ astner. Viscoelasticty with physics-augmented neural networks: model formulation and training methods without prescribed internal variables.Computational Mechanics, 74(6):1279–1301, 2024. 36

  26. [26]

    A mechanics-informed neural network framework for data-driven nonlinear viscoelasticity

    Faisal As’ad and Charbel Farhat. A mechanics-informed neural network framework for data-driven nonlinear viscoelasticity. InAIAA SCITECH 2023 Forum, page 0949, 2023

  27. [27]

    Theory and implementation of inelastic constitutive artificial neural networks.Computer Methods in Applied Mechanics and Engineering, 428:117063, 2024

    Hagen Holthusen, Lukas Lamm, Tim Brepols, Stefanie Reese, and Ellen Kuhl. Theory and implementation of inelastic constitutive artificial neural networks.Computer Methods in Applied Mechanics and Engineering, 428:117063, 2024

  28. [28]

    Rausch, Francisco Sahli Costabal, and Adrian Buganza Tepole

    Vahidullah Ta¸ c, Manuel K. Rausch, Francisco Sahli Costabal, and Adrian Buganza Tepole. Data-driven anisotropic finite viscoelasticity using neural ordinary differential equations. Computer Methods in Applied Mechanics and Engineering, 411:116046, 2023

  29. [29]

    Abdolazizi, Kevin Linka, and Christian J

    Kian P. Abdolazizi, Kevin Linka, and Christian J. Cyron. Viscoelastic constitutive arti- ficial neural networks (vcanns) – a framework for data-driven anisotropic nonlinear finite viscoelasticity.Journal of Computational Physics, 499:112704, 2024

  30. [30]

    Califano and J

    F. Califano and J. Ciambella. Enhancing nonlinear viscoelastic modeling of elastomers through neural networks: A deep rheological element.Mechanics of Materials, 212:105525, 2026

  31. [31]

    Kalina, J¨ org Brummund, and Markus K¨ astner

    Karl A. Kalina, J¨ org Brummund, and Markus K¨ astner. A physics-augmented neural net- work framework for finite strain incompressible viscoelasticity.Computer Methods in Ap- plied Mechanics and Engineering, 455:118892, 2026

  32. [32]

    Convex neural networks learn generalized standard material models.Journal of the Mechanics and Physics of Solids, 200:106103, 2025

    Moritz Flaschel, Paul Steinmann, Laura De Lorenzis, and Ellen Kuhl. Convex neural networks learn generalized standard material models.Journal of the Mechanics and Physics of Solids, 200:106103, 2025

  33. [33]

    Boes, J.-W

    B. Boes, J.-W. Simon, and H. Holthusen. Accounting for plasticity: An extension of inelastic constitutive artificial neural networks.European Journal of Mechanics - A/Solids, 117:105998, 2026

  34. [34]

    van der Velden, B

    T. van der Velden, B. Boes, T. Brepols, E. Kuhl, and H. Holthusen. A note on constitutive artificial neural networks for finite strain plasticity.Mechanics Research Communications, 155:104707, 2026

  35. [35]

    Jadoon, Knut Andreas Meyer, and Jan N

    Asghar A. Jadoon, Knut Andreas Meyer, and Jan N. Fuhg. Automated model discovery of finite strain elastoplasticity from uniaxial experiments.Computer Methods in Applied Mechanics and Engineering, 435:117653, 2025

  36. [36]

    Abdolazizi, Hagen Holthusen, Christian J

    Chenyi Ji, Kian P. Abdolazizi, Hagen Holthusen, Christian J. Cyron, and Kevin Linka. Inelastic constitutive kolmogorov–arnold networks: a generalized framework for automated discovery of interpretable inelastic material models.arXiv preprint arXiv:2602.17750, 2026

  37. [37]

    Z. Liu, Y. Wang, S. Vaidya, et al. KAN: Kolmogorov–Arnold Networks.arXiv preprint arXiv:2404.19756, 2024

  38. [38]

    Thakolkaran, Y

    P. Thakolkaran, Y. Guo, S. Saini, M. Peirlinck, B. Alheit, and S. Kumar. Can KAN CANs? Input-convex Kolmogorov–Arnold Networks (KANs) as hyperelastic constitutive artificial neural networks (CANs).Computer Methods in Applied Mechanics and Engineering, 2025

  39. [39]

    Califano and J

    F. Califano and J. Ciambella. Enforcing physics in hyperelasticity modeling using Kolmogorov–Arnold networks. InProceedings of the XXVI AIMETA Conference 2024, pages 920–928. Springer, Cham, 2026

  40. [40]

    Germain, Q

    P. Germain, Q. S. Nguyen, and P. Suquet. Continuum thermodynamics.Journal of Applied Mechanics, 50(4b):1010–1020, 1983. 37

  41. [41]

    G. K. Batchelor. An introduction to fluid dynamics.Cambridge University Press, 1967

  42. [42]

    W. Ostwald. ¨Uber die Geschwindigkeitsfunktion der Viskosit¨ at disperser Systeme. I. Kolloid-Zeitschrift, 36(2):99–117, 1925

  43. [43]

    E. C. Bingham.Fluidity and Plasticity. McGraw-Hill, 1922

  44. [44]

    H. Eyring. Viscosity, plasticity, and diffusion as examples of absolute reaction rates.The Journal of Chemical Physics, 4(4):283–291, 1936

  45. [45]

    Viscoplastic simple shear at finite strains.Proc

    F Califano and J Ciambella. Viscoplastic simple shear at finite strains.Proc. R. Soc. A Math. Phys. Eng. Sci., 479(2280), dec 2023

  46. [46]

    Boyd and L

    S. Boyd and L. Vandenberghe.Convex Optimization. Cambridge University Press, 2004

  47. [47]

    Tyrrell Rockafellar.Convex Analysis

    R. Tyrrell Rockafellar.Convex Analysis. Princeton University Press, Princeton, NJ, 1970

  48. [48]

    B. Amos, L. Xu, and J. Z. Kolter. Input convex neural networks. InProceedings of the 34th International Conference on Machine Learning, pages 146–155, 2017

  49. [49]

    Dugas, Y

    C. Dugas, Y. Bengio, F. B´ elisle, C. Nadeau, and R. Garcia. Incorporating second-order functional knowledge for better option pricing. InProceedings of the 13th International Conference on Neural Information Processing Systems (NIPS’00), pages 451–457. MIT Press, 2000

  50. [50]

    La Cava, B

    W. La Cava, B. Burlacu, M. Virgolin, et al. Contemporary symbolic regression methods and their relative performance.Advances in Neural Information Processing Systems, 2021

  51. [51]

    Kommenda, B

    M. Kommenda, B. Burlacu, G. Kronberger, and M. Affenzeller. Parameter identification for symbolic regression using nonlinear least squares.Genetic Programming and Evolvable Machines, 21:471–501, 2020

  52. [52]

    Springer Science & Business Media, 1996

    Heinz Werner Engl, Martin Hanke, and Andreas Neubauer.Regularization of Inverse Problems, volume 375. Springer Science & Business Media, 1996

  53. [53]

    R. H. Byrd, P. Lu, J. Nocedal, and C. Zhu. A limited memory algorithm for bound constrained optimization.SIAM Journal on Scientific Computing, 16(5):1190–1208, 1995

  54. [54]

    C. Zhu, R. H. Byrd, P. Lu, and J. Nocedal. Algorithm 778: L-BFGS-B: Fortran subrou- tines for large-scale bound-constrained optimization.ACM Transactions on Mathematical Software, 23(4):550–560, 1997

  55. [55]

    Nelder and Roger Mead

    John A. Nelder and Roger Mead. A simplex method for function minimization.The Computer Journal, 7(4):308–313, 1965

  56. [56]

    Implementing the nelder–mead simplex algorithm with adaptive parameters.Computational Optimization and Applications, 51(1):259–277, 2012

    Fuchang Gao and Lixing Han. Implementing the nelder–mead simplex algorithm with adaptive parameters.Computational Optimization and Applications, 51(1):259–277, 2012

  57. [57]

    Colton and R

    D. Colton and R. Kress.Inverse Acoustic and Electromagnetic Scattering Theory. Springer, New York, 3rd edition, 2013

  58. [58]

    A. Wirgin. The inverse crime, 2004. arXiv preprint math-ph/0401050

  59. [59]

    Metropolis and S

    N. Metropolis and S. Ulam. The Monte Carlo method.Journal of the American Statistical Association, 44(247):335–341, 1949

  60. [60]

    G. E. Uhlenbeck and L. S. Ornstein. On the theory of the Brownian motion.Physical Review, 36(5):823–841, 1930. 38

  61. [61]

    Food and Drug Administration

    U.S. Food and Drug Administration. Artificial intelligence/machine learning (AI/ML)- based software as a medical device (SaMD) action plan. Technical report, FDA, 2021

  62. [62]

    American Society of Mechanical Engineers, 2018

    ASME.V&V 40–2018: Assessing Credibility of Computational Modeling through Verifi- cation and Validation: Application to Medical Devices. American Society of Mechanical Engineers, 2018. 39