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.

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

Stability of Block LU Factorization for Block Tridiagonal Block H-matrices.

Necessary and sufficient condition for quadratic stability of interval matrices

Numerical stability of the symplectic $LL^T$ factorization

On Parametric Factorization of Bi-Orthogonal Laurent Polynomial Wavelet Filters

Block incomplete LU factorization for block-tridiagonal M-matrices

Verified Error Bounds for Matrix Decompositions

Stability and Fast Algorithms of Incomplete LU Factorization with Zero-Fill for Nine-Diagonal Matrices

D-stability Analysis of Polynomial Matrices and Polynomial Matrix Polytopes by Using LMI Approach

Necessary And Sufficient Conditions For Existence of the LU Factorization of an Arbitrary Matrix

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

Quadratic transformation and matrix biorthogonal polynomials: an $\mathcal{LU}$ factorization approach

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

Single-step triangular splitting iteration method for a class of complex symmetric linear system

High Performance Block Incomplete LU Factorization

Butterfly factorization by algorithmic identification of rank-one blocks

PERTURBATION THEORY FOR THE LU AND QR FACTORIZATIONS

Algebraic method for LU decomposition of dual quaternion matrix and its corresponding structure-preserving algorithm

The stability region for Schur stable trinomials with general complex coefficients

Block Maps and Fourier Analysis

New Block Triangular Preconditioners for Saddle Point Linear Systems with Highly Singular (1,1) Blocks