Towards an efficient tile matrix inversion of symmetric positive definite matrices on multicore architectures

Emmanuel Agullo (1, 2, 3) , Henricus Bouwmeester (4) , Jack Dongarra (3) , Jakub Kurzak (3) , Julien Langou (4) , Lee Rosenberg (4)

1 HiePACS - High-End Parallel Algorithms for Challenging Numerical Simulations
2 LaBRI - Laboratoire Bordelais de Recherche en Informatique
3 ICL - Innovative Computing Laboratory [Knoxville]
4 Department of Mathematical and Statistical Sciences

The algorithms in the current sequential numerical linear algebra libraries (e.g. LAPACK) do not parallelize well on multicore architectures. A new family of algorithms, the tile algorithms, has recently been introduced. Previous research has shown that it is possible to write efficient and scalable tile algorithms for performing a Cholesky factorization, a (pseudo) LU factorization, a QR factorization, and computing the inverse of a symmetric positive denite matrix. In this extended abstract, we revisit the computation of the inverse of a symmetric positive denite matrix. We observe that, using a dynamic task scheduler, it is relatively painless to translate existing LAPACK code to obtain a ready-to-be-executed tile algorithm. However we demonstrate that, for some variants, non trivial compiler techniques (array renaming, loop reversal and pipelining) need then to be applied to further increase the parallelism of the application. We present preliminary experimental results.

Domains: Distributed, Parallel, and Cluster Computing [cs.DC]

Emmanuel Agullo, Henricus Bouwmeester, Jack Dongarra, Jakub Kurzak, Julien Langou, et al.. Towards an efficient tile matrix inversion of symmetric positive definite matrices on multicore architectures. 9th International Meeting on High Performance Computing for Computational Science (VecPar'10), Jun 2010, Berkeley, United States. ⟨inria-00548906⟩