Publication at Faculty of Pharmacy in Hradec Králové |
2014
Abstract
This paper deals with estimating the condition number of triangular matrices in the Euclidean norm. The two main incremental methods, based on the work of Bischof and the later work of Duff and Vomel, are compared.
The paper presents new theoretical results revealing their similarities and differences. As typical in condition number estimation, there is no universal always-winning strategy, but theoretical and experimental arguments show that the clearly preferable approach is the algorithm of Duff and Vomel when appropriately applied to both the triangular matrix itself and its inverse.
This leads to a highly accurate incremental condition number estimator.