Structure Preserving Algorithm for the Logarithm of Symplectic Matrices

Authors : Bahar Arslan

Pages : 133-141

Doi:10.36753/mathenot.868902

View : 18 | Download : 19

Publication Date : 2021-09-30

Article Type : Research

Abstract :The current algorithms use either the full form or the Schur decomposition of the matrix in the inverse scaling and squaring method to compute the matrix logarithm. The inverse scaling and squaring method consists of two main calculations: taking a square root and evaluating the Padé approximants. In this work, we suggest using the structure preserving iteration as an alternative to Denman-Beavers iteration for taking a square root. Numerical experiments show that while using the structure preserving square root iteration in the inverse scaling and squaring method preserves the Hamiltonian structure of matrix logarithm, Denman-Beavers iteration and Schur decomposition cause a structure loss.
Keywords : Matrix functions, matrix logarithm, symplectic matrix, Hamiltonian matrix, inverse scaling and squaring method