A simple OpenMP scheme for parallel iteration solvers in finite element analysis

Research output: Contribution to journalArticlepeer-review

13 Citations (Scopus)

Abstract

This study develops an OpenMP scheme to parallel the preconditioned conjugate gradient methods (PCG) in shared memory computers. The proposed method is simple and systematic, so a minor change in traditional PCG methods may produce effective parallelism. At first, the global stiffness matrix is renumbered in order to produce a parallel three-line form matrix, and a subroutine only needs to be called once in the finite element analysis. Several basic OpenMP commands are then added into the traditional incomplete Cholesky factorization (ILU) and symmetric successive over-relaxation (SSOR) codes to make the procedures of matrix multiplication, decomposition, forward substitution, and backward substitution fully parallel.

Original languageEnglish
Pages (from-to)91-108
Number of pages18
JournalCMES - Computer Modeling in Engineering and Sciences
Volume64
Issue number1
Publication statusPublished - 2010 Oct 11

All Science Journal Classification (ASJC) codes

  • Software
  • Modelling and Simulation
  • Computer Science Applications

Fingerprint Dive into the research topics of 'A simple OpenMP scheme for parallel iteration solvers in finite element analysis'. Together they form a unique fingerprint.

Cite this