A dissection solver with kernel detection for symmetric finite element matrices on shared memory computers - Sorbonne Université
Preprints, Working Papers, ... Year : 2013

A dissection solver with kernel detection for symmetric finite element matrices on shared memory computers

Abstract

A direct solver for symmetric sparse matrices from finite element problems is presented. The solver is supposed to work as a local solver of domain decomposition methods for hybrid parallelization on cluster systems of multi-core CPUs, then it is required to run on shared memory computers and to have an ability of kernel detection. Symmetric pivoting with a given threshold factorizes a matrix with a decomposition introduced by a nested bisection and selects suspicious null pivots from the threshold. The Schur complement constructed from the suspicious null pivots is examined by a factorization with 1×1 and 2×2 pivoting and by a robust kernel detection algorithm based on measurements of residuals with orthogonal projections onto supposed image spaces. A static data structure from the nested bisection and a block sub- structure for Schur complements on all bisection-levels can well employ level 3 BLAS routines. Asynchronous task execution for each block can reduce idling time of processors drastically, then the solver has high parallel efficiency. Competitive performance of the developed solver to Intel Pardiso on shared memory computers is shown by numerical experiments.
Fichier principal
Vignette du fichier
SuzukiRoux.pdf (1.42 Mo) Télécharger le fichier
Origin Files produced by the author(s)

Dates and versions

hal-00816916 , version 1 (23-04-2013)
hal-00816916 , version 2 (30-10-2013)
hal-00816916 , version 3 (04-04-2014)

Identifiers

  • HAL Id : hal-00816916 , version 1

Cite

Atsushi Suzuki, François-Xavier Roux. A dissection solver with kernel detection for symmetric finite element matrices on shared memory computers. 2013. ⟨hal-00816916v1⟩

Collections

ONERA
796 View
612 Download

Share

More