Global and localized parallel preconditioning techniques for large scale solid Earth simulations

We investigate and compare a few parallel preconditioning techniques in the iterative solution of large sparse linear systems arising from solid Earth simulation with and without using contact information in the domain partitioning process. Previous studies are focused on using static or matrix pattern based incomplete LU (ILU) preconditioners in a localized preconditioner implementation. Our current studies are concerned with preconditioner performance for solving two different cases with and without known contact information. For the cases with contact information, we use localized threshold value based incomplete LU (ILUT) preconditioner to improve efficiency. For the cases without contact information, we use a global sparse approximate inverse preconditioner with a static sparsity pattern to achieve robustness. Numerical results from simulating ground motion on a parallel supercomputer are given to compare the effectiveness of these parallel preconditioning techniques.