A quantitative performance analysis for Stokes solvers at the extreme scale

by   Björn Gmeiner, et al.

This article presents a systematic quantitative performance analysis for large finite element computations on extreme scale computing systems. Three parallel iterative solvers for the Stokes system, discretized by low order tetrahedral elements, are compared with respect to their numerical efficiency and their scalability running on up to 786 432 parallel threads. A genuine multigrid method for the saddle point system using an Uzawa-type smoother provides the best overall performance with respect to memory consumption and time-to-solution. The largest system solved on a Blue Gene/Q system has more than ten trillion (1.1 · 10 ^13) unknowns and requires about 13 minutes compute time. Despite the matrix free and highly optimized implementation, the memory requirement for the solution vector and the auxiliary vectors is about 200 TByte. Brandt's notion of "textbook multigrid efficiency" is employed to study the algorithmic performance of iterative solvers. A recent extension of this paradigm to "parallel textbook multigrid efficiency" makes it possible to assess also the efficiency of parallel iterative solvers for a given hardware architecture in absolute terms. The efficiency of the method is demonstrated for simulating incompressible fluid flow in a pipe filled with spherical obstacles.


page 16

page 17

page 18


Parallel solution of saddle point systems with nested iterative solvers based on the Golub-Kahan Bidiagonalization

We present a scalability study of Golub-Kahan bidiagonalization for the ...

Accelerating linear solvers for Stokes problems with C++ metaprogramming

The efficient solution of large sparse saddle point systems is very impo...

High-order matrix-free incompressible flow solvers with GPU acceleration and low-order refined preconditioners

We present a matrix-free flow solver for high-order finite element discr...

Textbook efficiency: massively parallel matrix-free multigrid for the Stokes system

We employ textbook multigrid efficiency (TME), as introduced by Achi Bra...

Performance of an iterative wavelet reconstructor for the Multi-conjugate Adaptive Optics RelaY of ESO's ELT

The Multi-conjugate Adaptive Optics RelaY (MAORY) is one of the key Adap...

Analytical Estimation of the Scalability of Iterative Numerical Algorithms on Distributed Memory Multiprocessors

This article presents a new high-level parallel computational model name...

Recent Developments in Iterative Methods for Reducing Synchronization

On modern parallel architectures, the cost of synchronization among proc...