REVIEW 4 major objections 5 minor 37 references
A Kolmogorov High Order Deep Neural Network for High Frequency Partial Differential Equations in High Dimensions
T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This paper claims that a neural network built from the Kolmogorov superposition theorem approximates high-frequency, high-dimensional PDE solutions with parameter count polynomial in the dimension, breaking the curse of dimensionality for…
desk verdict A real architectural fix for HOrderDNN's exponential basis count, backed by broad experiments, but the headline CoD claim is conditional on a class K_L that the paper never checks for its own test problems. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the Kolmogorov Superposition Theorem in Lorentz's form, which represents any continuous $d$-variate function as a sum over $2d+1$ univariate inner functions $\varphi_q$ composed with a single outer function $g$: $f(x_1,\ldots,x_d)=\sum_{q=0}^{2d}g(\sum_{i=1}^d\lambda_i\varphi_q(x_i))$. The paper's architecture makes this concrete: the inner functions are learned by a high-order network $h_p$ (polynomial basis of degree $p$, applied to each coordinate with shared weights), and the outer function by a fully connected network $G_{\mathrm{NN}}$, giving the whole ansatz $k_p(x)=G_2\circ G_{\mathrm{NN}}\circ\sigma\circ G_1([h_p(x_1),\ldots,h_p(x_d)]^T)$. The approximation argument then splits the error into two parts: Jackson's theorem bounds the polynomial approximation of the H\"older-$\alpha$ inner functions by $p^{-\alpha}$, and a linear-spline (ReLU) or tanh network bound of the Lipschitz outer function by $n^{-1}$; the product of the two error terms with the $2d+1$ sums produces the $O(d^2(\cdot))$ rate.
What would settle it
Take a ten-dimensional target function whose Kolmogorov outer function is explicitly non-Lipschitz, for example by constructing the KST representation with $g(z)=\sqrt{z}$, and train K-HOrderDNN with increasing $p$; if the error continues to decay at roughly $p^{-0.30}$ or faster, the Lipschitz assumption is not the operative constraint, while if the error stalls or decays much slower, the result is confined to the dense subclass $\mathcal{K}_L$ and does not extend to arbitrary continuous targets.
Extended reading notes
Core claim
K-HOrderDNN replaces the $(p+1)^d$ tensor-product basis functions of HOrderDNN with $d(p+1)$ univariate basis functions arranged through the Kolmogorov–Lorentz representation $f(x)=\sum_{q=0}^{2d}g(\sum_{i=1}^d\lambda_i\varphi_q(x_i))$. The inner functions $\varphi_q$ are approximated by a shared high-order subnetwork $h_p$ applied coordinate-wise, and the outer function $g$ by a fully connected subnetwork; the construction mirrors the proof of Theorem 2, which shows that when $g$ is Lipschitz the error is bounded by a constant times $d^2(1/n+1/p^\alpha)$ with $\alpha=\log_{10}2$. The paper claims this breaks the curse of dimensionality for the dense subset $\mathcal{K}_L$ of continuous functions, and demonstrates numerically that in dimensions 10, 20, and 50 K-HOrderDNN($p>1$) remains trainable and accurate while HOrderDNN($p>1$) becomes intractable.
Load-bearing premise
The theoretical guarantee of breaking the curse of dimensionality holds only for functions whose Kolmogorov outer function is Lipschitz continuous, and the paper does not establish that the high-frequency PDE solutions used in the experiments actually lie in that dense subclass.
Editorial extensions
If this is right
- For target functions whose KST outer function is Lipschitz, the parameter count $2n+d+(p+1)(2d+1)$ is polynomial in the dimension $d$, so the method stays trainable as $d$ grows.
- K-HOrderDNN($p>1$) inherits the high-order basis advantage of HOrderDNN for oscillatory problems but removes the exponential $(p+1)^d$ blow-up, so the order $p$ can be increased in high dimensions without an explosion in parameters.
- In the ten-dimensional Helmholtz benchmark the reported relative L2 error drops from $9.93\times10^{-1}$ (HOrderDNN $p=1$ and PINN) to $4.40\times10^{-3}$ (K-HOrderDNN $p=7$), suggesting orders-of-magnitude gains for high-frequency high-dimensional wave problems.
- With tanh activation the same CoD-avoidance result holds (Theorem 4), making the construction compatible with the second-derivative losses used in PINN training.
Reading between the lines
- The reported experiments show faster-than-$p^{-0.30}$ decay, which the paper attributes to smoothness; measuring the empirical $p$-exponent across problems with controlled smoothness would separate the worst-case guarantee from the smooth case, and is a natural test the paper does not perform.
- Because the KST inner functions are universal (independent of $f$), the shared subnetwork $h_p$ could potentially be pre-trained on one target and fine-tuned for another high-dimensional PDE at nearly constant cost, an extension the paper leaves implicit.
- The paper's comparison with KAN shows comparable accuracy at roughly twelve times faster training, but it does not test highly oscillatory outer functions; a natural extension would be to check whether the architecture retains its advantage when $g$ itself carries the high-frequency content.
- The method inherits the Monte Carlo sampling error of the PINN loss, and the paper suggests insufficient sample size at $d=50$; a hybrid with quadrature-based losses, as used in tensor neural networks, is an untested direction that could tighten the gap between theory and practice at larger dimensions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proposes the Kolmogorov high-order deep neural network (K-HOrderDNN), a variant of HOrderDNN in which the tensor-product high-order layer, whose size is (p+1)^d, is replaced by a KST-inspired construction: a shared high-order univariate subnetwork h_p approximates the K-inner functions, and a fully connected subnetwork approximates the K-outer function. The paper proves approximation rates for a special class of ReLU K-HOrderDNNs (Theorem 2) and a tanh version (Theorem 4), giving parameter counts 2n+d+(p+1)(2d+1) and O(d^2(n^{-1}+p^{-alpha})) error for the Lipschitz-outer class K_L, where alpha=log_10 2; Theorem 3 extends the ReLU result to arbitrary continuous outer functions via the modulus of continuity. Numerical experiments fit high-frequency functions and solve 2D, 10D, 20D, and 50D Poisson and Helmholtz problems, reporting strong empirical gains for K-HOrderDNN(p>1) over PINN and HOrderDNN(p=1), with parameter counts avoiding the exponential (p+1)^d growth. The central limitations are that the CoD-free statement is conditional on membership in K_L, which is not verified for the PDE solutions tested, and the proof of Theorem 4 contains a missing lemma.
Significance. If the theorems were fully established and the K_L-membership gap were closed, the paper would provide a useful explicit mechanism: a dense subclass of C([0,1]^d) approximable with polynomial-in-d parameter count and an algebraic rate O(d^2(n^{-1}+p^{-alpha})). The paper is commendable for stating precise parameter counts, defining the approximation class cleanly, and running a broad numerical comparison against PINN, HOrderDNN, KAN, and TNN, including high-dimensional examples; the factor of roughly 225 improvement in Table 10 is striking. The main value is the architecture and its empirical behavior, while the theoretical CoD claim is narrower than the abstract and conclusion suggest. The paper is also honest in Section 4.2.5 that the faster-than-theory p-convergence is outside the analyzed KST regime. No code is provided, but the experimental setup is described in enough detail to be replicable.
major comments (4)
- [Section 3 (Theorem 2) and Section 4.3-4.4] The CoD-free guarantee is stated for f in K_L, the class whose KST outer function is Lipschitz. The paper does not verify that any of the PDE solutions in Section 4 belongs to K_L. The solutions in Eqs. (4.5)-(4.8) are sums and products of trigonometric polynomials; since the KST inner functions are fixed Holder functions with exponent alpha=log_10 2, the Lipschitz property of the outer function is not automatic and has not been checked. Consequently the polynomial-in-d error bound does not, as stated, apply to the target problems in Tables 6, 7, 9, and 10, and the strong high-dimensional results there remain empirical. Moreover, density of K_L in C([0,1]^d) does not transfer a single uniform CoD-free estimate, because the constant C_g in the bound can grow as the approximating member of K_L approaches the target.
- [Appendix E, Theorem 4] The proof of Theorem 4 asserts, after Eq. (E8), that |phi_q(x)-bar_L_q(x)| <= 3|phi_q(x)-L_q(x)| "by applying the triangle inequality, the definition of bar_L_q, and Lemma E3." No Lemma E3 is stated or proved anywhere in the manuscript. The needed estimate can in fact be obtained from Lemma E2 with w large enough, which suggests that E3 is a typo, but as it stands the proof of the tanh result is incomplete. In addition, the proof of Lemma E2 contains a limit error: the quantity 1 - w tanh(1/w) tends to 0 as w tends to infinity, not as w tends to 0+.
- [Section 3 and Section 4.2.5] The parameter counts in Theorems 2 and 4 are for the special classes K^sigma_{p,n} (one-hidden-layer ReLU spline outer network) and K^{sigma1}_{p,n} (outer tanh network with layer widths at most N-1 and 6N), not for the implemented K-HOrderDNN, whose g-subnetwork has depth gd at least 2 and width gw. Section 4.2.5 identifies n with gw, but the theoretical n is the number of knots in a linear spline, not the width parameter of a deep fully connected GNN. For example, for the d=2, p=9, hw=45, gw=90 configuration reported in Table 4, Theorem 4's count is (2d+6N+14)(N-1)+13+(p+1)(2d+1), which equals 49,725 if N=90 and 1,575 if N=15; neither matches the reported 14,136 parameters. Thus the proved parameter-efficiency rates are not directly tied to the network used in the experiments.
- [Section 4.2.5, Tables 8 and 10] The experiments converge in p much faster than the theoretical rate p^{-alpha} with alpha=log_10 2 approx 0.301. In Table 10, K-HOrderDNN(p=7) reaches 4.40E-03, whereas the worst-case p-term in the bound at p=7 is only about 0.55 times its value at p=1, so the theoretical rate cannot explain even one order of magnitude of the observed improvement. Section 4.2.5 explicitly attributes the faster decay to the potentially higher smoothness of the K-inner function, i.e., to behavior outside the analyzed KST construction. Statements such as "This observation aligns well with Theorem 2" in Section 4.2.1 therefore overstate the theory's explanatory power for the reported accuracies; either add a theorem for smoother inner functions or confine the theory claim to the n-rate.
minor comments (5)
- [Abstract and Section 2.3] The abstract says K-HOrderDNN reduces the number of introduced basis functions to d(p+1), while the definition and Theorem 2 use a shared set of p+1 univariate basis functions with (p+1)(2d+1) trainable coefficients. These counts should be reconciled to avoid confusion.
- [Section 3, definition of K^{sigma1}_{p,n}] The tanh class is denoted K^{sigma1}_{p,n}, but its definition and Theorem 4 use N for the outer network size and never define n; the notation should be unified, for example K^{sigma1}_{p,N}.
- [Section 4.1] The sentence "the depths of PINN and HOrderDNN match the depth of K-HOrderDNN" is ambiguous because K-HOrderDNN has two subnetworks h_p and g with separate depths hd and gd; the rule L = hd + gd + 1 should be stated explicitly before it is used in the experiments.
- [Table 3] Several entries in Table 3 are backslashes without an explanatory note; please add a note that HOrderDNN(p>1) was not executed for d=20 and d=50 because the parameter count is prohibitive.
- [Section 2.3, Eq. (2.9)] The roles of G1 and G2 in Eq. (2.9) are not specified in terms of input and output dimensions; clarifying these layers would help connect the definition to the theoretical class K^sigma_{p,n} in Section 3.
Circularity Check
No significant circularity: the CoD bound is derived from external KST, Jackson, and spline approximation lemmas; self-citations are architectural precedents, not load-bearing.
full rationale
The central derivation chain is self-contained. Theorem 2 combines the external Kolmogorov Superposition Theorem (Theorem 1), Jackson's polynomial approximation bound (Lemma 1), and a standard linear-spline estimate (Lemma 2, from Lai and Shen [26]) to bound the error of the special network class Kσ_{p,n}, whose definition (Eqs. (3.1)-(3.3)) is exactly the KST formula with polynomial-clipped inner approximants and a ReLU-spline outer approximant. The hypothesis class K_L is a regularity condition on the KST outer function, not on the network, and the stated error O(d^2(n^{-1}+p^{-α})) with parameter count 2n+d+(p+1)(2d+1) is proved rather than assumed. The citations to the authors' earlier HOrderDNN work ([15], [16]) supply the high-order transformation layer as an architectural ingredient, but the CoD argument does not reduce to those papers. The numerical results in Section 4 are benchmark comparisons against PINN and HOrderDNN, not fitted parameters relabeled as predictions. Section 4.2.5 explicitly admits that the observed p-convergence is faster than the theoretical O(p^{-α}) and attributes this to possibly higher smoothness of the K-inner functions; that is an acknowledged scope limitation, as is the unverified membership of the Section 4 PDE solutions in K_L. These are correctness or applicability concerns, not circularity.
Assumptions & free parameters
free parameters (4)
- p (polynomial degree/order) =
1,3,5,7,9,40 in experiments
- n (spline knots, approximated by gw) =
gw = 45 to 315 in experiments
- lambda_i (KST arithmetic coefficients) =
learned during training
- All trainable DNN weights in hp and GNN =
not reported (Xavier init, Adam)
assumptions (5)
- standard math Kolmogorov Superposition Theorem (Theorem 1) holds for every continuous f on [0,1]^d with inner functions in Lip(log_10 2).
- standard math The class K_L of functions with Lipschitz KST outer function is dense in C([0,1]^d).
- domain assumption PDE solutions in the experiments belong to K_L or have controlled outer-function modulus of continuity.
- domain assumption Training by Adam finds a network close to the best approximation in K_sigma_p,n.
- ad hoc to paper The unstated bound |phi_q - bar_L_q| <= 3|phi_q - L_q| (missing Lemma E3) is true.
Cite this review
Pith. "Pith review of A Kolmogorov High Order Deep Neural Network for High Frequency Partial Differential Equations in High Dimensions." pith.science (2026). https://pith.science/paper/6BSIV4WD
@misc{pith2026250201938,
author = {Pith},
title = {Pith review of: A Kolmogorov High Order Deep Neural Network for High Frequency Partial Differential Equations in High Dimensions},
year = {2026},
howpublished = {\url{https://pith.science/paper/6BSIV4WD}},
note = {Machine review of arXiv:2502.01938}
}
abstract
This paper proposes a Kolmogorov high order deep neural network (K-HOrderDNN) for solving high-dimensional partial differential equations (PDEs), which improves the high order deep neural networks (HOrderDNNs). HOrderDNNs have been demonstrated to outperform conventional DNNs for high frequency problems by introducing a nonlinear transformation layer consisting of $(p+1)^d$ basis functions. However, the number of basis functions grows exponentially with the dimension $d$, which results in the curse of dimensionality (CoD). Inspired by the Kolmogorov superposition theorem (KST), which expresses a multivariate function as superpositions of univariate functions and addition, K-HOrderDNN utilizes a HOrderDNN to efficiently approximate univariate inner functions instead of directly approximating the multivariate function, reducing the number of introduced basis functions to $d(p+1)$. We theoretically demonstrate that CoD is mitigated when target functions belong to a dense subset of continuous multivariate functions. Extensive numerical experiments show that: for high-dimensional problems ($d$=10, 20, 50) where HOrderDNNs($p>1$) are intractable, K-HOrderDNNs($p>1$) exhibit remarkable performance. Specifically, when $d=10$, K-HOrderDNN($p=7$) achieves an error of 4.40E-03, two orders of magnitude lower than that of HOrderDNN($p=1$) (see Table 10); for high frequency problems, K-HOrderDNNs($p>1$) can achieve higher accuracy with fewer parameters and faster convergence rates compared to HOrderDNNs (see Table 8).
Figures
Figures from the paper (14 more)
Reference graph
Works this paper leans on
-
[1]
J. Han, A. Jentzen, W. E, Solving high-dimensional partial differential equations using deep learning, Proceedings of the National Academy of Sciences 115 (34) (2018) 8505–8510
2018
-
[2]
Yu, et al., The deep ritz method: a deep learning-based numerical algorithm for solving variational problems, Communications in Mathematics and Statistics 6 (1) (2018) 1–12
B. Yu, et al., The deep ritz method: a deep learning-based numerical algorithm for solving variational problems, Communications in Mathematics and Statistics 6 (1) (2018) 1–12
2018
- [3]
-
[4]
Y. Qian, Y. Zhang, Y. Huang, S. Dong, Physics-informed neural networks for approximating dynamic (hyperbolic) pdes of second order in time: Error analysis and algorithms, Journal of Computational Physics 495 (2023) 112527. 32
work page 2023
-
[5]
Sirignano, K
J. Sirignano, K. Spiliopoulos, Dgm: A deep learning algorithm for solving partial differential equations, Journal of computational physics 375 (2018) 1339–1364
2018
-
[6]
Y. Zang, G. Bao, X. Ye, H. Zhou, Weak adversarial networks for high-dimensional partial differential equations, Journal of Computational Physics 411 (2020) 109409
2020
-
[7]
Z. J. Xu, Understanding training and generalization in deep learning by fourier analysis, arXiv preprint arXiv:1808.04295 (2018)
arXiv 2018
-
[8]
Z.-Q. J. Xu, Y. Zhang, T. Luo, Y. Xiao, Z. Ma, Frequency principle: Fourier analysis sheds light on deep neural networks, Communications in Computational Physics 28 (5) (2020) 1746–1767. doi:https://doi.org/10.4208/cicp.OA-2020-0085
Show all 37 references
-
[9]
W. Cai, X. Li, L. Liu, A phase shift deep neural network for high frequency approximation and wave problems, SIAM Journal on Scientific Computing 42 (5) (2020) A3285–A3312
2020
-
[10]
Cai, Z.-Q
W. Cai, Z.-Q. J. Xu, Multi-scale deep neural networks for solving high dimensional pdes, arXiv preprint arXiv:1910.11710 (2019)
2019 arXiv
-
[11]
Z. Liu, W. Cai, Z.-Q. J. Xu, Multi-scale deep neural network (mscalednn) for solving poisson- boltzmann equation in complex domains, arXiv preprint arXiv:2007.11207 (2020)
2020 arXiv
-
[12]
Li, Z.-Q
X.-A. Li, Z.-Q. J. Xu, L. Zhang, A dnn-based algorithm for multi-scale elliptic problems, arXiv e-prints (2020) arXiv–2009
2020
-
[14]
Bai, G.-R
J. Bai, G.-R. Liu, A. Gupta, L. Alzubaidi, X.-Q. Feng, Y. Gu, Physics-informed radial basis network (pirbn): A local approximating neural network for solving nonlinear partial dif- ferential equations, Computer Methods in Applied Mechanics and Engineering 415 (2023) 116290. do...
2023
-
[15]
Chang, K
Z. Chang, K. Li, X. Zou, X. Xiang, High order deep neural network for solving high fre- quency partial differential equations, Communications in Computational Physics 31 (2) (2022) 370–397
2022
-
[16]
Chang, K
Z. Chang, K. Li, X. Xiang, High order deep domain decomposition method for solving high frequency interface problems, Advances in Applied Mathematics and Mechanics 15 (6) (2023) 1602–1630. doi:10.4208/aamm.OA-2022-0006
2023 doi
-
[17]
A. N. Kolmogorov, On the representation of continuous functions of many variables by superposition of continuous functions of one variable and addition, Dokl.akad.nauk Sssr 114 (5) (1957) 953–956
1957
-
[18]
Lorentz, Metric entropy, widths, and superpositions of functions, The American Mathe- matical Monthly 69 (6) (1962) 469–485
G. Lorentz, Metric entropy, widths, and superpositions of functions, The American Mathe- matical Monthly 69 (6) (1962) 469–485
1962
-
[19]
Lorentz, Approximation of functions.-holt, rinehart and wilson, Inc., New York (1966)
G. Lorentz, Approximation of functions.-holt, rinehart and wilson, Inc., New York (1966)
1966
-
[20]
Doss, On the representation of the continuous functions of two variables by mean of addition and continuous functions of one variable, Colloquium Mathematicum 10 (2) (2017) 249–259
R. Doss, On the representation of the continuous functions of two variables by mean of addition and continuous functions of one variable, Colloquium Mathematicum 10 (2) (2017) 249–259. doi:10.4064/cm-10-2-249-259
2017 doi
-
[21]
S. A. Morris, Hilbert 13: Are there any genuine continuous multivariate real-valued func- tions?, Bulletin of the American Mathematical Society (2020) 107–118
2020
-
[22]
D. A. Sprecher, A numerical implementation of kolmogorov’s superpositions ii, Neural Net- works 10 (3) (1997) 447–457. doi:10.1016/s0893-6080(96)00073-1
1997 doi
-
[23]
Braun, M
J. Braun, M. Griebel, On a constructive proof of kolmogorov’s superposition theorem, Con- structive Approximation 30 (3) (2009) 653–675. doi:10.1007/s00365-009-9054-2
2009 doi
-
[24]
D. A. Sprecher, S. Draghici, Space-filling curves and kolmogorov superposition-based neural networks, Neural Networks 15 (1) (2002) 57–67. doi:10.1016/s0893-6080(01)00107-1. 33
2002 doi
-
[25]
Igelnik, N
B. Igelnik, N. Parikh, Kolmogorov’s spline network, IEEE transactions on neural networks 14 (4) (2003) 725–733
2003
-
[26]
M.-J. Lai, Z. Shen, The kolmogorov superposition theorem can break the curse of dimen- sionality when approximating high dimensional functions (2022). arXiv:2112.09963
2022
-
[27]
K ˚urkov´a, Kolmogorov’s theorem is relevant, Neural Computation (1991) 617–622 doi: 10.1162/neco.1991.3.4.617
V . K ˚urkov´a, Kolmogorov’s theorem is relevant, Neural Computation (1991) 617–622 doi: 10.1162/neco.1991.3.4.617
1991 doi
-
[28]
K ˚urkov´a, Kolmogorov’s theorem and multilayer neural networks, Neural Networks (1992) 501–506doi:10.1016/0893-6080(92)90012-8
V . K ˚urkov´a, Kolmogorov’s theorem and multilayer neural networks, Neural Networks (1992) 501–506doi:10.1016/0893-6080(92)90012-8
1992 doi
-
[29]
Schmidt-Hieber, The kolmogorov-arnold representation theorem revisited, Neural net- works 137 (2021) 119–126
J. Schmidt-Hieber, The kolmogorov-arnold representation theorem revisited, Neural net- works 137 (2021) 119–126
2021
-
[30]
Montanelli, H
H. Montanelli, H. Yang, Error bounds for deep relu networks using the kolmogorov–arnold superposition theorem, Neural Networks 129 (2020) 1–6
2020
-
[31]
Wang, Numerical Approximation, 2nd Edition, Beijing Higher Education Press, Beijing, 2012
R. Wang, Numerical Approximation, 2nd Edition, Beijing Higher Education Press, Beijing, 2012
2012
-
[32]
Maczuga, M
P . Maczuga, M. Paszy ´nski, Influence of activation functions on the convergence of physics- informed neural networks for 1d wave equation, in: International Conference on Computa- tional Science, Springer, 2023, pp. 74–88
2023
-
[33]
S. Wang, Y. Teng, P . Perdikaris, Understanding and mitigating gradient flow pathologies in physics-informed neural networks, SIAM Journal on Scientific Computing 43 (5) (2021) A3055–A3081
2021
-
[34]
S. Zeng, Z. Zhang, Q. Zou, Adaptive deep neural networks methods for high-dimensional partial differential equations, Journal of Computational Physics 463 (2022) 111232
2022
-
[35]
Y. Wang, Z. Lin, Y. Liao, H. Liu, H. Xie, Solving high-dimensional partial differential equa- tions using tensor neural network and a posteriori error estimators, Journal of Scientific Computing 101 (3) (2024) 1–29
2024
-
[36]
Z. Liu, Y. Wang, S. Vaidya, F. Ruehle, J. Halverson, M. Soljaˇci´c, T. Y. Hou, M. Tegmark, Kan: Kolmogorov-arnold networks, arXiv preprint arXiv:2404.19756 (2024)
2024 arXiv
-
[37]
Y. Li, Z. Lin, Y. Wang, H. Xie, Tensor neural network interpolation and its applications, arXiv preprint arXiv:2404.07805 (2024)
2024 arXiv
-
[38]
De Ryck, S
T. De Ryck, S. Lanthaler, S. Mishra, On the approximation of functions by tanh neural net- works, Neural Networks 143 (2021) 732–750. 34 Appendix A. Proof of Theorem 3 Proof. Let g be the K-outer function associated with f , whose smoothness is characterized by the standard mo...
2021
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.