We consider the solution of large-scale Lyapunov and Stein equations. For Stein equations, the well-known Smith method will be adapted, with A_k = A^{2^k} not explicitly computed but in the recursive form A_k = A^{2^k}, and the fast growing but diminishing components in the approximate solutions truncated. Lyapunov equations will be first treated with the Cayley transform before the Smith method is applied. For algebraic equations with numerically low-ranked solutions of dimension <i>n</i>, the resulting algorithms are of an efficient <i>O</i>(<i>n</i>) computational complexity and memory requirement per iteration and converge essentially quadratically. An application in the estimation of a lower bound of the condition number for continuous-time algebraic Riccati equations is presented, as well as some numerical results.

We undertake a systematic study of the infinitesimal geometry of the Thurston metric, showing that the topology, convex geometry and metric geometry of the tangent and cotangent spheres based at any marked hyperbolic surface representing a point in Teichmüller space can recover the marking and geometry of this marked surface. We then translate the results concerning the infinitesimal structures to global geometric statements for the Thurston metric, most notably deriving rigidity statements for the Thurston metric analogous to the celebrated Royden theorem.

We develop a new method for deriving local laws for a large class of random matrices. It is applicable to many matrix models built from sums and products of deterministic or independent random matrices. In particular, it may be used to obtain local laws for matrix ensembles that are \emph{anisotropic} in the sense that their resolvents are well approximated by deterministic matrices that are not multiples of the identity. For definiteness, we present the method for sample covariance matrices of the form Q:=TXX∗T∗, where T is deterministic and X is random with independent entries. We prove that with high probability the resolvent of Q is close to a deterministic matrix, with an optimal error bound and down to optimal spectral scales. As an application, we prove the edge universality of Q by establishing the Tracy-Widom-Airy statistics of the eigenvalues of Q near the soft edges. This result applies in the single-cut and multi-cut cases. Further applications include the distribution of the eigenvectors and an analysis of the outliers and BBP-type phase transitions in finite-rank deformations; they will appear elsewhere. We also apply our method to Wigner matrices whose entries have arbitrary expectation, i.e. we consider W+A where W is a Wigner matrix and a Hermitian deterministic matrix. We prove the anisotropic local law for W+A and use it to establish edge universality.

No abstract uploaded!

No abstract uploaded!