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.

Optimal unit triangular factorization of symplectic matrices

Factorization and Symplectic Uniton Numbers for Harmonic Maps into Symplectic Groups

Constructive Factorizations of Imitons Via Singular Darboux Transformations

Interpolative Butterfly Factorization

Numerical stability of the symplectic $LL^T$ factorization

Matrix factorizations for domestic triangle singularities

Zero-sum Triangles for Involutory, Idempotent, Nilpotent and Unipotent Matrices

A Lower-Upper-Lower Block Triangular Decomposition with Minimal Off-Diagonal Ranks

Fifty Three Matrix Factorizations: A systematic approach

Butterfly factorization with error guarantees

Joint Tensor Factorization and Outlying Slab Suppression With Applications

Complexity estimates for triangular hierarchical matrix algorithms

Matrix factorizations and pentagon maps

Symmetrically Factorizable Groups and Set-Theoretical Solutions of the Pentagon Equation

On Some New Approximate Factorization Methods for Block Tridiagonal Matrices Suitable for Vector and Parallel Processors

Stability Analysis of Block LU Factorization for Complex Symmetric Block Tridiagonal Matrices

An efficient numerical algorithm for solving linear systems with cyclic tridiagonal coefficient matrices

Parallelization of incomplete factorization preconditioning of block tridiagonal linear systems with 1-D domain decomposition

Butterfly factorization via randomized matrix-vector multiplications

Sparse Factorization of Large Square Matrices

Communication-optimal parallel and sequential QR and LU factorizations: theory and practice