Symmetric permutations for I-Matrices to delay and avoid small pivots during factorization

In this article, we present several new permutations for I-matrices making these more suitable for incomplete LDU-factorization preconditioners used in solving linear sys-tems by iterative methods. A general matrix can be transformed by row permutation as well as row and columns scaling into an I-matrix, i.e. a matrix having elements of modulus 1 on the diagonal and elements of modulus of no more than 1 elsewhere. Re-ordering rows and columns by the same permutation clearly preserves I-matrices. In this article, we consider such reordering techniques which make the permuted matrix more suitable for an incomplete LDU-factorization preconditioner than the original I-matrix. We use a multilevel ILUC, an incomplete LDU-factorization preconditioner using Crout’s implementation of Gaussian elimination without pivoting to test these reorderings. The combination of I-matrix preprocessing with the various algorithms presented here and the multilevel incomplete LDU-factorizations forms a powerful preconditioning method for unsymmetric, highly indefinite problems