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.
|Number of pages||18|
|Journal||CMES - Computer Modeling in Engineering and Sciences|
|Publication status||Published - 2010 Oct 11|
All Science Journal Classification (ASJC) codes
- Modelling and Simulation
- Computer Science Applications