MPI parallel implementation of CBF preconditioning for 3D elasticity problems Original Research Article

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.