A Simple OpenMP Scheme for Parallel Iteration Solvers in Finite Element Analysis
نویسنده
چکیده
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.
منابع مشابه
OpenMP / MPI Hybrid vs. Flat MPI on the Earth Simulator: Parallel Iterative Solvers for Finite Element Method
An efficient parallel iterative method for finite element method has been developed for symmetric multiprocessor (SMP) cluster architectures with vector processors such as the Earth Simulator. The method is based on a three-level hybrid parallel programming model, including message passing for inter-SMP node communication, loop directives by OpenMP for intra-SMP node parallelization and vectori...
متن کاملRuntime affinity optimization for hybrid MPI + OpenMP solvers
In the physical sciences, it is not uncommon to encounter parallel numerical codes in which computation is partitioned between solution of multiple fundamentally different problems, related in a pipelined manner through their state variables. Such a multi-phase scheme might prove advantageous by reducing the “spin up” and “wind down” time that would otherwise be incurred when using separate cod...
متن کاملParallel Iterative Solvers for Ill-Conditioned Problems with Reordering
1. Preconditioned Iterative Solvers with Multicoloring In the previous work [1], author developed an efficient parallel iterative solver for finite-element applications on the Earth Simulator (ES) [2] using multi-level hybrid parallel programming model with MPI and OpenMP. The method employs three-level hybrid parallel programming model for SMP cluster architectures, consisting of MPI, OpenMP a...
متن کاملViscous Models Comparison in Water Impact of Twin 2D Falling Wedges Simulation by Different Numerical Solvers
In this paper, symmetric water entry of twin wedges is investigated for deadrise angle of 30 degree. Three numerical simulation of a symmetric impact, considering rigid body dynamic equations of motion in two-phase flow is presented. The two-phase flow around the wedges is solved by Finite Element based on Finite Volume method (FEM-FVM) which is used in conjunction with Volume of Fluid (VOF) sc...
متن کاملFEMT, Open Source Tools for Solving Large Systems of Equations in Parallel
We present a new open source library and tools for solving large linear systems of equations with sparse matrices resulting from simulations with finite element, finite volume and finite differences. FEMT is a multi-platform software (Windows, GNU/Linux, Mac OS), released as open source (GNU LGPL). FEMT can run from laptops up to clusters of computers. It has been programmed in modern standard ...
متن کامل