Related papers: Approximating the Perfect Sampling Grids for Compu…
An important problem in applications is the approximation of a function $f$ from a finite set of randomly scattered data $f(x_j)$. A common and powerful approach is to construct a trigonometric least squares approximation based on the set…
The classical orthogonal polynomials (Hermite, Laguerre and Jacobi) are involved in a vast number of applications in physics and engineering. When large degrees $n$ are needed, the use of recursion to compute the polynomials is not a good…
The computation of the matrix exponential is a ubiquitous operation in numerical mathematics, and for a general, unstructured $n\times n$ matrix it can be computed in $\mathcal{O}(n^3)$ operations. An interesting problem arises if the input…
We study the limiting spectral measure of large symmetric random matrices of linear algebraic structure. For Hankel and Toeplitz matrices generated by i.i.d. random variables $\{X_k\}$ of unit variance, and for symmetric Markov matrices…
We find two convergent series expansions for Legendre's first incomplete elliptic integral $F(\lambda,k)$ in terms of recursively computed elementary functions. Both expansions are valid at every point of the unit square $0<\lambda,k<1$.…
We present an approximation scheme for functions in three dimensions, that requires only their samples on the Cartesian grid, under the assumption that the functions are sufficiently concentrated in both space and frequency. The scheme is…
In this paper we develop algorithms for approximating matrix multiplication with respect to the spectral norm. Let A\in{\RR^{n\times m}} and B\in\RR^{n \times p} be two matrices and \eps>0. We approximate the product A^\top B using two…
This paper is concerned with the asymptotic behavior of the free energy for a class of Hermitean random matrix models, with odd degree polynomial potential, in the large N limit. It continues an investigation initiated and developed in a…
We formulate and partially prove a general conjecture providing necessary and sufficient conditions for the reality of the asymptotic spectrum of an arbitrary real banded block Toeplitz matrix. Additionally we present numerical experiments…
The current study investigates the asymptotic spectral properties of a finite difference approximation of nonlocal Helmholtz equations with a Caputo fractional Laplacian and a variable coefficient wave number $\mu$, as it occurs when…
We present an explicit solution of the eigen-spectrum Toeplitz matrix $C_{ij}= e^{- \kappa |i-j|}$ with $0\leq i,j \leq N$ and apply it to find analytically the plasma modes of a layered assembly of 2-dimensional electron gas. The solution…
A detailed analysis of the remainder obtained by truncating the Euler series up to the $n$th-order term is presented. In particular, by using an approach recently proposed by Weniger, asymptotic expansions of the remainder, both in inverse…
We present a sublinear time algorithm for computing a near optimal low-rank approximation to any positive semidefinite (PSD) Toeplitz matrix $T\in \mathbb{R}^{d\times d}$, given noisy access to its entries. In particular, given entrywise…
In this paper we derive uniform asymptotic expansions for the partial sums of the exponential series. We indicate how this information will be used in a later publication to obtain full and explicitly computable asymptotic expansions with…
Spectral properties of Toeplitz operators and their finite truncations have long been central in operator theory. In the finite dimensional, non-normal setting, the spectrum is notoriously unstable under perturbations. Random perturbations…
In this paper we consider symmetric, positive semidefinite (SPSD) matrix $A$ and present two algorithms for computing the $p$-Schatten norm $\|A\|_p$. The first algorithm works for any SPSD matrix $A$. The second algorithm works for…
The problem of recovering the configuration of points from their partial pairwise distances, referred to as the Euclidean Distance Matrix Completion (EDMC) problem, arises in a broad range of applications, including sensor network…
The paper pursues three objectives. Firstly, we provide an expanded version of spectral analysis of self-adjoint Toeplitz operators, initially built by M. Rosenblum in the 1960's. We offer some improvements to Rosenblum's approach: for…
Multivariate global polynomial approximations - such as polynomial chaos or stochastic collocation methods - are now in widespread use for sensitivity analysis and uncertainty quantification. The pseudospectral variety of these methods uses…
This paper studies Graphical SLOPE for precision matrix estimation, with emphasis on its ability to recover both sparsity and clusters of edges with equal or similar strength. In a fixed-dimensional regime, we establish that the root-$n$…