Iterative solution methods
MPI parallel implementation of CBF preconditioning for 3D elasticity problems
Mathematics and Computers in Simulation - Special issue from IMACS sponsored conference: “Modelling '98”
MPI: The Complete Reference
MPI solver for 3D elasticity problems
Mathematics and Computers in Simulation - MODELLING 2001 - Second IMACS conference on mathematical modelling and computational methods in mechanics, physics, biomechanics and geodynamics
| Hi-index | 0.00 |
In this paper we consider numerical solution of 3D linear elasticity equations described by a coupled system of second order elliptic partial differential equations. This system is discretized by trilinear parallelepipedal finite elements. Preconditioned Conjugate Gradient iterative method is used for solving large-scale linear algebraic systems arising after the Finite Element Method (FEM) discretization of the problem. The displacement decomposition technique is applied at the first step to construct a preconditioner using the decoupled block diagonal part of the original matrix. Then circulant block factorization is used to precondition thus obtained block diagonal matrix. Since both preconditioning techniques, displacement decomposition and circulant block factorization, are highly parallelizable, a portable parallel FEM code utilizing MPI for communication is implemented. Results of numerical tests performed on a number of modern parallel computers using real life engineering problems from the geosciences (geomechanics in particular) are reported and discussed.