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. |