Related papers: Reconstruction of tridiagonal matrices from spectr…
In this article, we study several reconstruction methods for the inverse source problem of photoacoustic tomography (PAT) with spatially variable sound speed and damping. The backbone of these methods is the adjoint operators, which we…
The paper studies the global convergence of the block Jacobi me\-thod for symmetric matrices. Given a symmetric matrix $A$ of order $n$, the method generates a sequence of matrices by the rule $A^{(k+1)}=U_k^TA^{(k)}U_k$, $k\geq0$, where…
Transforming a matrix over a field to echelon form, or decomposing the matrix as a product of structured matrices that reveal the rank profile, is a fundamental building block of computational exact linear algebra. This paper surveys the…
We study stability of spectral types for semi-infinite self-adjoint tridiagonal matrices under random decaying perturbations. We show that absolutely continuous spectrum associated with bounded eigenfunctions is stable under Hilbert-Schmidt…
We introduce and analyze parallelizable algorithms to compress and accurately reconstruct finite simplicial complexes that have non-trivial automorphisms. The compressed data -- called a complex of groups -- amounts to a functor from (the…
We present theory and practice for robust implementations of bivariate Jacobi set and Reeb space algorithms. Robustness is a fundamental topic in computational geometry that deals with the issues of numerical errors and degenerate cases in…
The problem to be studied in this work is within the context of coefficient identification problems for the wave equation. More precisely, we consider the problem of reconstruction of the refractive index (or equivalently, the dielectric…
In this paper, we study complex Jacobi matrices obtained by the Christoffel and Geronimus transformations at a nonreal complex number, including the properties of the corresponding sequences of orthogonal polynomials. We also present some…
We describe a recursive algorithm that decomposes an algebraic set into locally closed equidimensional sets, i.e. sets which each have irreducible components of the same dimension. At the core of this algorithm, we combine ideas from the…
We study structural properties of the Lyapunov exponent $\gamma$ and the density of states $k$ for ergodic (or just invariant) Jacobi matrices in a general framework. In this analysis, a central role is played by the function…
In a recent paper, a new method was proposed to find the common invariant subspaces of a set of matrices. This paper invstigates the more general problem of putting a set of matrices into block triangular or block-diagonal form…
For tridiagonal block Jacobi operators, the standard transfer operator techniques only work if the off-diagonal entries are invertible. Under suitable assumptions on the range and kernel of these off-diagonal operators which assure a…
We consider the inverse spectral problem for periodic Jacobi matrices in terms of the vertical slits on the quasi-momentum domain plus the Dirichlet eigenvalues, i.e., the Marchenko-Ostrovsky mapping. Moreover, we show that the gradients of…
We show that on a two-dimensional compact nontrapping Riemannian manifold with strictly convex boundary, a piecewise constant function can be recovered from its integrals over geodesics. We adapt the injectivity proof which uses variations…
In this work, a new algorithm for solving symmetric indefinite systems of linear equations is presented. It factorizes the matrix into the form LDLt using Jacobi rotations in order to increase the pivot's absolute value. Furthermore, Rook's…
We consider the Hamburger, Stieltjes and Hausdorff moment problems, that are problems of the construction of a Borel measure supported on a real line, on a half-line or on an interval $(0,1)$, from a prescribed set of moments. We propose a…
In this paper a recursive algorithm is presented for evaluating multivariate Pad\'e approximants (of the rectangular type described in the work of Lutterodt) which is analogous to the Jacobi formula for univariate Pad\'e approximants. This…
We present here necessary and sufficient conditions for the invertibility of circulant and symmetric matrices that depend on three parameters and moreover, we explicitly compute the inverse. The techniques we use are related with the…
A real persymmetric Jacobi matrix of order $n$ whose eigenvalues are $2k^2$ $(k=0, ..., n-1)$ is presented, with entries given as explicit functions of $n$. Besides the possible use for testing forward and inverse numerical algorithms, such…
This paper proposes a harmonic Lanczos bidiagonalization method for computing some interior singular triplets of large matrices. It is shown that the approximate singular triplets are convergent if a certain Rayleigh quotient matrix is…