REVIEW 2 major objections 4 minor 50 references
Near-optimal Delta-convex Estimation of Lipschitz Functions
T0 review · 2 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read A tractable estimator matches the optimal error rate for Lipschitz regression, adapting to the intrinsic dimension of the data.
desk verdict A genuinely new convex-optimization estimator for Lipschitz regression with a plausible near-minimax proof; the main vulnerability is a load-bearing external covering lemma whose local justification in the paper is suspect. 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 object is the feature map φ▷(x, x̂) = (x−x̂, ‖x−x̂‖▷), which turns each affine piece into a 'cone' term; the function class F▷(X̂) is the max over such terms centered at the cover points. The key identity is that the classical Lipschitz extension f̂(x)=max_{x̂∈X_ε}(f(x̂)−λ‖x−x̂‖) approximates any λ-Lipschitz f to within O(λε) on an ε-cover and is itself Lipschitz with constant O(λ). DCF then solves a convex second-order-cone program over these parameters, with constraints enforced only at the cover centers, and a refinement step that prunes unused terms and re-centers the intercept.
What would settle it
Generate n i.i.d. samples from a subgaussian distribution whose support is a 1-dimensional manifold embedded in R^10, run the adaptive farthest-point clustering as defined in Algorithm 1, and check whether the returned number of centers K is O(n^{1/3}) and the covering radius is O(√(K/n)). A single example where K grows faster than n^{d_*/(2+d_*)} or the covering radius exceeds the stated bound would refute the cited lemma and collapse Theorem 1; alternatively, a constructed distribution with d_*≪d where the estimator's squared error decays slower than n^{-2/(2+d_*)} times a polylog factor wou
Extended reading notes
Core claim
The central claim is that functions of the form max over centers of (affine term plus norm term) — a subclass of delta-convex functions — are universal approximators of Lipschitz functions at uniform error O(λε) while preserving the Lipschitz constant up to a constant factor. Combined with an adaptive farthest-point clustering that produces K ≈ n^{d_*/(2+d_*)} centers covering the data at resolution ε ≈ n^{-1/(2+d_*)}, this yields an empirical-risk-minimization estimator whose squared error is within logarithmic factors of the minimax lower bound. The proof decomposes the error into approximation and estimation terms, controls the former with the covering argument and the latter with metric-
Load-bearing premise
The rate proof relies on a previously published covering lemma, cited without proof here, that the adaptive farthest-point clustering always returns K = O(n^{d_*/(2+d_*)}) centers whose covering radius is O(R_X√(K/n)); if this data-dependent covering guarantee fails on some subgaussian distribution, the n^{-2/(2+d_*)} rate does not follow.
Editorial extensions
If this is right
- If Theorem 1 holds, Lipschitz regression can be solved at the minimax rate by a polynomial-time algorithm that adapts to low intrinsic dimension without any external model-selection procedure.
- The near-minimax guarantee also covers max-min-affine variants and the symmetric difference representation F^∆, so practitioners can choose richer function classes without losing the rate.
- The same machinery gives a convex regression estimator with a PAC guarantee for all n and no knowledge of the Lipschitz constant, improving on earlier results that hold only for large n or high dimension.
- The O((1+d I{▷≠2})) factor identifies Euclidean norm features (▷=2) as the most parameter-efficient choice, with other norms costing a dimension factor.
Reading between the lines
- The analysis suggests that the covering-resolution tradeoff is the real obstacle to adaptivity: any clustering method that returns O(n^{d_*/(2+d_*)}) centers with covering radius O(n^{-1/(2+d_*)}) would slot into the same rate proof, so other data-dependent quantizations could plausibly match DCF.
- The proof's exclusion of quadratic features (‖x−x̂‖²) points to an open avenue: a locally Lipschitz regularizer that keeps the quadratic coefficient independent of ε would likely extend near-minimax rates to smooth regression functions.
- The documented sensitivity of empirical performance to the slope penalty θ₂ suggests that data-dependent selection of this regularization parameter is a natural next target; the theory only requires θ₂ to lie in a range, not a specific value.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a polynomial-time estimator (delta-convex fitting, DCF) for nonparametric regression of an unknown Lipschitz function in the random-design setting, under squared loss and subgaussian covariate/noise assumptions. The estimator combines adaptive farthest-point clustering (AFPC) to select center points, a convex optimization problem over a class of piecewise affine/"delta-convex" functions, a regularization that avoids knowing the true Lipschitz constant, and an optional local refinement. The main theoretical claim (Theorem 1) is a PAC bound of order O(d(1+d I{▷≠2}) n^{-2/(2+d_*)} β), with β polylogarithmic, matching Stone's minimax lower bound up to logarithmic factors in terms of the intrinsic dimension d_*. The paper also gives approximation results for max-min-affine and other delta-convex classes and extends the framework to convex regression. Experiments compare DCF with k-NN, Nadaraya-Watson, random forests, and XGBoost on three public datasets.
Significance. If the main theorem holds, this is a genuinely valuable contribution: it is, to my knowledge, the first tractable estimator for general Lipschitz regression that is proven to achieve the adaptive near-minimax rate in the random-design setting without knowledge of the Lipschitz constant, while also admitting a convex initialization. The approximation construction based on McShane's extension and the use of AFPC to adapt to intrinsic dimension are elegant, and the paper explicitly provides a PAC bound for all n rather than only asymptotic convergence. The experiments are a useful complement, and the code is made available. However, two load-bearing points in the proof currently prevent me from accepting the paper as written: the adaptive covering lemma is cited without proof and with an incorrect local justification, and an algebraic simplification in the concentration step appears invalid and would, as written, break the n^{-2/(2+d_*)} rate. Both appear repairable, but they need to be fixed before the central claim is established.
major comments (2)
- [§2.1 and §4.3 (Lemma 8)] The paper states that the AFPC stopping rule implies 'K−1 < k̂(X_K) ≤ K' and uses this to justify ε² ≈ R² K/n. This inequality is not guaranteed: adding the final farthest point can sharply reduce ε_n(X_K) (e.g., two tight clusters), making k̂(X_K) much smaller than K−1. The useful upper bound ε = O(R√(K/n)) does follow directly from K ≥ k̂(X_K), and the complexity bound K ≤ k* requires K−1 < k̂(X_{K−1}) together with a covering/doubling argument. Neither is supplied. Since Theorem 1's rate enters exactly through Lemma 8, the manuscript must either prove Lemma 8 (especially the d* part) or provide a correct derivation; the current local justification is incorrect.
- [§4.4.2, Lemma 16, Eq. (20)] With δ_n = r_σ√(dK/n), the term r_σ δ_n in the concentration bound equals r_σ²√(dK/n) up to constants. The displayed simplification in Eq. (20) claims this is O(dKσ² ln(β2/γ)/n), which is false when dK/n < 1 — exactly the regime relevant for the near-minimax rate, where dK/n → 0. If this term is retained, the bound on E_approx in Eq. (23) would contain a √(K/n) term rather than K/n, destroying the n^{-2/(2+d_*)} rate. The proof can likely be repaired by choosing δ_n = r_σ dK/n (or a similar smaller scale), since δ appears only logarithmically in Lemma 15; but as written the central rate is not established.
minor comments (4)
- [Throughout §4] The prose consistently refers to lemmas as 'Theorem' (e.g., 'Theorem 4', 'Theorem 8', 'Theorem 16'). The cross-references should be corrected to Lemma numbers.
- [§2.1] The sentence 'The stopping condition of AFPC ensures that K−1 < k̂(X_K) ≤ K' is false as noted in Major 1. It should be rephrased, for example as 'K ≥ k̂(X_K) and K−1 < k̂(X_{K−1})'.
- [Proof of Theorem 1, after Eq. (33)] The simplification '(r_ρ β_0)/(ηδ) = O(n²√d n)' is ambiguous and appears to be a typo; it should be O(n²√d ln n) (or, if using ln n ≤ n, written with the explicit substitution).
- [Various] Typos and stylistic issues: 'Cauchy-Schwartz' should be 'Cauchy-Schwarz'; 'strenghten' → 'strengthen'; 'straightforwadly' → 'straightforwardly'; 'implmentation' → 'implementation'; 'samples sizes' → 'sample sizes'.
Circularity Check
No significant circularity: the near-minimax rate follows from a genuine upper-bound proof with independent geometric and concentration lemmas.
full rationale
Theorem 1's rate is derived by bounding approximation error via Theorem 2 (McShane extension over an epsilon-cover) and estimation error via ERM/concentration inequalities, then plugging in the AFPC covering guarantee (Lemma 8). Lemma 8 is cited from the author's prior work, but it is a parameter-free clustering statement about AFPC's output (K and epsilon), not a restatement of the regression rate, and it is externally checkable. The statistical analysis—error decomposition (11), the 'basic inequality' (16), and concentration lemmas (11, 17)—is carried out in the paper. The regularization (7) is chosen as a function of K and n rather than fitted to the target loss, and the proof does not tune parameters to match the minimax rate. The only notable burden is reliance on self-cited AFPC and prior concentration results; these are independent published results and do not make the derivation circular. A separate correctness concern: Section 2.1's sentence 'The stopping condition of AFPC (Algorithm 1) ensures that K−1 < khat(X_K) ≤ K' is not an obvious consequence of the loop condition, since the loop condition is evaluated before the final insertion and the covering radius can drop sharply; this affects the informal justification of Lemma 8 but is a proof-rigor issue, not circularity.
Assumptions & free parameters
free parameters (4)
- θ_0 =
O((R_Yn/max{1,R_Xn}) ln n); experiments use (R_Yn/R_Xn) ln n
- θ_1 =
Θ(max{1,R_Xn^2}(dK/n))
- θ_2 =
0≤θ_2≤θ_1/K; experiments use R_Xn^2/n or (R_Xn/n)^2
- θ_3 =
1≤θ_3=O(ln n); experiments use ln n
assumptions (4)
- standard math McShane extension theorem: every λ-Lipschitz f on X extends to R^d via sup_{\hat x∈X} f(\hat x)-λ\|x-\hat x\|.
- domain assumption Doubling-dimension covering bound: N_{\|·\|}(X_n, ε) ≤ max{1,(4R_Xn/ε)^{d_*}} where d_* is the intrinsic dimension.
- domain assumption Subgaussian model: E[e^{\|X-E[X]\|^2/ρ^2}] ≤ 2 and E[e^{(f^*(X)-Y)^2/σ^2}|X] ≤ 2 a.s.
- domain assumption AFPC guarantee (Balázs 2022, Lemma 4.2): AFPC output \hat X has K=O(n^{d_*/(2+d_*)}) a.s. and ε_n(\hat X)=O(R_Xn√(K/n)).
Cite this review
Pith. "Pith review of Near-optimal Delta-convex Estimation of Lipschitz Functions." pith.science (2026). https://pith.science/paper/RITXYIKZ
@misc{pith2026251115615,
author = {Pith},
title = {Pith review of: Near-optimal Delta-convex Estimation of Lipschitz Functions},
year = {2026},
howpublished = {\url{https://pith.science/paper/RITXYIKZ}},
note = {Machine review of arXiv:2511.15615}
}
read the original abstract
This paper presents a tractable algorithm for estimating an unknown Lipschitz function from noisy observations and establishes an upper bound on its convergence rate. The approach extends max-affine methods from convex shape-restricted regression to the more general Lipschitz setting. A key component is a nonlinear feature expansion that maps max-affine functions into a subclass of delta-convex functions, which act as universal approximators of Lipschitz functions while preserving their Lipschitz constants. Leveraging this property, the estimator attains the minimax convergence rate (up to logarithmic factors) with respect to the intrinsic dimension of the data under squared loss and subgaussian distributions in the random design setting. The algorithm integrates adaptive partitioning to capture intrinsic dimension, a penalty-based regularization mechanism that removes the need to know the true Lipschitz constant, and a two-stage optimization procedure combining a convex initialization with local refinement. The framework is also straightforward to adapt to convex shape-restricted regression. Experiments demonstrate competitive performance relative to other theoretically justified methods, including nearest-neighbor and kernel-based regressors.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
A parallelizable dual smoothing method for large scale convex regression problems
Necdet Serhat Aybat and Zi Wang. A parallelizable dual smoothing method for large scale convex regression problems. arXiv preprint arXiv:1608.02227, 2016
arXiv 2016
-
[2]
An algorithm for the estimation of a regression function by continuous piecewise linear functions
Adil Bagirov, Conny Clausen, and Michael Kohler. An algorithm for the estimation of a regression function by continuous piecewise linear functions. Computational Optimization and Applications, 45 0 (1): 0 159--179, 2010
2010
-
[3]
Bagirov, Sona Taheri, Napsu Karmitsa, Nargiz Sultanova, and Soodabeh Asadi
Adil M. Bagirov, Sona Taheri, Napsu Karmitsa, Nargiz Sultanova, and Soodabeh Asadi. Robust piecewise linear L1 -regression via nonsmooth DC optimization. Optimization Methods and Software, 37 0 (4): 0 1289--1309, 2022
2022
-
[4]
Convex Regression: Theory, Practice, and Applications
G\'abor Bal\'azs. Convex Regression: Theory, Practice, and Applications. PhD thesis, University of Alberta, 2016
2016
-
[5]
Adaptively partitioning max-affine estimators for convex regression
G\'abor Bal\'azs. Adaptively partitioning max-affine estimators for convex regression. In Gustau Camps-Valls, Francisco J. R. Ruiz, and Isabel Valera, editors, Proceedings of The 25th International Conference on Artificial Intelligence and Statistics, volume 151 of Proceedings of Machine Learning Research, pages 860--874. PMLR, 2022
2022
-
[6]
Near-optimal max-affine estimators for convex regression
G\'abor Bal\'azs, Andr\'as György, and Cs aba Szepesv\'ari. Near-optimal max-affine estimators for convex regression . In Guy Lebanon and S. V. N. Vishwanathan, editors, Proceedings of the Eighteenth International Conference on Artificial Intelligence and Statistics, volume 38 of Proceedings of Machine Learning Research, pages 56--64, San Diego, Californi...
2015
-
[7]
Chaining bounds for empirical risk minimization
G\'abor Bal\'azs, Andr\'as Gy\"orgy, and Cs aba Szepesv\'ari. Chaining bounds for empirical risk minimization. arXiv preprint arXiv:1609.01872v1, 2016
arXiv 2016
-
[8]
Bartlett and Shahar Mendelson
Peter L. Bartlett and Shahar Mendelson. Empirical minimization. Probability Theory and Related Fields, 135 0 (3): 0 311--334, 2006
2006
Show all 50 references
-
[9]
Bartlett, St\'ephane Boucheron, and G\'abor Lugosi
Peter L. Bartlett, St\'ephane Boucheron, and G\'abor Lugosi. Model selection and error estimation. Machine Learning, 48 0 (1): 0 85--113, 2002
2002
-
[10]
Multivariate distributionally robust convex regression under absolute error loss
Jose Blanchet, Peter W Glynn, Jun Yan, and Zhengqing Zhou. Multivariate distributionally robust convex regression under absolute error loss. In H. Wallach, H. Larochelle, A. Beygelzimer, F. d Alch\' e -Buc, E. Fox, and R. Garnett, editors, Advances in Neural Information Proces...
2019
-
[11]
Concentration Inequalities: A nonasymptotic theory of independence
St\'ephane Boucheron, G\'abor Lugosi, and Pascal Massart. Concentration Inequalities: A nonasymptotic theory of independence. Oxford University Press, 2013
2013
-
[12]
Convex Optimization
Stephen Boyd and Lieven Vandenberghe. Convex Optimization. Cambridge University Press, 2004
2004
-
[13]
Random forests
Leo Breiman. Random forests. Machine Learning, 45 0 (1): 0 5--32, 2001
2001
-
[14]
XGB oost: A S calable T ree B oosting S ystem
Tianqi Chen and Carlos Guestrin. XGB oost: A S calable T ree B oosting S ystem. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, KDD '16, page 785–794. Association for Computing Machinery, 2016. doi:10.1145/2939672.2939785
2016
-
[15]
Subgradient regularized multivariate convex regression at scale
Wenyu Chen and Rahul Mazumder. Subgradient regularized multivariate convex regression at scale. SIAM Journal on Optimization, 34 0 (3): 0 2350--2377, 2024
2024
-
[16]
Random projection trees and low dimensional manifolds
Sanjoy Dasgupta and Yoav Freund. Random projection trees and low dimensional manifolds. In Proceedings of the Fortieth Annual ACM Symposium on Theory of Computing, STOC '08, page 537–546, New York, NY, USA, 2008. Association for Computing Machinery
2008
-
[17]
DeVore, Ralph Howard, and Charles Micchelli
Ronald A. DeVore, Ralph Howard, and Charles Micchelli. Optimal nonlinear approximation. Manuscripta Mathematica, 63 0 (4): 0 469--478, 1989
1989
-
[18]
Gonzalez
Teofilo F. Gonzalez. Clustering to minimize the maximum intercluster distance. Theoretical Computer Science, 38: 0 293--306, 1985
1985
-
[19]
M axout N etworks
Ian Goodfellow, David Warde-Farley, Mehdi Mirza, Aaron Courville, and Yoshua Bengio. M axout N etworks. In Sanjoy Dasgupta and David McAllester, editors, Proceedings of the 30th International Conference on Machine Learning, volume 28 of Proceedings of Machine Learning Research...
2013
-
[20]
Goulart and Yuwen Chen
Paul J. Goulart and Yuwen Chen. Clarabel: An interior-point solver for conic programs with quadratic objectives. arXiv preprint arXiv:2405.12762, 2024
2024 arXiv
-
[21]
Gupta, R
A. Gupta, R. Krauthgamer, and J.R. Lee. Bounded geometries, fractals, and low-distortion embeddings. In 44th Annual IEEE Symposium on Foundations of Computer Science, 2003. Proceedings., pages 534--543, 2003
2003
-
[22]
A distribution-free theory of nonparametric regression
L\'aszl\'o Gy\"orfi, Michael Kohler, Adam Krzyzak, and Harro Walk. A distribution-free theory of nonparametric regression. Springer, 2002
2002
-
[23]
Qiyang Han and Jon A. Wellner. Multivariate convex regression: Global risk bounds and adaptation. arXiv preprint arXiv:1601.06844v1, 2016
2016 arXiv
-
[24]
On functions representable as a difference of convex functions
Philip Hartman. On functions representable as a difference of convex functions. Pacific Journal of Mathematics, 9 0 (3): 0 707--713, 1959
1959
-
[25]
Extension of L ipschitz functions
J.-B Hiriart-Urruty. Extension of L ipschitz functions. Journal of Mathematical Analysis and Applications, 77 0 (2): 0 539--554, 1980
1980
-
[26]
Generalized differentiability, duality and optimization for problems dealing with differences of convex functions
Jean-Baptiste Hiriart-Urruty. Generalized differentiability, duality and optimization for problems dealing with differences of convex functions. Convexity and Duality in Optimization, pages 37--70, 1985
1985
-
[27]
Hochbaum and David B
Dorit S. Hochbaum and David B. Shmoys. A best possible heuristic for the k -center problem. Mathematics of Operations Research, 10 0 (2): 0 180--184, 1985
1985
-
[28]
The curse of dimension in nonparametric regression
Samory Kpotufe. The curse of dimension in nonparametric regression. PhD thesis, University of California, 2010
2010
-
[29]
A tree-based regressor that adapts to intrinsic dimension
Samory Kpotufe and Sanjoy Dasgupta. A tree-based regressor that adapts to intrinsic dimension. Journal of Computer and System Sciences, 78 0 (5): 0 1496--1515, 2011
2011
-
[30]
Kulkarni and S.E
S.R. Kulkarni and S.E. Posner. Rates of convergence of nearest neighbor estimation under arbitrary sampling. IEEE Transactions on Information Theory, 41 0 (4): 0 1028--1039, 1995
1995
-
[31]
Convex regression in multidimensions: Suboptimality of least squares estimators
Gil Kur, Fuchang Gao, Adityanand Guntuboyina, and Bodhisattva Sen. Convex regression in multidimensions: Suboptimality of least squares estimators . The Annals of Statistics, 52 0 (6): 0 2791--2815, 2024
2024
-
[32]
On convergence rates of convex regression in multiple dimensions
Eunji Lim. On convergence rates of convex regression in multiple dimensions. INFORMS Journal on Computing, 26 0 (3): 0 616--628, 2014
2014
-
[33]
Convex regression with a penalty
Eunji Lim. Convex regression with a penalty. arXiv preprint arXiv:2509.19788v1, 2025
2025
-
[34]
E. J. McShane . Extension of range of functions. Bulletin of the American Mathematical Society, 40 0 (12): 0 837--842, 1934
1934
-
[35]
E. A. Nadaraya. On estimating regression. Theory of Probability & Its Applications, 9 0 (1): 0 141--142, 1964
1964
-
[36]
Nesterov
Yu. Nesterov. Smooth minimization of non-smooth functions. Mathematical Programming, 103 0 (1): 0 127--152, 2005
2005
-
[37]
Interior-Point Polynomial Algorithms in Convex Programming
Yurii Nesterov and Arkadii Nemirovskii. Interior-Point Polynomial Algorithms in Convex Programming. Society for Industrial and Applied Mathematics, 1994
1994
-
[38]
Jorge Nocedal and Stephen J. Wright. Numerical Optimization. Springer, 2006
2006
-
[39]
Max-min representation of piecewise linear functions
Sergei Ovchinnikov. Max-min representation of piecewise linear functions. Contributions to Algebra and Geometry, 43 0 (1): 0 297--302, 2002
2002
-
[40]
Scikit-learn: Machine learning in python
Fabian Pedregosa, Ga \"e l Varoquaux, Alexandre Gramfort, Vincent Michel, Bertrand Thirion, Olivier Grisel, Mathieu Blondel, Peter Prettenhofer, Ron Weiss, Vincent Dubourg, Jake Vanderplas, Alexandre Passos, David Cournapeau, Matthieu Brucher, Matthieu Perrot, and \'E douard D...
2011
-
[41]
Torrecilla, Miguel Carbajo Berrocal, Pablo Marcos Manchón, and Alberto Suárez
Carlos Ramos-Carreño, José L. Torrecilla, Miguel Carbajo Berrocal, Pablo Marcos Manchón, and Alberto Suárez. scikit-fda: A P ython P ackage for F unctional D ata A nalysis . Journal of Statistical Software, 109 0 (2): 0 1--37, 2024
2024
-
[42]
Piecewise linear regression via a difference of convex functions
Ali Siahkamari, Aditya Gangrade, Brian Kulis, and Venkatesh Saligrama. Piecewise linear regression via a difference of convex functions. In Hal Daumé III and Aarti Singh, editors, Proceedings of the 37th International Conference on Machine Learning, volume 119 of Proceedings o...
2020
-
[43]
Charles J. Stone. Optimal global rates of convergence for nonparametric regression. The Annals of Statistics, 10 0 (4): 0 1040--1053, 1982
1982
-
[44]
Least squares estimation of weakly convex functions
Sun Sun and Yaoliang Yu. Least squares estimation of weakly convex functions. In Kamalika Chaudhuri and Masashi Sugiyama, editors, Proceedings of the Twenty-Second International Conference on Artificial Intelligence and Statistics, volume 89 of Proceedings of Machine Learning ...
2019
-
[45]
Fitting piecewise linear continuous functions
Alejandro Toriello and Juan Pablo Vielma. Fitting piecewise linear continuous functions. European Journal of Operational Research, 219 0 (1): 0 86--95, 2012
2012
-
[46]
Empirical Processes in M-Estimation
Sara van de Geer. Empirical Processes in M-Estimation. Cambridge University Press, 2000
2000
-
[47]
Strong and weak convexity of sets and functions
Jean-Philippe Vial. Strong and weak convexity of sets and functions. Mathematics of Operations Research, 8 0 (2): 0 231--259, 1983
1983
-
[48]
Wainwright
Martin J. Wainwright. High-Dimensional Statistics: A Non-Asymptotic Viewpoint. Cambridge University Press, 2019
2019
-
[49]
Geoffrey S. Watson. Smooth regression analysis. Sankhyā: The Indian Journal of Statistics, Series A (1961-2002), 26 0 (4): 0 359--372, 1964
1961
-
[50]
C. T. Zahn. Black box maximization of circular coverage. Journal of Research of the National Bureau of Standards -- B. Mathematics and Mathematical Physics, 66B 0 (4): 0 181--216, 1962
1962
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.