On the Risk of Small Networks for Physics-Knowledgeable Studying

Introduction

within the interval of 2017-2019, physics-informed neural networks (PINNs) have been a very fashionable space of analysis within the scientific machine studying (SciML) group [1,2]. PINNs are used to resolve strange and partial differential equations (PDEs) by representing the unknown answer discipline with a neural community, and discovering the weights and biases (parameters) of the community by minimizing a loss perform primarily based on the governing differential equation. For instance, the unique PINNs strategy penalizes the sum of pointwise errors of the governing PDE, whereas the Deep Ritz technique minimizes an “power” useful whose minimal enforces the governing equation [3]. One other various is to discretize the answer with a neural community, then assemble the weak type of the governing equation utilizing adversarial networks [4] or polynomial take a look at features [5,6]. Whatever the selection of physics loss, neural community discretizations have been efficiently used to investigate quite a few methods ruled by PDEs. From the Navier-Stokes equations [7] to conjugate warmth switch [8] and elasticity [9], PINNs and their variants have confirmed themselves a worthy addition to the computational scientist’s toolkit.

As is the case with all machine studying issues, an important ingredient of acquiring strong and correct options with PINNs is hyperparameter tuning. The analyst has a lot freedom in developing an answer technique—as mentioned above, the selection of physics loss perform is just not distinctive, neither is the optimizer, strategy of boundary situation enforcement, or neural community structure. For instance, whereas ADAM has traditionally been the go-to optimizer for machine studying issues, there was a surge of curiosity in second-order Newton-type strategies for physics-informed issues [10,11]. Different research have in contrast methods for implementing boundary circumstances on the neural community discretization [12]. Greatest practices for the PINNs structure have primarily been investigated by the selection of activation features. To fight the spectral bias of neural networks [13], sinusoidal activation features have been used to raised symbolize high-frequency answer fields [14,15]. In [16], quite a few normal activation features have been in contrast on compressible fluid move issues. Activation features with partially learnable options have been proven to enhance answer accuracy in [17]. Whereas most PINNs depend on multi-layer perceptron (MLP) networks, convolutional networks have been investigated in [18] and recurrent networks in [19].

The research referenced above are removed from an exhaustive listing of works investigating the selection of hyperparameters for physics-informed coaching. Nevertheless, these research show that the loss perform, optimizer, activation perform, and primary class of community structure (MLP, convolutional, recurrent, and many others.) have all obtained consideration within the literature as attention-grabbing and necessary elements of the PINN answer framework. One hyperparameter that has seen comparatively little scrutiny is the dimensions of the neural community discretizing the answer discipline. In different phrases, to the perfect of our information, there aren’t any printed works that ask the next query: what number of parameters ought to the physics-informed community encompass? Whereas this query is in some sense apparent, the group’s lack of curiosity in it’s not shocking—there isn’t any value to pay in answer accuracy for an overparameterized community. In actual fact, overparameterized networks can present useful regularization of the answer discipline, as is seen with the phenomenon of double descent [20]. Moreover, within the context of data-driven classification issues, overparameterized networks have been proven to result in smoother loss landscapes [21]. As a result of answer accuracy offers no incentive to drive the dimensions of the community down, and since the optimization drawback may very well favor overparameterization, many authors use very giant networks to symbolize PDE options.

Whereas the accuracy of the answer solely stands to achieve by growing the parameter depend, the computational price of the answer does scale with the community measurement. On this research, we take three examples from the PINNs literature and present that networks with orders of magnitude fewer parameters are able to satisfactorily reproducing the outcomes of the bigger networks. The conclusion from these three examples is that, within the case of low-frequency answer fields, small networks can acquire correct options with decreased computational price. We then present a counterexample, the place regression to a fancy oscillatory perform repeatedly advantages from growing the community measurement. Thus, our suggestion is as follows: the variety of parameters in a PINN needs to be as few as attainable, however no fewer.

Examples

The primary three examples are impressed by issues taken from the PINNs literature. In these works, giant networks are used to acquire the PDE answer, the place the dimensions of the community is measured by the variety of parameters. Whereas totally different community architectures could carry out in a different way with totally different parameter counts, we use this metric as a proxy for community complexity, unbiased of the structure. In our examples, we incrementally lower the parameter depend of a multilayer perceptron community till the error with a reference answer begins to extend. This level represents a decrease restrict on the community measurement for the actual drawback, and we examine the variety of parameters at this level to the variety of parameters used within the authentic paper. In every case, we discover that the networks from the literature are overparameterized by at the least an order of magnitude. Within the fourth instance, we clear up a regression drawback to point out how small networks can fail to symbolize oscillatory fields, which acts as a caveat to our findings.

Part discipline fracture

The part discipline mannequin of fracture is a variational strategy to fracture mechanics, which concurrently finds the displacement and harm fields by minimizing a suitably outlined power useful [22]. Our research is predicated on the one-dimensional instance drawback given in [23], which makes use of the Deep Ritz technique to find out the displacement and harm fields that decrease the fracture power useful. This power useful is given by

[PiBig(u(x),alpha(x)Big) = Pi^u + Pi^{alpha}=int_0^1 frac{1}{2}(1-alpha)^2 Big(frac{partial u}{partial x}Big)^2 + frac{3}{8}Big( alpha + ell^2 Big(frac{partial alpha}{partial x}Big)^2 Big) dx,]

the place ( x ) is the spatial coordinate, ( u(x) ) is the displacement, ( alpha(x)in[0,1] ) is the crack density, and ( ell ) is a size scale figuring out the width of smoothed cracks. The power useful includes two elements ( Pi^u ) and ( Pi^{alpha} ), that are the elastic and fracture energies respectively. As within the cited work, we take (ell=0.05). The displacement and part discipline are discretized with a single neural community ( N: mathbb R rightarrow mathbb R^2 ) with parameters (boldsymbol theta). The issue is pushed by an utilized tensile displacement on the appropriate finish, which we denote ( U ). Boundary circumstances are constructed into the 2 fields with

[ begin{bmatrix}
u(x ; boldsymbol theta) alpha (x; boldsymbol theta)
end{bmatrix} = begin{bmatrix} x(1-x) N_1(x;boldsymbol theta) + Ux x(1-x) N_2(x; boldsymbol theta)
end{bmatrix} ,]

the place (N_i ) refers back to the (i)-th output of the community and the Dirichlet boundary circumstances on the crack density are used to suppress cracking on the edges of the area. In [23], a 4 hidden-layer MLP community with a width of fifty is used to symbolize the 2 answer fields. If we neglect the bias on the remaining layer, this corresponds to (7850) trainable parameters. For all of our research, we use a two-hidden layer community with hyperbolic tangent activation features and no bias on the output layer, as, in our expertise, these networks suffice to symbolize any answer discipline of curiosity. If each hidden layers have width (M), the full variety of parameters on this community is (M^2+5M). When (M=86), we acquire (7826) trainable parameters. Within the absence of an analytical answer, we use this as the big community reference answer to which the smaller networks are in contrast.

To generate the answer fields, we decrease the full potential power utilizing ADAM optimization with a studying fee of (1 occasions 10^{-3}). Complete fracture of the bar is noticed round (U=0.6). As within the paper, we compute the elastic and fracture energies over a variety of utilized displacements, the place ADAM is run for (3500) epochs at every displacement increment to acquire the answer fields. These “loading curves” are used to match the efficiency of networks of various sizes. Our experiment is carried out with (8) totally different community sizes, every comprising (20) increments of the utilized displacement to construct the loading curves. See Determine 1 for the loading curves computed with the totally different community sizes. Solely when there are (|boldsymbol theta|=14) parameters, which corresponds to a community of width (2), will we see a divergence from the loading curves of the big community reference answer. The smallest community that performs nicely has (50) parameters, which is (157times) smaller than the community used within the paper. Determine 2 confirms that this small community is able to approximating the discontinuous displacement discipline, in addition to the localized harm discipline.

Burgers’ equation

We now research the impact of community measurement on a standard mannequin drawback from fluid mechanics. Burgers’ equation is steadily used to check numerical answer strategies due to its nonlinearity and tendency to type sharp options. The viscous Burgers’ equation with homogeneous Dirichlet boundaries is given by

[frac{partial u}{partial t} + ufrac{partial u}{partial x} = nu frac{partial^2 u}{partial x^2}, quad u(x,0) = u_0(x), quad u(-1,t)=u(1,t) = 0,]

the place (xin[-1,1]) is the spatial area, (tin[0,T]) is the time coordinate, (u(x,t)) is the rate discipline, (nu) is the viscosity, and (u_0(x)) is the preliminary velocity profile. In [24], a neural community discretization of the rate discipline is used to acquire an answer to the governing differential equation. Their community comprises (3) hidden layers with (64) neurons in every layer, comparable to (8576) trainable parameters. Once more, we use a two hidden-layer community that has (5M+M^2) trainable parameters the place (M) is the width of every layer. If we take (M=90), we acquire (8550) trainable parameters in our community. We take the answer from this community to be the reference answer (u_{textual content{ref}}(x,t)), and compute the discrepancy between velocity fields from smaller networks. We do that with an error perform given by

[ E Big( u(x,t)Big)= frac{int_{Omega}| u(x,t) – u_{text{ref}}(x,t)| dOmega}{int_{Omega}| u_{text{ref}}(x,t)| dOmega},]

the place (Omega = [-1,1] occasions [0,T]) is the computational area. To unravel Burgers’ equation, we undertake the usual PINNs strategy and decrease the squared error of the governing equation:

[ underset{boldsymbol theta}{text{argmin }} L(boldsymbol theta), quad L(boldsymbol theta) = frac{1}{2} int_{Omega} Big(frac{partial u}{partial t} + ufrac{partial u}{partial x} – nu frac{partial^2 u}{partial x^2}Big)^2 dOmega.]

The rate discipline is discretized with the assistance of an MLP community (N (x,t;boldsymbol theta)), and the boundary and preliminary circumstances are built-in with a distance function-type strategy [25]:

[ u(x,t;boldsymbol theta) = (1+x)(1-x)Big(t N(x,t; boldsymbol theta) + u_0(x)(1-t/T)Big). ]

On this drawback, we take the viscosity to be (nu=0.01) and the ultimate time to be (T=2). The preliminary situation is given by (u_0(x) = – sin(pi x)), which results in the well-known shock sample at (x=0). We run ADAM optimization for (1.5 occasions 10^{4}) epochs with a studying fee of (1.5 occasions 10^{-3}) to resolve the optimization drawback at every community measurement. By sweeping over (8) community sizes, we once more search for the parameter depend at which the answer departs from the reference answer. Be aware that we confirm our reference answer in opposition to a spectral solver to make sure the accuracy of our implementation. See Determine 3 for the outcomes. All networks with (|boldsymbol theta|geq 150) parameters present roughly equal efficiency. As such, the unique community is overparameterized by an element of (57).

Neohookean hyperelasticity

On this instance, we take into account the nonlinearly elastic deformation of a dice underneath a prescribed displacement. The pressure power density of a 3D hyperelastic stable is given by the compressible Neohookean mannequin [26] as

[PsiBig( mathbf{u}(mathbf{X}) Big) = frac{ell_1}{2}Big( I_1 – 3 Big) – ell_1 ln J + frac{ell_2}{2} Big( ln J Big)^2 ,]

the place (ell_1) and (ell_2) are materials properties which we take as constants. The pressure power makes use of the next definitions:

[ mathbf{F} = mathbf{I} + frac{partial mathbf{u}}{partial mathbf{X}},
I_1 = mathbf{F} : mathbf{F},
J = det(mathbf{F}),]

the place (mathbf{u}) is the displacement discipline, (mathbf{X}) is the place within the reference configuration, and (mathbf F) is the deformation gradient tensor. The displacement discipline is obtained by minimizing the full potential power, given by

[ PiBig( mathbf{u}(mathbf{X}) Big) = int_{Omega} PsiBig( mathbf{u}(mathbf{X}) Big) – mathbf{b} cdot mathbf{u} dOmega – int_{partial Omega} mathbf{t} cdot mathbf{u} dS,]

the place (Omega) is the undeformed configuration of the physique, (mathbf{b}) is a volumetric power, and (mathbf{t}) is an utilized floor traction. Our investigation into the community measurement is impressed by [27], wherein the Deep Ritz technique is used to acquire a minimal of the hyperelastic complete potential power useful. Nevertheless, we choose to make use of the Neohookean mannequin of the pressure power, versus the Lopez-Pamies mannequin they make use of. As within the cited work, we take the undeformed configuration to be the unit dice (Omega=[0,1]^3) and we topic the dice to a uniaxial pressure state. To implement this pressure state, we apply a displacement (U) within the (X_3) path on the highest floor of the dice. Curler helps, which zero solely the (X_3) element of the displacement, are utilized on the underside floor. All different surfaces are traction-free, which is enforced weakly by the chosen power useful. The boundary circumstances are glad robotically by discretizing the displacement as

[ begin{bmatrix}
u_1(mathbf{X}; boldsymbol theta)
u_2(mathbf{X}; boldsymbol theta)
u_3(mathbf{X}; boldsymbol theta)
end{bmatrix} = begin{bmatrix}
X_3 N_1(mathbf{X}; boldsymbol theta)
X_3 N_2(mathbf{X}; boldsymbol theta)
sin(pi X_3) N_3(mathbf{X}; boldsymbol theta) + UX_3
end{bmatrix},]

the place ( N_i) is the (i)-th element of the community output. Within the cited work, a six hidden-layer community of width (40) is used to discretize the three elements of the displacement discipline. This corresponds to (8480) trainable parameters. On condition that the community is a map ( N: mathbb R^3 rightarrow mathbb R^3), a two hidden-layer community of width (M) has (8M+M^2) trainable parameters when no bias is utilized on the output layer. Thus, if we take ( M=88 ), our community has (8448) trainable parameters. We’ll take this community structure to be the big community reference.

In [27], the connection between the conventional element of the primary Piola-Kirchhoff stress tensor (mathbf{P}) within the path of the utilized displacement and the corresponding element of the deformation gradient (mathbf{F}) was computed to confirm their Deep Ritz implementation. Right here, we research the connection between this tensile stress and the utilized displacement (U). The primary Piola-Kirchhoff stress tensor is obtained with the pressure power density as

[ mathbf{P} = frac{ partial Psi}{partial mathbf{F}} = ell_1( mathbf{F} – mathbf{F}^{-T} ) + ell_2 mathbf{F}^{-T}log J.]

Given the unit dice geometry and the uniaxial stress/pressure state, the deformation gradient is given by

[ mathbf{F} = begin{bmatrix}
1 & 0 & 0 0 & 1 & 0 0 & 0 & 1+U
end{bmatrix}.]

With these two equations, we compute the tensile stress (P_{33}) and the pressure power as a perform of the utilized displacement to be

[ begin{aligned} P_{33} = ell_1Big( 1 + U – frac{1}{1+U}Big) + ell_2frac{log(1+U)}{1+U}, Pi = frac{ell_1}{2}(2+(1+U)^2-3) – ell_1 log(1+U) + frac{ell_2}{2}(log(1+U))^2. end{aligned}]

These analytical options can be utilized to confirm our implementation of the hyperelastic mannequin, in addition to to gauge the efficiency of various measurement networks. Utilizing the neural community mannequin, the tensile stress and the pressure power are computed at every utilized displacement with:

[ P_{33} = int_{Omega} ell_1( {mathbf{F}} – {mathbf{F}}^{-T} ) + ell_2 {mathbf{F}}^{-T}log J dOmega, quad Pi = int_{Omega} PsiBig( {mathbf{u}}(mathbf{X})Big) dOmega,]

the place the displacement discipline is constructed from parameters obtained from the Deep Ritz technique. To compute the stress, we common over your entire area, on condition that we count on a relentless stress state. On this instance, the fabric parameters are set at (ell_1=1) and (ell_2=0.25). We iterate over (8) community sizes and take (10) load steps at every measurement to acquire the stress and pressure power as a perform of the utilized displacement. See Determine 4 for the outcomes. All networks precisely reproduce the pressure power and stress loading curves. This consists of even the community of width (2), with solely (20) trainable parameters. Thus, the unique community has (424times) extra parameters than essential to symbolize the outcomes of the tensile take a look at.

A counterexample

Within the fourth and remaining instance, we clear up a regression drawback to point out the failure of small networks to suit high-frequency features. The one-dimensional regression drawback is given by

[ underset{boldsymbol theta}{text{argmin }} L(boldsymbol theta), quad L(boldsymbol theta) = frac{1}{2}int_0^1Big( v(x) – N(x;boldsymbol theta) Big)^2 dx,]

the place ( N) is a two hidden-layer MLP community and (v(x)) is the goal perform. On this instance, we take (v(x)=sin^5(20pi x)). We iterate over (5) totally different community sizes and report the converged loss worth (L) as an error measure. We practice utilizing ADAM optimization for (5 occasions 10^4) epochs and with a studying fee of (5 occasions 10^{-3}). See Determine 5 for the outcomes. In contrast to the earlier three examples, the goal perform is sufficiently complicated that giant networks are required to symbolize it. The converged error decreases monotonically with the parameter depend. We additionally time the coaching process at every community measurement, and be aware the dependence of the run time (in seconds) on the parameter depend. This instance illustrates that representing oscillatory features requires bigger networks, and that the parameter depend drives up the price of coaching.

Conclusion

Whereas tuning hyperparameters governing the loss perform, optimization course of, and activation perform is frequent within the PINNs group, it’s much less frequent to tune the community measurement. With three instance issues taken from the literature, now we have proven that very small networks typically suffice to symbolize PDE options, even when there are discontinuities and/or different localized options. See Desk 1 for a abstract of our outcomes on the opportunity of utilizing small networks. To qualify our findings, we then introduced the case of regression to a high-frequency goal perform, which required a lot of parameters to suit precisely. Thus, our conclusions are as follows: answer fields which don’t oscillate can typically be represented by small networks, even after they comprise localized options comparable to cracks and shocks. As a result of the price of coaching scales with the variety of parameters, smaller networks can expedite coaching for physics-informed issues with non-oscillatory answer fields. In our expertise, such answer fields seem frequently in sensible issues from warmth conduction and static stable mechanics. By shrinking the dimensions of the community, these issues and others symbolize alternatives to render PINN options extra computationally environment friendly, and thus extra aggressive with conventional approaches such because the finite ingredient technique.

References

[1] Justin Sirignano and Konstantinos Spiliopoulos. DGM: A deep studying algorithm for fixing partial differential equations. Journal of Computational Physics, 375:1339–1364, December 2018. arXiv:1708.07469 [q-fin].

[2] M. Raissi, P. Perdikaris, and G.E. Karniadakis. Physics-informed neural networks: A deep studying framework for fixing ahead and inverse issues involving nonlinear partial differential equations. Journal of Computational Physics, 378:686–707, February 2019.

[3] Weinan E and Bing Yu. The Deep Ritz technique: A deep learning-based numerical algorithm for fixing variational issues, September 2017. arXiv:1710.00211 [cs].

[4] Yaohua Zang, Gang Bao, Xiaojing Ye, and Haomin Zhou. Weak Adversarial Networks for Excessive-dimensional Partial Differential Equations. Journal of Computational Physics, 411:109409, June 2020. arXiv:1907.08272 [math].

[5] Reza Khodayi-Mehr and Michael M. Zavlanos. VarNet: Variational Neural Networks for the Answer of Partial Differential Equations, December 2019. arXiv:1912.07443 [cs].

[6] E. Kharazmi, Z. Zhang, and G. E. Karniadakis. Variational Physics-Knowledgeable Neural Networks For Fixing Partial Differential Equations, November 2019. arXiv:1912.00873 [cs].

[7] Xiaowei Jin, Shengze Cai, Hui Li, and George Em Karniadakis. NSFnets (Navier-Stokes Stream nets): Physics-informed neural networks for the incompressible Navier-Stokes equations. Journal of Computational Physics, 426:109951, February 2021.

[8] Shengze Cai, Zhicheng Wang, Sifan Wang, Paris Perdikaris, and George Em Karniadakis. Physics-Knowledgeable Neural Networks for Warmth Switch Issues. Journal of Warmth Switch, 143(060801), April 2021.

[9] Min Liu, Zhiqiang Cai, and Karthik Ramani. Deep Ritz technique with adaptive quadrature for linear elasticity. Pc Strategies in Utilized Mechanics and Engineering, 415:116229, October 2023.

[10] Sifan Wang, Ananyae Kumar Bhartari, Bowen Li, and Paris Perdikaris. Gradient Alignment in Physics-informed Neural Networks: A Second-Order Optimization Perspective, September 2025. arXiv:2502.00604 [cs].

[11] Jorge F. City, Petros Stefanou, and Jose A. Pons. Unveiling the optimization means of physics knowledgeable neural networks: How correct and aggressive can PINNs be? Journal of Computational Physics, 523:113656, February 2025.

[12] Conor Rowan, Kai Hampleman, Kurt Maute, and Alireza Doostan. Boundary situation enforcement with PINNs: a comparative research and verification on 3D geometries, December 2025. arXiv:2512.14941[math].

[13] Nasim Rahaman, Aristide Baratin, Devansh Arpit, Felix Draxler, Min Lin, Fred A. Hamprecht, Yoshua Bengio, and Aaron Courville. On the Spectral Bias of Neural Networks, Might 2019. arXiv:1806.08734 [stat].

[14] Mirco Pezzoli, Fabio Antonacci, and Augusto Sarti. Implicit neural illustration with physics-informed neural networks for the reconstruction of the early a part of room impulse responses. In Proceedings of the tenth Conference of the European Acoustics Affiliation Discussion board Acusticum 2023, pages 2177–2184, January 2024.

[15] Shaghayegh Fazliani, Zachary Frangella, and Madeleine Udell. Enhancing Physics-Knowledgeable Neural Networks By Function Engineering, June 2025. arXiv:2502.07209 [cs].

[16] Duong V. Dung, Nguyen D. Music, Pramudita S. Palar, and Lavi R. Zuhal. On The Alternative of Activation Capabilities in Physics-Knowledgeable Neural Community for Fixing Incompressible Fluid Flows. In AIAA SCITECH 2023 Discussion board. American Institute of Aeronautics and Astronautics. Eprint: https://arc.aiaa.org/doi/pdf/10.2514/6.2023-1803.

[17] Honghui Wang, Lu Lu, Shiji Music, and Gao Huang. Studying Specialised Activation Capabilities for Physics-informed Neural Networks. Communications in Computational Physics, 34(4):869–906, June 2023. arXiv:2308.04073 [cs].

[18] Zhao Zhang, Xia Yan, Piyang Liu, Kai Zhang, Renmin Han, and Sheng Wang. A physics-informed convolutional neural community for the simulation and prediction of two-phase Darcy flows in heterogeneous porous media. Journal of Computational Physics, 477:111919, March 2023.

[19] Pu Ren, Chengping Rao, Yang Liu, Jianxun Wang, and Hao Solar. PhyCRNet: Physics-informed Convolutional-Recurrent Community for Fixing Spatiotemporal PDEs. Pc Strategies in Utilized Mechanics and Engineering, 389:114399, February 2022. arXiv:2106.14103 [cs].

[20] Mikhail Belkin, Daniel Hsu, Siyuan Ma, and Soumik Mandal. Reconciling trendy machine-learning apply and the classical bias–variance trade-off. Proceedings of the Nationwide Academy of Sciences, 116(32):15849–15854, August 2019. Writer: Proceedings of the Nationwide Academy of Sciences.

[21] Arthur Jacot, Franck Gabriel, and Cl´ement Hongler. Neural Tangent Kernel: Convergence and Generalization in Neural Networks, February 2020. arXiv:1806.07572 [cs].

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

[23] M. Manav, R. Molinaro, S. Mishra, and L. De Lorenzis. Part-field modeling of fracture with physics-informed deep studying. Pc Strategies in Utilized Mechanics and Engineering, 429:117104, September 2024.

[24] Xianke Wang, Shichao Yi, Huangliang Gu, Jing Xu, and Wenjie Xu. WF-PINNs: fixing ahead and inverse issues of burgers equation with steep gradients utilizing weak-form physics-informed neural networks. Scientific Studies, 15(1):40555, November 2025. Writer: Nature Publishing Group.

[25] N. Sukumar and Ankit Srivastava. Precise imposition of boundary circumstances with distance features in physics-informed deep neural networks. Pc Strategies in Utilized Mechanics and Engineering, 389:114333, February 2022.

[26] Javier Bonet and Richard D. Wooden. Nonlinear Continuum Mechanics for Finite Component Evaluation. Cambridge College Press, Cambridge, 2 version, 2008.

[27] Diab W. Abueidda, Seid Koric, Rashid Abu Al-Rub, Corey M. Parrott, Kai A. James, and Nahil A. Sobh. A deep studying power technique for hyperelasticity and viscoelasticity. European Journal of Mechanics – A/Solids, 95:104639, September 2022. arXiv:2201.08690 [cs].