Parallelisation of BICG-STAB algorithm for finite element computations

Date

2003-09-03

Department

Program

Citation of Original Publication

Bedřich Sousedík, Jaroslav Novotný: Parallelization of the BICG-STAB algorithm for finite element computations, Proceedings of Software and Algorithms of Numerical Mathematics (15.), Hejnice, Czech Republic, September 8-12, 2003, Institute of Thermomechanics AS CR

Rights

This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.

Subjects

Abstract

We describe parallelization and implementation of BICG-STAB method for systems of unsymmetric linear equations arising in nite element discretization. Global matrix is not constructed, we use element by element approach. All element matrices together with addressing vectors are distributed among processors and read into memory. Then multiplication of element matrices by a global vector is performed simultaneously on different processors. Two odes were developed: using MPI library for distributed memory systems and using OpenMP for shared memory systems. Good scalability is demonstrated on an example with about 30000 unknowns.