Parallel Performance of an Iterative Solver Based on the Golub-Kahan Bidiagonalization

We present an iterative method based on a generalization of the Golub-Kahan bidiagonalization for solving indefinite matrices with a 2 \(\times \) 2 block structure. We focus in particular on our recent implementation of the algorithm using the parallel numerical library PETSc. Since the algorithm is a nested solver, we investigate different choices for parallel inner solvers and show its strong scalability for two Stokes test problems. The algorithm is found to be scalable for large sparse problems.

[1]  Dominique Orban,et al.  Iterative Solution of Symmetric Quasi-Definite Linear Systems , 2017 .

[2]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[3]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[4]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[5]  M. Saunders.  Computing projections with LSQR , 1997 .

[6]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[7]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[8]  Barbara I. Wohlmuth,et al.  A quantitative performance study for Stokes solvers at the extreme scale , 2016, J. Comput. Sci..

[9]  Gene H. Golub,et al.  Calculating the singular values and pseudo-inverse of a matrix , 2007, Milestones in Matrix Computation.

[10]  Gene H. Golub,et al.  On Solving Block-Structured Indefinite Linear Systems , 2003, SIAM J. Sci. Comput..

[11]  S. MacLachlan,et al.  On iterative methods for the incompressible Stokes problem , 2011 .

[12]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[13]  Andrew J. Wathen,et al.  Performance and analysis of saddle point preconditioners for the discrete steady-state Navier-Stokes equations , 2002, Numerische Mathematik.

[14]  Christiaan M. Klaij.  On the stabilization of finite volume methods with co-located variables for incompressible flow , 2015, J. Comput. Phys..

[15]  Achi Brandt,et al.  Multigrid Techniques: 1984 Guide with Applications to Fluid Dynamics, Revised Edition , 2011 .

[16]  Ulrich Rüde,et al.  An iterative generalized Golub-Kahan algorithm for problems in structural mechanics , 2018, ArXiv.

[17]  Mario Arioli,et al.  Generalized Golub-Kahan Bidiagonalization and Stopping Criteria , 2013, SIAM J. Matrix Anal. Appl..

[18]  Jack Dongarra,et al.  Applied Mathematics Research for Exascale Computing , 2014 .

[19]  M. Saunders.  Solution of sparse rectangular systems using LSQR and CRAIG , 1995 .

[20]  Ulrich Rüde,et al.  Towards Textbook Efficiency for Parallel Multigrid , 2015 .

[21]  Howard C. Elman,et al.  Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics , 2014 .