A dimension-oblivious domain decomposition method based on space-filling curves

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Griebel, Michael, Schweitzer, Marc Alexander, Troska, Lukas
Format: Preprint
Published: 2021
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866914518368518144
author Griebel, Michael
Schweitzer, Marc Alexander
Troska, Lukas
author_facet Griebel, Michael
Schweitzer, Marc Alexander
Troska, Lukas
contents In this paper we present an algebraic dimension-oblivious two-level domain decomposition solver for discretizations of elliptic partial differential equations. The proposed parallel solver is based on a space-filling curve partitioning approach that is applicable to any discretization, i.e. it directly operates on the assembled matrix equations. Moreover, it allows for the effective use of arbitrary processor numbers independent of the dimension of the underlying partial differential equation while maintaining optimal convergence behavior. This is the core property required to attain a sparse grid based combination method with extreme scalability which can utilize exascale parallel systems efficiently. Moreover, this approach provides a basis for the development of a fault-tolerant solver for the numerical treatment of high-dimensional problems. To achieve the required data redundancy we are therefore concerned with large overlaps of our domain decomposition which we construct via space-filling curves. In this paper, we propose our space-filling curve based domain decomposition solver and present its convergence properties and scaling behavior. The results of numerical experiments clearly show that our approach provides optimal convergence and scaling behavior in arbitrary dimension utilizing arbitrary processor numbers.
format Preprint
id arxiv_https___arxiv_org_abs_2110_11211
institution arXiv
publishDate 2021
record_format arxiv
spellingShingle A dimension-oblivious domain decomposition method based on space-filling curves
Griebel, Michael
Schweitzer, Marc Alexander
Troska, Lukas
Numerical Analysis
65F50, 65Y05, 65M55, 65N55
In this paper we present an algebraic dimension-oblivious two-level domain decomposition solver for discretizations of elliptic partial differential equations. The proposed parallel solver is based on a space-filling curve partitioning approach that is applicable to any discretization, i.e. it directly operates on the assembled matrix equations. Moreover, it allows for the effective use of arbitrary processor numbers independent of the dimension of the underlying partial differential equation while maintaining optimal convergence behavior. This is the core property required to attain a sparse grid based combination method with extreme scalability which can utilize exascale parallel systems efficiently. Moreover, this approach provides a basis for the development of a fault-tolerant solver for the numerical treatment of high-dimensional problems. To achieve the required data redundancy we are therefore concerned with large overlaps of our domain decomposition which we construct via space-filling curves. In this paper, we propose our space-filling curve based domain decomposition solver and present its convergence properties and scaling behavior. The results of numerical experiments clearly show that our approach provides optimal convergence and scaling behavior in arbitrary dimension utilizing arbitrary processor numbers.
title A dimension-oblivious domain decomposition method based on space-filling curves
topic Numerical Analysis
65F50, 65Y05, 65M55, 65N55
url https://arxiv.org/abs/2110.11211