A Comparison of Algorithms for Solving Symmetric Indefinite Systems of Linear Equations

A brief description of numerically stable methods for solving symmetric indefinite systems of linear equations Ax = b is presented. The computer implementation of the methods is described and some numerical experiments are reported, comparing program performance on two different computers. The comparison leads to the conclusion that the algorithms proposed by Aasen and Bunch-Kaufman appear to be substantially superior to alternate methods.