Related papers: Fast Toeplitz eigenvalue computations, joining int…
We provide fast algorithms for overconstrained $\ell_p$ regression and related problems: for an $n\times d$ input matrix $A$ and vector $b\in\mathbb{R}^n$, in $O(nd\log n)$ time we reduce the problem $\min_{x\in\mathbb{R}^d} \|Ax-b\|_p$ to…
For two-loop two-point diagrams with arbitrary masses, an algorithm to derive the asymptotic expansion at large external momentum squared is constructed. By using a general theorem on asymptotic expansions of Feynman diagrams, the…
This paper studies two structured approximation problems: (1) Recovering a corrupted low-rank Toeplitz matrix and (2) recovering the range of a Fourier matrix from a single observation. Both problems are computationally challenging because…
Let $T_N$ denote an $N\times N$ Toeplitz matrix with finite, $N$ independent symbol ${\bf a}$. For $E_N$ a noise matrix satisfying mild assumptions (ensuring, in particular, that $N^{-1/2}\|E_N\|_{{\rm HS}}\to_{N\to\infty} 0$ at a…
One of the main applications of free probability is to show that for appropriately chosen independent copies of $d$ random matrix models, any noncommutative polynomial in these $d$ variables has a spectral distribution that converges…
A common claim of evolutionary computation methods is that they can achieve good results without the need for human intervention. However, one criticism of this is that there are still hyperparameters which must be tuned in order to achieve…
For compact self-adjoint operators in Hilbert spaces, two algorithms are proposed to provide fully computable a posteriori error estimate for eigenfunction approximation. Both algorithms apply well to the case of tight clusters and multiple…
An algorithm for computing an analytic function of a matrix $A$ is described. The algorithm is intended for the case where $A$ has some close eigenvalues, and clusters (subsets) of close eigenvalues are separated from each other. This…
We consider the asymptotic behavior of the eigenvalues of Toeplitz matrices with rational symbol as the size of the matrix goes to infinity. Our main result is that the weak limit of the normalized eigenvalue counting measure is a…
We study asymptotic properties of expectation propagation (EP) -- a method for approximate inference originally developed in the field of machine learning. Applied to generalized linear models, EP iteratively computes a multivariate…
This article focuses on the fluctuations of linear eigenvalue statistics of $T_{n\times p}T'_{n\times p}$, where $T_{n\times p}$ is an $n\times p$ Toeplitz matrix with real, complex or time-dependent entries. We show that as $n \rightarrow…
Asymptotically, we analytically derive the form of eigenvectors for two Fisher-Hartwig symbols besides those which were previously investigated in a $2016$ work due to Movassagh and Kadanoff, in which the authors characterized the…
We develop and analyze a single-loop algorithm for minimizing the sum of a Lipschitz differentiable function $f$, a prox-friendly proper closed function $g$ (with a closed domain on which $g$ is continuous) and the composition of another…
We present general prescriptions for the asymptotic expansion of massive multi-loop Feynman integrals near threshold. As in the case of previously known prescriptions for various limits of momenta and masses, the terms of the threshold…
Many models of interest in the natural and social sciences have no closed-form likelihood function, which means that they cannot be treated using the usual techniques of statistical inference. In the case where such models can be…
Let $\bigl\{X_k\bigr\}_{k \in \mathbb{Z}} \in \mathbb{L}^2(\mathcal{T})$ be a stationary process with associated lag operators ${\boldsymbol{\cal C}}_h$. Uniform asymptotic expansions of the corresponding empirical eigenvalues and…
We develop a theory of linear Mahler systems in several variables from the perspective of transcendence and algebraic independence, which also includes the possibility of dealing with several systems associated with sufficiently independent…
We introduce and analyze a continuous primal-dual dynamical system in the context of the minimization problem $f(x)+g(Ax)$, where $f$ and $g$ are convex functions and $A$ is a linear operator. In this setting, the trajectories of the…
We investigate real eigenvalues of real elliptic Ginibre matrices of size $n$, indexed by the parameter of asymmetry $\tau \in [0,1]$. In both the strongly and weakly non-Hermitian regimes, where $\tau \in [0,1)$ is fixed or…
We provide in this work an algorithm for approximating a very broad class of symmetric Toeplitz matrices to machine precision in $\mathcal{O}(n \log n)$ time with applications to fitting time series models. In particular, for a symmetric…