5,255 research outputs found
Efficient Multigrid Preconditioners for Atmospheric Flow Simulations at High Aspect Ratio
Many problems in fluid modelling require the efficient solution of highly
anisotropic elliptic partial differential equations (PDEs) in "flat" domains.
For example, in numerical weather- and climate-prediction an elliptic PDE for
the pressure correction has to be solved at every time step in a thin spherical
shell representing the global atmosphere. This elliptic solve can be one of the
computationally most demanding components in semi-implicit semi-Lagrangian time
stepping methods which are very popular as they allow for larger model time
steps and better overall performance. With increasing model resolution,
algorithmically efficient and scalable algorithms are essential to run the code
under tight operational time constraints. We discuss the theory and practical
application of bespoke geometric multigrid preconditioners for equations of
this type. The algorithms deal with the strong anisotropy in the vertical
direction by using the tensor-product approach originally analysed by B\"{o}rm
and Hiptmair [Numer. Algorithms, 26/3 (2001), pp. 219-234]. We extend the
analysis to three dimensions under slightly weakened assumptions, and
numerically demonstrate its efficiency for the solution of the elliptic PDE for
the global pressure correction in atmospheric forecast models. For this we
compare the performance of different multigrid preconditioners on a
tensor-product grid with a semi-structured and quasi-uniform horizontal mesh
and a one dimensional vertical grid. The code is implemented in the Distributed
and Unified Numerics Environment (DUNE), which provides an easy-to-use and
scalable environment for algorithms operating on tensor-product grids. Parallel
scalability of our solvers on up to 20,480 cores is demonstrated on the HECToR
supercomputer.Comment: 22 pages, 6 Figures, 2 Table
Distributed Hierarchical SVD in the Hierarchical Tucker Format
We consider tensors in the Hierarchical Tucker format and suppose the tensor
data to be distributed among several compute nodes. We assume the compute nodes
to be in a one-to-one correspondence with the nodes of the Hierarchical Tucker
format such that connected nodes can communicate with each other. An
appropriate tree structure in the Hierarchical Tucker format then allows for
the parallelization of basic arithmetic operations between tensors with a
parallel runtime which grows like , where is the tensor dimension.
We introduce parallel algorithms for several tensor operations, some of which
can be applied to solve linear equations directly in the
Hierarchical Tucker format using iterative methods like conjugate gradients or
multigrid. We present weak scaling studies, which provide evidence that the
runtime of our algorithms indeed grows like . Furthermore, we present
numerical experiments in which we apply our algorithms to solve a
parameter-dependent diffusion equation in the Hierarchical Tucker format by
means of a multigrid algorithm
On the validity of the local Fourier analysis
Local Fourier analysis (LFA) is a useful tool in predicting the convergence
factors of geometric multigrid methods (GMG). As is well known, on rectangular
domains with periodic boundary conditions this analysis gives the exact
convergence factors of such methods. In this work, using the Fourier method, we
extend these results by proving that such analysis yields the exact convergence
factors for a wider class of problems
- …