Scaled Companion Linearizations for Coefficientwise Backward Stability
Abstract
A central goal in solving a polynomial eigenvalue problem is to compute eigenpairs that are exact for a nearby polynomial, with each coefficient perturbation bounded by a modest multiple of unit roundoff times that coefficient's norm. A backward stable eigensolver for a linearization does not automatically provide this coefficientwise relative guarantee. Scaling can improve accuracy, but a single scaling may not work across a widely spread spectrum. We construct a family of scaled companion linearizations that covers the nonzero complex plane. Each linearization is paired with a region and an eigenvector recovery rule for which the backward-error amplification factor is bounded at exact pencil eigenvectors, independently of coefficient norm disparities. Tropical roots guide the construction, and their multiplicities reduce the number of required solves. These results provide a basis for coefficientwise relative backward stability across the spectrum.
Related
- Forward-Backward McKean-Vlasov System with Jumps and a Monotone Source With à Control-Shifted Interface : Well-Posedness, Stability, and Yosida Rates
- Invariant subspaces for free linearizations of Lipschitz maps
- Cokernels of random p-adic orthogonal matrices and their linearizations
- Interlacing on the Unit Circle via Coefficientwise Reciprocals
- Coefficientwise Hankel-total positivity from a universal cumulant lifting
BibTeX
@misc{perez2026scaled,
title = {{Scaled Companion Linearizations for Coefficientwise Backward Stability}},
author = {Javier Pérez},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.20285},
eprint = {2610.20285},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.20285},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Scaled Companion Linearizations for Coefficientwise Backward Stability
AU - Javier Pérez
PY - 2026
DA - 2026/10/14
PB - arXiv
DO - 10.48550/arxiv.2610.20285
UR - https://arxiv.org/abs/2610.20285
AB - A central goal in solving a polynomial eigenvalue problem is to compute eigenpairs that are exact for a nearby polynomial, with each coefficient perturbation bounded by a modest multiple of unit roundoff times that coefficient's norm. A backward stable eigensolver for a linearization does not automatically provide this coefficientwise relative guarantee. Scaling can improve accuracy, but a single scaling may not work across a widely spread spectrum. We construct a family of scaled companion linearizations that covers the nonzero complex plane. Each linearization is paired with a region and an eigenvector recovery rule for which the backward-error amplification factor is bounded at exact pencil eigenvectors, independently of coefficient norm disparities. Tropical roots guide the construction, and their multiplicities reduce the number of required solves. These results provide a basis for coefficientwise relative backward stability across the spectrum.
ER -