Large dimensional Gram type matrices are common objects in high-dimensional statistics and machine learning. In this paper, we study the limiting distribution of the edge eigenvalues for a general class of high-dimensional Gram type random matrices, including separable sample covariance matrices, sparse sample covariance matrices, bipartite stochastic block model and random Gram matrices with general variance profiles. Specifically, we prove that under (almost) sharp moment conditions and certain tractable regularity assumptions, the edge eigenvalues, i.e., the largest few eigenvalues of non-spiked Gram type random matrices or the extremal bulk eigenvalues of spiked Gram type random matrices, satisfy the Tracy-Widom distribution asymptotically. Our results can be used to construct adaptive, accurate and powerful statistics for high-dimensional statistical inference. In particular, we propose data-dependent statistics to infer the number of signals under general noise structure, test the one-sided sphericity of separable matrix, and test the structure of bipartite stochastic block model. Numerical simulations show strong support of our proposed statistics. The core of our proof is to establish the edge universality and Tracy-Widom distribution for a rectangular Dyson Brownian motion with regular initial data. This is a general strategy to study the edge statistics for high-dimensional Gram type random matrices without exploring the specific independence structure of the target matrices. It has potential to be applied to more general random matrices that are beyond the ones considered in this paper.

In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method detects critical numerical fluxes which may lead to negative density and pressure, and then imposes a simple flux limiter combining the high-order numerical flux with the first-order Lax-Friedrichs flux to satisfy a sufficient condition for preserving positivity. Though an extra time-step size condition is required to maintain the formal order of accuracy, it is less restrictive than those in previous works. A number of numerical examples suggest that this method, when applied on an essentially non-oscillatory scheme, can be used to prevent positivity failure when the flow involves vacuum or near vacuum and very strong discontinuities.

No abstract uploaded!

We investigate the minimal and isoperimetric surface problems in a large class of sub-Riemannian manifolds, the so-called Vertically Rigid spaces. We construct an adapted connection for such spaces and, using the variational tools of Bryant, Griffiths and Grossman, derive succinct forms of the Euler-Lagrange equations for critical points for the associated variational problems. Using the Euler-Lagrange equations, we show that minimal and isoperimetric surfaces satisfy a constant horizontal mean curvature conditions away from characteristic points. Moreover, we use the formalism to construct a horizontal second fundamental form, II0, for vertically rigid spaces and, as a first application, use II0 to show that minimal surfaces cannot have points of horizontal positive curvature and that minimal surfaces in Carnot groups cannot be locally strictly horizontally geometrically convex. We note that the convexity condition is distinct from others currently in the literature.

No abstract uploaded!