## 2009.51: A Framework for Analyzing Nonlinear Eigenproblems and Parametrized Linear Systems

2009.51:
Laurence Grammont, Nicholas J. Higham and Françoise Tisseur
(2011)
*A Framework for Analyzing Nonlinear Eigenproblems
and Parametrized Linear Systems.*
Linear Algebra and its Applications, 435 (3).
pp. 623-640.
ISSN 0024-3795

*This is the latest version of this eprint.*

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 264 Kb |

DOI: 10.1016/j.laa.2009.12.038

## Abstract

Associated with an $n\times n$ matrix polynomial of degree $\ell$, $P(\lambda) = \sum_{j=0}^\ell \lambda^j A_j$, are the eigenvalue problem $P(\lambda)x = 0$ and the linear system problem $P(\omega)x = b$, where in the latter case $x$ is to be computed for many values of the parameter $\omega$. Both problems can be solved by conversion to an equivalent problem $L(\lambda)z = 0$ or $L(\omega)z = c$ that is linear in the parameter $\lambda$ or $\omega$. This linearization process has received much attention in recent years for the eigenvalue problem, but it is less well understood for the linear system problem. We develop a framework in which more general versions of both problems can be analyzed, based on one-sided factorizations connecting a general nonlinear matrix function $N(\lambda)$ to a simpler function $M(\lambda)$, typically a polynomial of degree 1 or 2. Our analysis relates the solutions of the original and lower degree problems and in the linear system case indicates how to choose the right-hand side $c$ and recover the solution $x$ from $z$. For the eigenvalue problem this framework includes many special cases studied in the literature, including the vector spaces of pencils $\mathbb{L}_1(P)$ and $\mathbb{L}_2(P)$ recently introduced by Mackey, Mackey, Mehl, and Mehrmann and a class of rational problems. We use the framework to investigate the conditioning and stability of the parametrized linear system $P(\omega)x = b$ and thereby study the effect of scaling, both of the original polynomial and of the pencil $L$. Our results identify situations in which scaling can potentially greatly improve the conditioning and stability and our numerical results show that dramatic improvements can be achieved in practice.

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

Uncontrolled Keywords: | nonlinear eigenvalue problem, polynomial eigenvalue problem, rational eigenvalue problem, linearization, quadratization, parametrized linear system, backward error, scaling, companion form, CICADA |

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

MIMS number: | 2009.51 |

Deposited By: | Nick Higham |

Deposited On: | 06 May 2011 |

### Available Versions of this Item

- A Framework for Analyzing Nonlinear Eigenproblems
and Parametrized Linear Systems (deposited 06 May 2011)
**[Currently Displayed]**

Download Statistics: last 4 weeks

Repository Staff Only: edit this item