Parallel DD-MIC(0) Preconditioning of Nonconforming Rotated Trilinear FEM Elasticity Systems

  • Authors:
  • Yavor Vutov

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

  • Venue:
  • Large-Scale Scientific Computing
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

A new parallel preconditioning algorithm for 3D nonconforming FEM elasticity systems is presented. The preconditioner is constructed in two steps. First, displacement decomposition of the stiffness matrix is used. Then MIC(0) factorization is applied to a proper auxiliary M-matrix to get an approximate factorization of the obtained block-diagonal matrix. The auxiliary matrix has a special block structure -- its diagonal blocks are diagonal matrices themselves. This allows the solution of the preconditioning system to be performed efficiently in parallel. Estimates for the parallel times, speedups and efficiencies are derived. The performed parallel tests are in total agreement with them. The robustness of the proposed algorithm is confirmed by the presented experiments solving problems with strong coefficient jumps.