A Parallel MIC(0) Preconditioning Technique for Elliptic Problems


Svetozar Margenov and Yavor Vutov
IPP, BAS, Sofia

In this talk we present a parallel preconditioner technique for solution of large scale second order 3D FEM elliptic systems. The problem is discretized by rotated trilinear non-conforming finite elements.
The algorithm is based on application of modified incomplete Cholesky factorization (MIC(0)) to a locally constructed modification B of the original stiffness matrix A. The matrix B preserves the robustness of the point-wise factorization and has a special block structure allowing parallelization. Two types of modification are considered which lead to algorithms with different parallel properties. Results from parallel experiments on several parallel computers are shown.