HSL_MA97 : a bit-compatible multifrontal code for sparse symmetric systems
JD Hogg, J. A. Scott · Science and Technology Facilities Council · 2011
The multifrontal method is widely used for the numerical solution of large sparse symmetric linear systems of equations. In this report, we discuss the design and implementation of a new multifrontal code HSL_MA97. Our motivation for the new code is discussed along with the key design features. HSL_MA97 is for real and complex problems and is designed to be efficient for both positive definite and indefinite systems. HSL_MA97 can be run in serial or in a shared memory environment using OpenMP. An important feature is that in parallel it offers bit-compatible solutions. Numerical results are presented for a range of problems from practical applications and comparisons are made with existing codes. Future plans for HSL_MA97 are outlined.