Wilkinson's inertia‐revealing factorization and its application to sparse matrices

Publisher: John Wiley & Sons Inc

E-ISSN: 1099-1506|25|2|nla.2130-nla.2130

ISSN: 1070-5325

Source: NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS (ELECTRONIC), Vol.25, Iss.2, 2018-03, pp. : n/a-n/a

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

SummaryWe propose a new inertia‐revealing factorization for sparse symmetric matrices. The factorization scheme and the method for extracting the inertia from it were proposed in the 1960s for dense, banded, or tridiagonal matrices, but they have been abandoned in favor of faster methods. We show that this scheme can be applied to any sparse symmetric matrix and that the fill in the factorization is bounded by the fill in the sparse QR factorization of the same matrix (but is usually much smaller). We describe our serial proof‐of‐concept implementation and present experimental results, studying the method's numerical stability and performance.