Related papers: A tridiagonal matrix-valued process with stochasti…
The problem of computing recurrence coefficients of sequences of rational functions orthogonal with respect to a discrete inner product is formulated as an inverse eigenvalue problem for a pencil of Hessenberg matrices. Two procedures are…
We consider nonparametric invariant density and drift estimation for a class of multidimensional degenerate resp. hypoelliptic diffusion processes, so-called stochastic damping Hamiltonian systems or kinetic diffusions, under anisotropic…
For any $\beta>0$, we provide a tridiagonal matrix model and compute the joint eigenvalue density of a random rank one non-Hermitian perturbation of Gaussian and Laguerre $\beta$-ensembles of random matrices.
We develop an iterative refinement method that improves the accuracy of a user-chosen subset of $k$ eigenvectors ($k\ll n$) of an $n\times n$ real symmetric matrix. Using an orthogonal matrix represented in compact WY form, the method…
We study random matrices acting on tensor product spaces which have been transformed by a linear block operation. Using operator-valued free probability theory, under some mild assumptions on the linear map acting on the blocks, we compute…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
We present efficient numerical techniques for calculation of eigenvalue distributions of random matrices in the beta-ensembles. We compute histograms using direct simulations on very large matrices, by using tridiagonal matrices with…
In this work we introduce a theory of stochastic integration for operator-valued integrands with respect to some classes of cylindrical martingale-valued measures in Hilbert spaces. The integral is constructed via the radonification of…
We study the expectation value of the logarithm of the partition function of large binary-to-binary lattice-gas Restricted Boltzmann Machines (RBMs) within a replica-symmetric ansatz, averaging over the disorder represented by the…
We develop an improved version of the stochastic semigroup approach to study the edge of $\beta$-ensembles pioneered by Gorin and Shkolnikov, and later extended to rank-one additive perturbations by the author and Shkolnikov. Our method is…
We study the off-diagonal matrix elements of observables that break the translational symmetry of a spin-chain Hamiltonian, and as such connect energy eigenstates from different total quasimomentum sectors. We consider quantum-chaotic and…
A randomized algorithm for computing a data sparse representation of a given rank structured matrix $A$ (a.k.a. an $H$-matrix) is presented. The algorithm draws on the randomized singular value decomposition (RSVD), and operates under the…
We obtain the asymptotic distribution of eigenvalues of real symmetric tridiagonal matrices as their dimension increases to infinity and whose diagonal and off-diagonal elements asymptotically change with the index n as J_{nt+i nt+i}\sim…
We study triangulation schemes for the joint kernel of a diffusion process with uniformly continuous coefficients and an adapted, non-resonant Abelian process. The prototypical example of Abelian process to which our methods apply is given…
The real symmetric tridiagonal eigenproblem is of outstanding importance in numerical computations; it arises frequently as part of eigensolvers for standard and generalized dense Hermitian eigenproblems that are based on a reduction to…
In this paper we prove that there exists an asymptotical diagonalization algorithm for a class of sparse Hermitian (or real symmetric) matrices if and only if the matrices become Hessenberg matrices after some permutation of rows and…
We propose the supervised hierarchical Dirichlet process (sHDP), a nonparametric generative model for the joint distribution of a group of observations and a response variable directly associated with that whole group. We compare the sHDP…
We present an extension of the summation-by-parts (SBP) framework to tensor-product spectral-element operators in collapsed coordinates. The proposed approach enables the construction of provably stable discretizations of arbitrary order…
We develop an efficient algorithm for sampling the eigenvalues of random matrices distributed according to the Haar measure over the orthogonal or unitary group. Our technique samples directly a factorization of the Hessenberg form of such…
In this paper, we answer a question posed by Kurt Johansson, to find a PDE for the joint distribution of the Airy Process. The latter is a continuous stationary process, describing the motion of the outermost particle of the Dyson Brownian…