An hybrid approach for the parallelization of a block iterative algorithm

  • Authors:
  • Carlos Balsa;Ronan Guivarch;Daniel Ruiz;Mohamed Zenadi

  • Affiliations:
  • CEsA, FEUP, Porto, Portugal;Université de Toulouse, INP, ENSEEIHT, IRIT;Université de Toulouse, INP, ENSEEIHT, IRIT;Université de Toulouse, INP, ENSEEIHT, IRIT

  • Venue:
  • VECPAR'10 Proceedings of the 9th international conference on High performance computing for computational science
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

The Cimmino method is a row projection method in which the original linear system is divided into subsystems. At every iteration, it computes one projection per subsystem and uses these projections to construct an approximation to the solution of the linear system. The usual parallelization strategy in block algorithms is to distribute the different blocks on the available processors. In this paper, we follow another approach where we do not perform explicitly this block distribution to processors within the code, but let the multi-frontal sparse solver MUMPS handle the data distribution and parallelism. The data coming from the subsystems defined by the block partition in the Block Cimmino method are gathered in an unique block diagonal sparse matrix which is analysed, distributed and factorized in parallel by MUMPS. Our target is to define a methodology for parallelism based only on the functionalities provided by general sparse solver libraries and how efficient this way of doing can be.