Abstract:In this paper we give the detailed error analysis of two algorithms $W_1$ and $W_2$ for computing the symplectic factorization of a symmetric positive definite and symplectic matrix $A \in \mathbb R^{2n \times 2n}$ in the form $A=LL^T$, where $L \in \mathbb R^{2n \times 2n}$ is a symplectic block lower triangular matrix. We prove that Algorithm $W_2$ is numerically stable for a broader class of symmetric positive definite matrices $A \in \mathbb R^{2n \times 2n}$. It means that Algorithm $W_2$ is producing the computed factors $\tilde L$ in floating-point arithmetic with machine precision $\mathcal{u}$ such that $||A-\tilde L {\tilde L}^T||_{2} = {\cal O}(\mathcal{u} ||{A}||_{2})$. On the other hand, Algorithm $W_1$ is unstable, in general, for symmetric positive definite and symplectic matrix $A$. In this paper we also give corresponding bounds for Algorithm $W_1$ that are weaker. We show that the factorization error depends on the condition number $\kappa_2(A_{11})$ of the principal submatrix $A_{11}$. Bounds for the loss of symplecticity of the lower block triangular matrices $L$ for both Algorithms $W_1$ and $W_2$ that hold in exact arithmetic for a broader class of symmetric positive definite matrices $A$ (but not necessarily symplectic) are also given. The tests performed in \textsl{MATLAB} illustrate that our error bounds for considered algorithms are reasonably sharp.

An effective criterion for a stable factorisation of strictly non-singular 2 × 2 matrix functions: use of the ExactMPF package

An effective criterion for a stable factorisation of strictly nonsingular 2x2 matrix functions. Utilisation of the ExactMPF package

Regular approximate factorization of a class of matrix-function with an unstable set of partial indices

Spectral Factorization of Rank-Deficient Polynomial Matrix-Functions

On Factor Left Prime Factorization Problems for Multivariate Polynomial Matrices

Avoiding breakdown in incomplete factorizations in low precision arithmetic

Refined multiplicative tensor product of matrix factorizations

Bauer's Spectral Factorization Method for Low Order Multiwavelet Filter Design

On incomplete Cholesky factorizations in half precision arithmetic

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

On Minor Left Prime Factorization Problem for Multivariate Polynomial Matrices

Combining Sparse Approximate Factorizations with Mixed-precision Iterative Refinement

Multivariate Polynomial Factorization by Interpolation Method

The numerical factorization of polynomials

Numerical stability of the symplectic $LL^T$ factorization

Exact Bivariate Polynomial Factorization Over A"E By Approximation Of Roots

Exact Bivariate Polynomial Factorization over ℚ by Approximation of Roots

Factorization of positive-semidefinite operators with absolutely summable entries

Towards Optimal Bounds for Implicit Factorization Problem

On multivariable matrix spectral factorization method