Quasi-matrix-free hybrid multigrid on dynamically adaptive Cartesian grids

07/03/2016
by   Marion Weinzierl, et al.
0

We present a family of spacetree-based multigrid realizations using the tree's multiscale nature to derive coarse grids. They align with matrix-free geometric multigrid solvers as they never assemble the system matrices which is cumbersome for dynamically adaptive grids and full multigrid. The most sophisticated realizations use BoxMG to construct operator-dependent prolongation and restriction in combination with Galerkin/Petrov-Galerkin coarse-grid operators. This yields robust solvers for nontrivial elliptic problems. We embed the algebraic, problem- and grid-dependent multigrid operators as stencils into the grid and evaluate all matrix-vector products in-situ throughout the grid traversals. While such an approach is not literally matrix-free---the grid carries the matrix---we propose to switch to a hierarchical representation of all operators. Only differences of algebraic operators to their geometric counterparts are held. These hierarchical differences can be stored and exchanged with small memory footprint. Our realizations support arbitrary dynamically adaptive grids while they vertically integrate the multilevel operations through spacetree linearization. This yields good memory access characteristics, while standard colouring of mesh entities with domain decomposition allows us to use parallel manycore clusters. All realization ingredients are detailed such that they can be used by other codes.

READ FULL TEXT

page 4

page 19

page 24

page 29

research
03/25/2019

Dynamically Adaptive FAS for an Additively Damped AFAC Variant

Multigrid solvers face multiple challenges on parallel computers. Two fu...
research
08/10/2020

Multigrid with Nonstandard Coarsening

We consider the numerical solution of Poisson's equation on structured g...
research
08/03/2023

Fundamental Data Structures for Matrix-Free Finite Elements on Hybrid Tetrahedral Grids

This paper presents efficient data structures for the implementation of ...
research
10/27/2022

A matrix-free ILU realization based on surrogates

Matrix-free techniques play an increasingly important role in large-scal...
research
07/15/2023

Reducing operator complexity in Algebraic Multigrid with Machine Learning Approaches

We propose a data-driven and machine-learning-based approach to compute ...
research
08/11/2023

A matrix-free parallel two-level deflation preconditioner for the two-dimensional Helmholtz problems

We propose a matrix-free parallel two-level-deflation preconditioner com...
research
01/11/2023

Study of Mach reflection in inviscid flows

In this paper, we study the Mach reflection phenomenon in inviscid flows...

Please sign up or login with your details

Forgot password? Click here to reset