Parallel MIC(0) preconditioning of 3D elliptic problems discretized by Rannacher-Turek finite elements

  • Authors:
  • P. Arbenz;S. Margenov;Y. Vutov

  • Affiliations:
  • Institute of Computational Science, ETH Zurich, Universitätstrasse 6, 8092 Zurich, Switzerland;Institute for Parallel Processing, Bulgarian Academy of Sciences, Acad. G. Bonchev Bl. 25A, 1113 Sofia, Bulgaria;Institute for Parallel Processing, Bulgarian Academy of Sciences, Acad. G. Bonchev Bl. 25A, 1113 Sofia, Bulgaria

  • Venue:
  • Computers & Mathematics with Applications
  • Year:
  • 2008

Quantified Score

Hi-index 0.09

Visualization

Abstract

Novel parallel algorithms for the solution of large FEM linear systems arising from second order elliptic partial differential equations in 3D are presented. The problem is discretized by rotated trilinear nonconforming Rannacher-Turek finite elements. The resulting symmetric positive definite system of equations Ax=f is solved by the preconditioned conjugate gradient algorithm. The preconditioners employed are obtained by the modified incomplete Cholesky factorization MIC(0) of two kinds of auxiliary matrices B that both are constructed as locally optimal approximations of A in the class of M-matrices. Uniform estimates for the condition number @k(B^-^1A) are derived. Two parallel algorithms based on the different block structures of the related matrices B are studied. The numerical tests confirm theory in that the algorithm scales as O(N^7^/^6) in the matrix order N.