On computing inverse entries of a sparse matrix in an out-of-core environment

Abstract : The inverse of an irreducible sparse matrix is structurally full, so that it is impractical to think of computing or storing it. However, there are several applications where a subset of the entries of the inverse is required. Given a factorization of the sparse matrix held in out-of-core storage, we show how to compute such a subset e ciently, by accessing only parts of the factors. When there are many inverse entries to compute, we need to guarantee that the overall computation scheme has reasonable memory requirements, while minimizing the cost of loading the factors. This leads to a partitioning problem that we prove is NP-complete. We also show that we cannot get a close approximation to the optimal solution in polynomial time. We thus need to develop heuristic algorithms, and we propose: (i) a lower bound on the cost of an optimum solution; (ii) an exact algorithm for a particular case; (iii) two other heuristics for a more general case; and (iv) hypergraph partitioning models for the most general setting. We illustrate the performance of our algorithms in practice using the MUMPS software package on a set of real-life problems as well as some standard test matrices. We show that our techniques can improve the execution time by a factor of 50. Key words. Sparse matrices, direct methods for linear systems and matrix inversion, multifrontal method, graphs and hypergraphs.
Type de document :
Article dans une revue
SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (4), pp.A1975--A1999
Liste complète des métadonnées

https://hal.inria.fr/hal-00763556
Contributeur : Equipe Roma <>
Soumis le : mardi 11 décembre 2012 - 10:01:43
Dernière modification le : mercredi 12 septembre 2018 - 17:46:04

Identifiants

  • HAL Id : hal-00763556, version 1

Citation

Patrick R. Amestoy, Iain S. Duff, Jean-Yves L'Excellent, Yves Robert, François-Henry Rouet, et al.. On computing inverse entries of a sparse matrix in an out-of-core environment. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (4), pp.A1975--A1999. 〈hal-00763556〉

Partager

Métriques

Consultations de la notice

226