Multiprecision algorithms for computing the matrix logarithm

Computation of the matrix logarithm in multiprecision.

You are now following this Submission

This code computes the logarithm and the Frechét derivative of a square matrix with non-negative eigenvalues in multiple precision. Both functions require the Advanpix Multiprecision Computing Toolbox for MATLAB (see www.advanpix.com) to work correctly.
The submission contains two functions:
* LOGM_MP computes the matrix logarithm and its Frechét derivative by using the inverse scaling and squaring algorithm and exploiting the complex Schur decomposition of the matrix.
* LOGM_MP_FULL computes the matrix logarithm without transforming the input matrix to upper triangular form.
A detailed description of the algorithms can be found in:
M. Fasi and N. J. Higham, Multiprecision Algorithms for Computing the Matrix Logarithm. Technical Report 2017.16, Manchester Institute for Mathematical Sciences, The University of Manchester, UK, May 2017.

Cite As

Massimiliano Fasi (2026). Multiprecision algorithms for computing the matrix logarithm (https://nl.mathworks.com/matlabcentral/fileexchange/63841-multiprecision-algorithms-for-computing-the-matrix-logarithm), MATLAB Central File Exchange. Retrieved .

Acknowledgements

Inspired by: The Matrix Function Toolbox

General Information

MATLAB Release Compatibility

  • Compatible with any release

Platform Compatibility

  • Windows
  • macOS
  • Linux
Version Published Release Notes Action
1.0.0.0

Add a picture.

Add diagonalization for normal matrices.
Fix performance bug in ALPHA(A,K,M).