Testing matrix function algorithms using identities

Deadman, Edvin and Higham, Nicholas J. (2014) Testing matrix function algorithms using identities. [MIMS Preprint]

This is the latest version of this item.

[thumbnail of DeadmanHigham.pdf] PDF
DeadmanHigham.pdf
Download (274kB)

Abstract

Algorithms for computing matrix functions are typically tested by comparing the forward error with the product of the condition number and the unit roundoff. The forward error is computed with the aid of a reference solution, typically computed at high precision. An alternative approach is to use functional identities such as the ``round trip tests'' $e^{\log A} = A$ and $(A^{1/p})^p = A$, as are currently employed in a SciPy test module. We show how a linearized perturbation analysis for a functional identity allows the determination of a maximum residual consistent with backward stability of the constituent matrix function evaluations. Comparison of this maximum residual with a computed residual provides a necessary test for backward stability. We also show how the actual linearized backward error for these relations can be computed. Our approach makes use of Fr\'echet derivatives and estimates of their norms. Numerical experiments show that the proposed approaches are able both to detect instability and to confirm stability.

Item Type: MIMS Preprint
Additional Information: To appear in ACM Transactions on Mathematical Software
Uncontrolled Keywords: matrix function, normwise relative error, functional identity, testing, forward error, backward error, SciPy, Python, MATLAB, NAG Library
Subjects: MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User: Dr Edvin Deadman
Date Deposited: 16 Jan 2015
Last Modified: 08 Nov 2017 18:18
URI: https://eprints.maths.manchester.ac.uk/id/eprint/2242

Available Versions of this Item

Actions (login required)

View Item View Item