## 2006.155: Perturbation Theory and Backward Error for $AX-XB=C$

2006.155:
Nicholas J. Higham
(1993)
*Perturbation Theory and Backward Error for $AX-XB=C$.*
BIT Numerical Mathematics, 33.
pp. 124-136.
ISSN 1572-9125

Full text available as:

PDF - Archive staff only - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 856 Kb |

DOI: 10.1007/BF01990348

## Abstract

Because of the special structure of the equations $AX-XB=C$ the usual relation for linear equations ``$\mbox{backward error} = \mbox{relative residual}$'' does not hold, and application of the standard perturbation result for $Ax=b$ yields a perturbation bound involving ${\rm sep}(A,B)^{-1}$ that is not always attainable. An expression is derived for the backward error of an approximate solution $Y$; it shows that the backward error can exceed the relative residual by an arbitrary factor. A sharp perturbation bound is derived and it is shown that the condition number it defines can be arbitrarily smaller than the ${\rm sep}(A,B)^{-1}$-based quantity that is usually used to measure sensitivity. For practical error estimation using the residual of a computed solution an ``LAPACK-style'' bound is shown to be efficiently computable and potentially much smaller than a sep-based bound. A Fortran~77 code has been written that solves the Sylvester equation and computes this bound, making use of LAPACK routines.

Item Type: | Article |
---|---|

Uncontrolled Keywords: | Sylvester equation - Lyapunov equation - backward error - perturbation bound - condition number - error estimate - LAPACK |

Subjects: | MSC 2000 > 15 Linear and multilinear algebra; matrix theory MSC 2000 > 65 Numerical analysis |

MIMS number: | 2006.155 |

Deposited By: | Miss Louise Stait |

Deposited On: | 28 June 2006 |

Download Statistics: last 4 weeks

Repository Staff Only: edit this item