HIGH-PERFORMANCE PCG SOLVERS FOR FEM STRUCTURAL ANALYSIS

The preconditioned conjugate gradient algorithm is a well-known and powerful method used to solve large sparse symmetric positive definite linear systems. Such systems are generated by the finite element discretization in structural analysis but users of finite elements in this context generally still rely on direct methods. It is our purpose in the present work to highlight the improvement brought forward by some new preconditioning techniques and show that the preconditioned conjugate gradient method performs better than efficient direct methods.