We introduce a new code for the direct solution of sparse symmetric linear equations that solves indefinite systems with 2x2 pivoting for stability. The code, called MA57, is in HSL 2002 and supersedes the well used HSL code MA27. We describe the user interface in some detail and emphasize some of the novel features of MA57. These include restart facilities, matrix modification, partial solution for matrix factors, solution of multiple right-hand sides, and iterative refinement and error analysis.
展开▼