Parallel incomplete factorization of 3D NC FEM elliptic systems

  • Authors:
  • Yavor Vutov

  • Affiliations:
  • Institute for Parallel Processing, Bulgarian Academy of Sciences

  • Venue:
  • NMA'06 Proceedings of the 6th international conference on Numerical methods and applications
  • Year:
  • 2006

Quantified Score

Hi-index 0.00

Visualization

Abstract

A new parallel preconditioner for solution of large scale second order 3D FEM elliptic systems is presented. The problem is discretized by rotated trilinear non-conforming finite elements. The algorithm is based on application of modified incomplete Cholesky factorisation (MIC(0)) to a locally constructed modification B of the original stiffness matrix A. The matrix B preserves the robustness of the point-wise factorisation and has a special block structure allowing parallelization. The performed numerical tests are in agreement with the derived estimates for the parallel times.