Related papers: Index of a matrix, complex logarithms, and multidi…

The computation of eigenvalues of large-scale matrices arising from finite element discretizations has gained significant interest in the last decade. Here we present a new algorithm based on slicing the spectrum that takes advantage of the…

Let M be an irreducible Riemannian symmetric space. The index i(M) of M is the minimal codimension of a (non-trivial) totally geodesic submanifold of M. The purpose of this note is to determine the index i(M) for all irreducible Riemannian…

In this paper, the integral $\pmatrix{\lambda_1 &\lambda_2 &\lambda_3\cr 0 &0 &0\cr}\, \int_0^\infty \, r^{\lambda_3+2}\, \exp{(-\alpha r^2)}\, j_{\lambda_1}(k_1r) \,j_{\lambda_2}(k_2r) \,dr$, where $k_1$, $k_2$ and $\alpha$ are positive,…

We propose a multi-index algorithm for the Monte Carlo (MC) discretization of a linear, elliptic PDE with affine-parametric input. We prove an error vs. work analysis which allows a multi-level finite-element approximation in the physical…

Consider an $n\times n$ matrix polynomial $P(\lambda)$ and a set $\Sigma$ consisting of $k \le n$ distinct complex numbers. In this paper, a (weighted) spectral norm distance from $P(\lambda)$ to the matrix polynomials whose spectra include…

Suppose that we are given a string $s$ of length $n$ over an alphabet $\{0,1,\ldots,n^{O(1)}\}$ and $\delta$ is the string complexity of $s$, a known compression measure. We describe an index on $s$ with $O(\delta\log\frac{n}{\delta})$…

The spectral properties of two special classes of Jacobi operators are studied. For the first class represented by the $2M$-dimensional real Jacobi matrices whose entries are symmetric with respect to the secondary diagonal, a new…

We extend our recent result [Cipolloni, Erd\H{o}s, Schr\"oder 2019] on the central limit theorem for the linear eigenvalue statistics of non-Hermitian matrices $X$ with independent, identically distributed complex entries to the real…

We investigate the problem of pointwise convergence of the family of non-linear integral operators: \begin{equation} L_\lambda(f,x) = \int_a^b \sum_{m=1}^N f^m(t) K_{\lambda ,m}(x,t) dt, \end{equation} where $\lambda $ is a real parameters,…

Let $\boldsymbol{\Sigma}_N$ be a $M \times N$ random matrix defined by $\boldsymbol{\Sigma}_N = \mathbf{B}_N + \sigma \mathbf{W}_N$ where $\mathbf{B}_N$ is a uniformly bounded deterministic matrix and where $\mathbf{W}_N$ is an independent…

We present a real symmetric tri-diagonal matrix of order $n$ whose eigenvalues are $\{2k \}_{k=0}^{n-1}$ which also satisfies the additional condition that its leading principle submatrix has a uniformly interlaced spectrum, $\{2l + 1…

In this paper we study higher-order difference equations which can be written as follows: $$ \mathbf{J} (y_0,y_1,...)^T = \lambda^N (y_0,y_1,...)^T, $$ where $\mathbf{J}$ is a $(2N+1)$-diagonal bounded banded matrix…

Matrix sensing has many real-world applications in science and engineering, such as system control, distance embedding, and computer vision. The goal of matrix sensing is to recover a matrix $A_\star \in \mathbb{R}^{n \times n}$, based on a…

We develop an accelerated gradient descent algorithm on the Grassmann manifold to compute the subspace spanned by a number of leading eigenvectors of a symmetric positive semi-definite matrix. This has a constant cost per iteration and a…

We give a class of examples of $A$-hypergeometric systems that display integrality of mirror maps. Specifically, these systems have solutions $F(\lambda_1,\dots,\lambda_N) = 1$ and $\log\lambda^l + G(\lambda_1,\dots,\lambda_N)$ (for certain…

We develop a potential-theoretic and functional framework for the fractional--logarithmic Laplacian $(-\Delta)^{s+\ln}$ and its inhomogeneous counterpart $(\lambda I-\Delta)^{s+\ln}$ with $\lambda>1$. Their inverses yield logarithmic…

In this work we revisit the problem of solving multi-matrix systems through numerical large $N$ methods. The framework is a collective, loop space representation which provides a constrained optimization problem, addressed through…

The use of realistic input models has gained popularity in the theory community. Assuming a realistic input model often precludes complicated hypothetical inputs, and the analysis yields bounds that better reflect the behaviour of…

Let $\a$ be a complex random variable with mean zero and bounded variance $\sigma^{2}$. Let $N_{n}$ be a random matrix of order $n$ with entries being i.i.d. copies of $\a$. Let $\lambda_{1}, ..., \lambda_{n}$ be the eigenvalues of…

We derive a numerical algorithm for evaluating the Riemannian logarithm on the Stiefel manifold with respect to the canonical metric. In contrast to the existing optimization-based approach, we work from a purely matrix-algebraic…