MPI parallel implementation of CBF preconditioning for 3D elasticity problems1This research has been supported by the Copernicus CP 940820 (HIPERGEOS) Project and by the Bulgarian NSF Grant I-701/97.1

The numerical solution of 3D linear elasticity equations is considered. The problem is described by a coupled system of second-order elliptic partial differential equations. This system is discretized by trilinear brick finite elements. The PCG iterative method is used for solving the large-scale linear algebraic systems arising after the FEM discretization of the problem. 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 for preconditioning of the obtained block-diagonal matrix.