Favati P., Leoncini M., Martinez A. On the robustness of Gaussian Elimination with Partial Pivoting. Technical report, 1999. |

Abstract (English) |
It has been recently shown that large growth factors might occur in Gaussian Elimination with Partial Pivoting (GEPP) also when solving some plausibly {em natural} systems. In this note we argue that this potential problem could be easily solved, with much smaller risk of failure, by very small (and low cost) modifications of the basic algorithm, thus confirming its inherent robustness. To this end, we first propose an informal model with the goal of providing further support to the comprehension of the stability properties of GEPP. We then report the results of numerical experiments that confirm the viewpoint embedded in the model. Basing on the previous observations, we finally propose a simple scheme that could be turned into (even more) accurate software for the solution of linear systems. | |

Subject | G.1.3 Numerical Linear Algebra 65F05 Direct methods for linear systems and matrix inversion |

Open access Restricted Private