New parallel symmetric SOR preconditioners by multi-type partitioning

A new parallel symmetric successive over-relaxation (PSSOR) preconditioner is proposed in this paper by the multi-type partition techniques introduced in SIAM J. Scientific Computing 20, 2006, pp. 1513–1533. In a general matrix expression, it is proved to be symmetric and positive-definite (SPD) if the coefficient matrix of a linear system is SPD. It is also proved to be equivalent to the SSOR preconditioner using the multi-type ordering. Thus, it works for the preconditioned conjugate gradient method (PCG) and can be analysed by the classic SOR theory. Numerical tests on an anisotropic model problem show that the PSSOR preconditioner can make the PCG have a faster rate of convergence and better parallel performance than the red–black SSOR preconditioner. They also confirm that the PSSOR preconditioner can have a rate of convergence that is nearly the same as the classic sequential SSOR preconditioner when the problem has large anisotropy.

[1]  Loyce Adams,et al.  New Parallel SOR Method by Domain Partitioning , 1999, SIAM J. Sci. Comput..

[2]  R. Grimes,et al.  On vectorizing incomplete factorization and SSOR preconditioners , 1988 .

[3]  Dexuan Xie,et al.  A New Block Parallel SOR Method and Its Analysis , 2005, SIAM J. Sci. Comput..

[4]  S. A. STOTLAND,et al.  Orderings for Parallel Conjugate Gradient Preconditioners , 1997, SIAM J. Sci. Comput..

[5]  Anthony Skjellum,et al.  Using MPI: portable parallel programming with the message-passing interface, 2nd Edition , 1999, Scientific and engineering computation series.

[6]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[7]  William Gropp,et al.  Skjellum using mpi: portable parallel programming with the message-passing interface , 1994 .

[8]  Th. Lippert,et al.  A parallel SSOR preconditioner for lattice QCD , 1997 .

[9]  Andreas Frommer,et al.  Block colouring schemes for the SOR method on local memory parallel computers , 1990, Parallel Comput..

[10]  Charbel Farhat.  Multiprocessors in computational mechanics , 1987 .

[11]  Th. Lippert,et al.  A parallel SSOR preconditioner for lattice QCD , 1996, hep-lat/9608066.

[12]  James M. Ortega,et al.  Multicoloring with lots of colors , 1989, ICS '89.

[13]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[14]  Henk A. van der Vorst,et al.  Large tridiagonal and block tridiagonal linear systems on vector and parallel computers , 1987, Parallel Comput..

[15]  Thomas Lippert,et al.  Parallel SSOR preconditioning for lattice QCD , 1999, Parallel Comput..

[16]  Gene H. Golub,et al.  Matrix computations , 1983 .

[17]  Rami Melhem,et al.  Toward Efficient Implementation of Preconditioned Conjugate Gradient Methods On Vector Supercomputers , 1987 .

[18]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .