Sparse Supernodal Solver Using Hierarchical Compression over Runtime System

Grégoire Pichon (1) , Eric Darve (2, 3) , Mathieu Faverge (4, 1) , Pierre Ramet (4, 1) , Jean Roman (1)

1 HiePACS - High-End Parallel Algorithms for Challenging Numerical Simulations
2 Department of Mechanical Engineering [Stanford]
3 ICME - Institute for Computational and Mathematical Engineering [Stanford]
4 LaBRI - Laboratoire Bordelais de Recherche en Informatique

In this talk, we present the PaStiX sparse supernodal solver, using hierarchical compression to reduce the burden on large blocks appearing during the nested dissection process. We compare the numerical stability, and the performance in terms of memory consumption and time to solution of different approaches by selecting when the compression of the factorized matrix occurs. In order to improve the efficiency of the sparse update kernel for both BLR (block low rank) and HODLR (hierarchically off-diagonal low-rank), we investigate the BDLR (boundary distance low-rank) method to preselect rows and columns in the low-rank approximation algorithm.

Mots clés: sparse direct solver; h-matrix

Domaines: Calcul parallèle, distribué et partagé [cs.DC]

Grégoire Pichon, Eric Darve, Mathieu Faverge, Pierre Ramet, Jean Roman. Sparse Supernodal Solver Using Hierarchical Compression over Runtime System. SIAM Conference on Computation Science and Engineering (CSE'17), Feb 2017, Atlanta, United States. ⟨hal-01421379⟩