Pub Date : 2024-06-22DOI: 10.1007/s00440-024-01288-y
Arnaud Debussche, Ruoyuan Liu, Nikolay Tzvetkov, Nicola Visciglia
We consider the nonlinear Schrödinger equation with multiplicative spatial white noise and an arbitrary polynomial nonlinearity on the two-dimensional full space domain. We prove global well-posedness by using a gauge-transform introduced by Hairer and Labbé (Electron Commun Probab 20(43):11, 2015) and constructing the solution as a limit of solutions to a family of approximating equations. This paper extends a previous result by Debussche and Martin (Nonlinearity 32(4):1147–1174, 2019) with a sub-quadratic nonlinearity.
{"title":"Global well-posedness of the 2D nonlinear Schrödinger equation with multiplicative spatial white noise on the full space","authors":"Arnaud Debussche, Ruoyuan Liu, Nikolay Tzvetkov, Nicola Visciglia","doi":"10.1007/s00440-024-01288-y","DOIUrl":"https://doi.org/10.1007/s00440-024-01288-y","url":null,"abstract":"<p>We consider the nonlinear Schrödinger equation with multiplicative spatial white noise and an arbitrary polynomial nonlinearity on the two-dimensional full space domain. We prove global well-posedness by using a gauge-transform introduced by Hairer and Labbé (Electron Commun Probab 20(43):11, 2015) and constructing the solution as a limit of solutions to a family of approximating equations. This paper extends a previous result by Debussche and Martin (Nonlinearity 32(4):1147–1174, 2019) with a sub-quadratic nonlinearity.\u0000</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"27 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141507312","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-15DOI: 10.1007/s00440-024-01285-1
Alexander Drewitz, Alexis Prévost, Pierre-François Rodriguez
For a large class of amenable transient weighted graphs G, we prove that the sign clusters of the Gaussian free field on G fall into a regime of strong supercriticality, in which two infinite sign clusters dominate (one for each sign), and finite sign clusters are necessarily tiny, with overwhelming probability. Examples of graphs belonging to this class include regular lattices such as ({mathbb {Z}}^d), for (dge 3), but also more intricate geometries, such as Cayley graphs of suitably growing (finitely generated) non-Abelian groups, and cases in which random walks exhibit anomalous diffusive behavior, for instance various fractal graphs. As a consequence, we also show that the vacant set of random interlacements on these objects, introduced by Sznitman (Ann Math 171(3):2039–2087, 2010), and which is intimately linked to the free field, contains an infinite connected component at small intensities. In particular, this result settles an open problem from Sznitman (Invent Math 187(3):645–706, 2012).
对于一大类可处理的瞬时加权图 G,我们证明了 G 上高斯自由场的符号簇属于强超临界状态,其中两个无限符号簇占主导地位(每个符号一个),而有限符号簇必然很小,具有压倒性概率。属于这一类的图的例子包括规则网格,如 ({mathbb {Z}}^d), for (dge 3), 但也包括更复杂的几何图形,如适当增长的(有限生成的)非阿贝尔群的卡莱图,以及随机漫步表现出异常扩散行为的情况,如各种分形图。因此,我们还证明了由 Sznitman(Ann Math 171(3):2039-2087,2010 年)引入的、与自由场密切相关的这些对象上的随机置换空集,在小强度下包含一个无限连通分量。特别是,这一结果解决了 Sznitman 提出的一个未决问题(Invent Math 187(3):645-706, 2012)。
{"title":"Geometry of Gaussian free field sign clusters and random interlacements","authors":"Alexander Drewitz, Alexis Prévost, Pierre-François Rodriguez","doi":"10.1007/s00440-024-01285-1","DOIUrl":"https://doi.org/10.1007/s00440-024-01285-1","url":null,"abstract":"<p>For a large class of amenable transient weighted graphs <i>G</i>, we prove that the sign clusters of the Gaussian free field on <i>G</i> fall into a regime of <i>strong supercriticality</i>, in which two infinite sign clusters dominate (one for each sign), and finite sign clusters are necessarily tiny, with overwhelming probability. Examples of graphs belonging to this class include regular lattices such as <span>({mathbb {Z}}^d)</span>, for <span>(dge 3)</span>, but also more intricate geometries, such as Cayley graphs of suitably growing (finitely generated) non-Abelian groups, and cases in which random walks exhibit anomalous diffusive behavior, for instance various fractal graphs. As a consequence, we also show that the vacant set of random interlacements on these objects, introduced by Sznitman (Ann Math 171(3):2039–2087, 2010), and which is intimately linked to the free field, contains an infinite connected component at small intensities. In particular, this result settles an open problem from Sznitman (Invent Math 187(3):645–706, 2012).</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"54 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-06-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141507313","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-07DOI: 10.1007/s00440-024-01292-2
Markus Tempelmayr
We give a novel characterization of the centered model in regularity structures which persists for rough drivers even as a mollification fades away. We present our result for a class of quasilinear equations driven by noise, however we believe that the method is robust and applies to a much broader class of subcritical equations. Furthermore, we prove that a convergent sequence of noise ensembles, satisfying uniformly a spectral gap assumption, implies the corresponding convergence of the associated models. Combined with the characterization, this establishes a universality-type result.
{"title":"Characterizing models in regularity structures: a quasilinear case","authors":"Markus Tempelmayr","doi":"10.1007/s00440-024-01292-2","DOIUrl":"https://doi.org/10.1007/s00440-024-01292-2","url":null,"abstract":"<p>We give a novel characterization of the centered model in regularity structures which persists for rough drivers even as a mollification fades away. We present our result for a class of quasilinear equations driven by noise, however we believe that the method is robust and applies to a much broader class of subcritical equations. Furthermore, we prove that a convergent sequence of noise ensembles, satisfying uniformly a spectral gap assumption, implies the corresponding convergence of the associated models. Combined with the characterization, this establishes a universality-type result.</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"64 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-06-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141550772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-06DOI: 10.1007/s00440-024-01278-0
Julien Chhor, Suzanne Sigalla, Alexandre B. Tsybakov
We study benign overfitting in the setting of nonparametric regression under mean squared risk, and on the scale of Hölder classes. We construct a local polynomial estimator of the regression function that is minimax optimal on a Hölder class with any given smoothness, and that is a continuous function interpolating the set of observations with high probability. The key element of the construction is the use of singular kernels. Moreover, we prove that adaptation to unknown smoothness is compatible with benign overfitting. Namely, we construct a continuous and interpolating local polynomial estimator attaining the minimax optimal rate in (L_2) adaptively to the unknown Hölder smoothness. Our results highlight the fact that interpolation can be fundamentally decoupled from bias-variance tradeoff in the problem of nonparametric regression.
{"title":"Benign overfitting and adaptive nonparametric regression","authors":"Julien Chhor, Suzanne Sigalla, Alexandre B. Tsybakov","doi":"10.1007/s00440-024-01278-0","DOIUrl":"https://doi.org/10.1007/s00440-024-01278-0","url":null,"abstract":"<p>We study benign overfitting in the setting of nonparametric regression under mean squared risk, and on the scale of Hölder classes. We construct a local polynomial estimator of the regression function that is minimax optimal on a Hölder class with any given smoothness, and that is a continuous function interpolating the set of observations with high probability. The key element of the construction is the use of singular kernels. Moreover, we prove that adaptation to unknown smoothness is compatible with benign overfitting. Namely, we construct a continuous and interpolating local polynomial estimator attaining the minimax optimal rate in <span>(L_2)</span> adaptively to the unknown Hölder smoothness. Our results highlight the fact that interpolation can be fundamentally decoupled from bias-variance tradeoff in the problem of nonparametric regression.</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"1 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141550773","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-19DOI: 10.1007/s00440-024-01282-4
Jiaoyang Huang, Colin McSwiggen
This paper studies the asymptotic behavior of several central objects in Dunkl theory as the dimension of the underlying space grows large. Our starting point is the observation that a recent result from the random matrix theory literature implies a large deviations principle for the hydrodynamic limit of radial Dunkl processes. Using this fact, we prove a variational formula for the large-N asymptotics of generalized Bessel functions, as well as a large deviations principle for the more general family of radial Heckman–Opdam processes. As an application, we prove a theorem on the asymptotic behavior of weight multiplicities of irreducible representations of compact or complex simple Lie algebras in the limit of large rank. The theorems in this paper generalize several known results describing analogous asymptotics for Dyson Brownian motion, spherical matrix integrals, and Kostka numbers.
本文研究了当底层空间维度变大时,邓克尔理论中几个中心对象的渐近行为。我们的出发点是观察到随机矩阵理论文献中的一个最新结果隐含了径向邓克尔过程流体力学极限的大偏差原理。利用这一事实,我们证明了广义贝塞尔函数大 N 渐近线的变分公式,以及更一般的径向 Heckman-Opdam 过程族的大偏差原理。作为应用,我们证明了紧凑或复杂简单李代数不可还原表示的权乘在大秩极限的渐近行为定理。本文中的定理概括了描述戴森布朗运动、球形矩阵积分和科斯特卡数的类似渐近的几个已知结果。
{"title":"Asymptotics of generalized Bessel functions and weight multiplicities via large deviations of radial Dunkl processes","authors":"Jiaoyang Huang, Colin McSwiggen","doi":"10.1007/s00440-024-01282-4","DOIUrl":"https://doi.org/10.1007/s00440-024-01282-4","url":null,"abstract":"<p>This paper studies the asymptotic behavior of several central objects in Dunkl theory as the dimension of the underlying space grows large. Our starting point is the observation that a recent result from the random matrix theory literature implies a large deviations principle for the hydrodynamic limit of radial Dunkl processes. Using this fact, we prove a variational formula for the large-<i>N</i> asymptotics of generalized Bessel functions, as well as a large deviations principle for the more general family of radial Heckman–Opdam processes. As an application, we prove a theorem on the asymptotic behavior of weight multiplicities of irreducible representations of compact or complex simple Lie algebras in the limit of large rank. The theorems in this paper generalize several known results describing analogous asymptotics for Dyson Brownian motion, spherical matrix integrals, and Kostka numbers.\u0000</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"28 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141146390","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-16DOI: 10.1007/s00440-024-01286-0
Dan Mikulincer, Yair Shenfeld
Contraction properties of transport maps between probability measures play an important role in the theory of functional inequalities. The actual construction of such maps, however, is a non-trivial task and, so far, relies mostly on the theory of optimal transport. In this work, we take advantage of the infinite-dimensional nature of the Gaussian measure and construct a new transport map, based on the Föllmer process, which pushes forward the Wiener measure onto probability measures on Euclidean spaces. Utilizing the tools of the Malliavin and stochastic calculus in Wiener space, we show that this Brownian transport map is a contraction in various settings where the analogous questions for optimal transport maps are open. The contraction properties of the Brownian transport map enable us to prove functional inequalities in Euclidean spaces, which are either completely new or improve on current results. Further and related applications of our contraction results are the existence of Stein kernels with desirable properties (which lead to new central limit theorems), as well as new insights into the Kannan–Lovász–Simonovits conjecture. We go beyond the Euclidean setting and address the problem of contractions on the Wiener space itself. We show that optimal transport maps and causal optimal transport maps (which are related to Brownian transport maps) between the Wiener measure and other target measures on Wiener space exhibit very different behaviors.
{"title":"The Brownian transport map","authors":"Dan Mikulincer, Yair Shenfeld","doi":"10.1007/s00440-024-01286-0","DOIUrl":"https://doi.org/10.1007/s00440-024-01286-0","url":null,"abstract":"<p>Contraction properties of transport maps between probability measures play an important role in the theory of functional inequalities. The actual construction of such maps, however, is a non-trivial task and, so far, relies mostly on the theory of optimal transport. In this work, we take advantage of the infinite-dimensional nature of the Gaussian measure and construct a new transport map, based on the Föllmer process, which pushes forward the Wiener measure onto probability measures on Euclidean spaces. Utilizing the tools of the Malliavin and stochastic calculus in Wiener space, we show that this Brownian transport map is a contraction in various settings where the analogous questions for optimal transport maps are open. The contraction properties of the Brownian transport map enable us to prove functional inequalities in Euclidean spaces, which are either completely new or improve on current results. Further and related applications of our contraction results are the existence of Stein kernels with desirable properties (which lead to new central limit theorems), as well as new insights into the Kannan–Lovász–Simonovits conjecture. We go beyond the Euclidean setting and address the problem of contractions on the Wiener space itself. We show that optimal transport maps and causal optimal transport maps (which are related to Brownian transport maps) between the Wiener measure and other target measures on Wiener space exhibit very different behaviors.\u0000</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"68 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-05-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141060038","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-15DOI: 10.1007/s00440-024-01281-5
Christian Mönch
We study a general class of percolation models in Euclidean space including long-range percolation, scale-free percolation, the weight-dependent random connection model and several other previously investigated models. Our focus is on the weak decay regime, in which inter-cluster long-range connection probabilities fall off polynomially with small exponent, and for which we establish several structural properties. Chief among them are the continuity of the bond percolation function and the transience of infinite clusters.
{"title":"Inhomogeneous long-range percolation in the weak decay regime","authors":"Christian Mönch","doi":"10.1007/s00440-024-01281-5","DOIUrl":"https://doi.org/10.1007/s00440-024-01281-5","url":null,"abstract":"<p>We study a general class of percolation models in Euclidean space including long-range percolation, scale-free percolation, the weight-dependent random connection model and several other previously investigated models. Our focus is on the <i>weak decay regime</i>, in which inter-cluster long-range connection probabilities fall off polynomially with small exponent, and for which we establish several structural properties. Chief among them are the continuity of the bond percolation function and the transience of infinite clusters.</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"131 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-05-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141060039","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-11DOI: 10.1007/s00440-024-01283-3
Sébastien Gouëzel, Jérôme Rousseau, Manuel Stadlbauer
We study the minimal distance between two orbit segments of length n, in a random dynamical system with sufficiently good mixing properties. This problem has already been solved in non-random dynamical system, and on average in random dynamical systems (the so-called annealed version of the problem): it is known that the asymptotic behavior for this question is given by a dimension-like quantity associated to the invariant measure, called correlation dimension (or Rényi entropy). We study the analogous quenched question, and show that the asymptotic behavior is more involved: two correlation dimensions show up, giving rise to a non-smooth behavior of the associated asymptotic exponent.
我们研究的是在具有足够好的混合特性的随机动力系统中,长度为 n 的两个轨道段之间的最小距离。这个问题已经在非随机动力学系统和随机动力学系统的平均值中得到了解决(即所谓的退火版问题):众所周知,这个问题的渐近行为是由一个与不变度量相关联的类似维度的量给出的,这个量被称为相关维度(或雷尼熵)。我们研究了类似的淬火问题,结果表明渐近行为更复杂:出现了两个相关维度,导致相关渐近指数的非平滑行为。
{"title":"Minimal distance between random orbits","authors":"Sébastien Gouëzel, Jérôme Rousseau, Manuel Stadlbauer","doi":"10.1007/s00440-024-01283-3","DOIUrl":"https://doi.org/10.1007/s00440-024-01283-3","url":null,"abstract":"<p>We study the minimal distance between two orbit segments of length <i>n</i>, in a random dynamical system with sufficiently good mixing properties. This problem has already been solved in non-random dynamical system, and on average in random dynamical systems (the so-called annealed version of the problem): it is known that the asymptotic behavior for this question is given by a dimension-like quantity associated to the invariant measure, called correlation dimension (or Rényi entropy). We study the analogous quenched question, and show that the asymptotic behavior is more involved: two correlation dimensions show up, giving rise to a non-smooth behavior of the associated asymptotic exponent.</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"8 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-05-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140929907","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-29DOI: 10.1007/s00440-024-01284-2
Orphée Collin, Giambattista Giacomin, Yueyun Hu
We consider the continuum version of the random field Ising model in one dimension: this model arises naturally as weak disorder scaling limit of the original Ising model. Like for the Ising model, a spin configuration is conveniently described as a sequence of spin domains with alternating signs (domain-wall structure). We show that for fixed centered external field and as spin-spin couplings become large, the domain-wall structure scales to a disorder dependent limit that coincides with the infinite disorder fixed point process introduced by D. S. Fisher in the context of zero temperature quantum Ising chains. In particular, our results establish a number of predictions that one can find in Fisher et al. (Phys Rev E 64:41, 2001). The infinite disorder fixed point process for centered external field is equivalently described in terms of the process of suitably selected extrema of a Brownian trajectory introduced and studied by Neveu and Pitman (in: Séminaire de probabilités XXIII. Lecture notes in mathematics, vol 1372, pp 239–247, 1989). This characterization of the infinite disorder fixed point is one of the important ingredients of our analysis.
我们考虑的是一维随机场伊辛模型的连续体版本:该模型是作为原始伊辛模型的弱无序缩放极限而自然产生的。与伊辛模型一样,自旋构型可以方便地描述为一系列符号交替的自旋畴(畴-壁结构)。我们的研究表明,对于固定中心的外场,当自旋-自旋耦合变得很大时,畴-壁结构会扩展到一个依赖于无序的极限,这与费雪(D. S. Fisher)在零温量子伊辛链中引入的无限无序定点过程相吻合。特别是,我们的结果证实了费舍尔等人(Phys Rev E 64:41, 2001)的一系列预测。以外部场为中心的无限无序定点过程等同于由 Neveu 和 Pitman 引入并研究的布朗轨迹的适当选择极值过程(in:Séminaire de probabilités XXIII.数学讲义,第 1372 卷,第 239-247 页,1989 年)。无限无序定点的这一特征是我们分析的重要内容之一。
{"title":"Infinite disorder renormalization fixed point for the continuum random field Ising chain","authors":"Orphée Collin, Giambattista Giacomin, Yueyun Hu","doi":"10.1007/s00440-024-01284-2","DOIUrl":"https://doi.org/10.1007/s00440-024-01284-2","url":null,"abstract":"<p>We consider the continuum version of the random field Ising model in one dimension: this model arises naturally as weak disorder scaling limit of the original Ising model. Like for the Ising model, a spin configuration is conveniently described as a sequence of spin domains with alternating signs (<i>domain-wall structure</i>). We show that for fixed centered external field and as spin-spin couplings become large, the domain-wall structure scales to a disorder dependent limit that coincides with the <i>infinite disorder fixed point</i> process introduced by D. S. Fisher in the context of zero temperature quantum Ising chains. In particular, our results establish a number of predictions that one can find in Fisher et al. (Phys Rev E 64:41, 2001). The infinite disorder fixed point process for centered external field is equivalently described in terms of the process of <i>suitably selected</i> extrema of a Brownian trajectory introduced and studied by Neveu and Pitman (in: Séminaire de probabilités XXIII. Lecture notes in mathematics, vol 1372, pp 239–247, 1989). This characterization of the infinite disorder fixed point is one of the important ingredients of our analysis.\u0000</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"44 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140838222","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-22DOI: 10.1007/s00440-024-01276-2
Han Huang, Konstantin Tikhomirov
The Gaussian elimination with partial pivoting (GEPP) is a classical algorithm for solving systems of linear equations. Although in specific cases the loss of precision in GEPP due to roundoff errors can be very significant, empirical evidence strongly suggests that for a typical square coefficient matrix, GEPP is numerically stable. We obtain a (partial) theoretical justification of this phenomenon by showing that, given the random (ntimes n) standard Gaussian coefficient matrix A, the growth factor of the Gaussian elimination with partial pivoting is at most polynomially large in n with probability close to one. This implies that with probability close to one the number of bits of precision sufficient to solve (Ax = b) to m bits of accuracy using GEPP is (m+O(log n)), which improves an earlier estimate (m+O(log ^2 n)) of Sankar, and which we conjecture to be optimal by the order of magnitude. We further provide tail estimates of the growth factor which can be used to support the empirical observation that GEPP is more stable than the Gaussian Elimination with no pivoting.
带部分支点的高斯消元法(GEPP)是求解线性方程组的经典算法。虽然在特定情况下,舍入误差会导致 GEPP 的精度损失非常大,但经验证据有力地表明,对于典型的平方系数矩阵,GEPP 在数值上是稳定的。我们通过证明给定随机(n/times n/)标准高斯系数矩阵 A,具有部分支点的高斯消元法的增长因子在 n 上最多为多项式大,概率接近于 1,从而为这一现象提供了(部分)理论依据。这意味着在概率接近于1的情况下,使用GEPP求解m比特精度的(Ax = b) 所需的精度比特数是(m+O(log n)),这改进了桑卡尔(Sankar)早先的估计值(m+O(log ^2 n)),我们猜想这个估计值在数量级上是最优的。我们进一步提供了增长因子的尾部估计值,可用于支持经验观察,即 GEPP 比无支点的高斯消除法更稳定。
{"title":"Average-case analysis of the Gaussian elimination with partial pivoting","authors":"Han Huang, Konstantin Tikhomirov","doi":"10.1007/s00440-024-01276-2","DOIUrl":"https://doi.org/10.1007/s00440-024-01276-2","url":null,"abstract":"<p>The Gaussian elimination with partial pivoting (GEPP) is a classical algorithm for solving systems of linear equations. Although in specific cases the loss of precision in GEPP due to roundoff errors can be very significant, empirical evidence strongly suggests that for a <i>typical</i> square coefficient matrix, GEPP is numerically stable. We obtain a (partial) theoretical justification of this phenomenon by showing that, given the random <span>(ntimes n)</span> standard Gaussian coefficient matrix <i>A</i>, the <i>growth factor</i> of the Gaussian elimination with partial pivoting is at most polynomially large in <i>n</i> with probability close to one. This implies that with probability close to one the number of bits of precision sufficient to solve <span>(Ax = b)</span> to <i>m</i> bits of accuracy using GEPP is <span>(m+O(log n))</span>, which improves an earlier estimate <span>(m+O(log ^2 n))</span> of Sankar, and which we conjecture to be optimal by the order of magnitude. We further provide tail estimates of the growth factor which can be used to support the empirical observation that GEPP is more stable than the Gaussian Elimination with no pivoting.\u0000</p>","PeriodicalId":20527,"journal":{"name":"Probability Theory and Related Fields","volume":"82 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140803415","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}