2nd September 2016
HSL_MA97 uses a direct method to solve large sparse symmetric linear systems of equations . This package optionally uses OpenMP and is designed to achieve bit-compatible results regardless of the number of threads used. It computes the sparse factorization
where (real symmetric or complex symmetric) or (complex Hermitian, where denotes the conjugate transpose of ), is a permutation matrix, is lower triangular, and is block diagonal with blocks of size and .
If is large, bit-compatible answers are not required, we recommend using HSL_MA86 or HSL_MA87 instead for better eﬃciency.