Related papers: A Jacobi-Type Eigensolver for Diagonally Dominant …
The hierarchical interpolative factorization for elliptic partial differential equations is a fast algorithm for approximate sparse matrix inversion in linear or quasilinear time. Its accuracy can degrade, however, when applied to strongly…
We consider multi-agent, convex optimization programs subject to separable constraints, where the constraint function of each agent involves only its local decision vector, while the decision vectors of all agents are coupled via a common…
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…
In this paper a two-sided, parallel Kogbetliantz-type algorithm for the hyperbolic singular value decomposition (HSVD) of real and complex square matrices is developed, with a single assumption that the input matrix, of order $n$, admits…
The aim of this paper is to study differential properties of orthogonal polynomials with respect to a discrete Jacobi-Sobolev bilinear form with mass point at $-1$ and/or $+1$. In particular, we construct the orthogonal polynomials using…
When the classical Hamburger moment problem has solutions, it has either exactly one solution or infinitely many solutions. Correspondingly, the moment problem is said to be either determinate or indeterminate. In terms of Jacobi operators,…
We prove that on the spectrum the integrated density of states (IDS for short) of periodic Jacobi matrices is related to the discriminant. The method is to count the number of generalized zeros of Bloch wave solutions.
We consider the dynamic problems for the discrete systems with discrete time associated with finite and semi-infinite Jacobi matrices. The result of the paper is a procedure of association of special Hilbert spaces of functions, namely de…
We present a new matrix-valued isospectral ordinary differential equation that asymptotically block-diagonalizes $n\times n$ zero-diagonal Jacobi matrices employed as its initial condition. This o.d.e.\ features a right-hand side with a…
In this paper, we compute the expectation of traces of powers of the hermitian matrix Jacobi process for a large enough but fixed size. To proceed, we first derive the semi-group density of its eigenvalues process as a bilinear series of…
We describe the relation between block Jacobi matrices and minimization problems for discrete time optimal control problems. Using techniques developed for the continuous case, we provide new algorithms to compute spectral invariants of…
Iterative methods with certified convergence for the computation of Gauss--Jacobi quadratures are described. The methods do not require a priori estimations of the nodes to guarantee its fourth-order convergence. They are shown to be…
In this paper, we propose a fast distributed solver for linear equations given by symmetric diagonally dominant M-Matrices. Our approach is based on a distributed implementation of the parallel solver of Spielman and Peng by considering a…
In this paper, we derive a handable expression for the Jacobi process semi group which is given by a bilinear series involving Jacobi polynomials. Our attempt uses a subordination of the considered process by means of a suitable random…
$k$-diagonal circulant matrices and cyclic banded matrices are widely used in numerical simulations and signal processing of circular linear systems. Algorithms that directly involve or specify linear or quadratic complexity for the…
In this paper we propose an efficiently preconditioned Newton method for the computation of the leftmost eigenpairs of large and sparse symmetric positive definite matrices. A sequence of preconditioners based on the BFGS update formula is…
For an arbitrary Hermitian period-$T$ Jacobi operator, we assume a perturbation by a Wigner-von Neumann type potential to devise subordinate solutions to the formal spectral equation for a (possibly infinite) real set, $S$, of the spectral…
Our main result asserts that a certain natural non-linear operator on Jacobi matrices built by a hyperbolic polynomial with real Julia set is a contraction in operator norm if the polynomial is sufficiently hyperbolic. This allows us to get…
We study the bit complexity of two methods, related to the Euclidean algorithm, for computing cubic and quartic analogs of the Jacobi symbol. The main bottleneck in such procedures is computation of a quotient for long division. We give…
The paper presents a comparative analysis of iterative numerical methods of Jacobi and Gauss-Seidel for solving systems of linear algebraic equations (SLAEs) with complex and real matrices. The ranges of convergence for both methods for…