Researcher profile

Martin J. Gander

Martin J. Gander contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

Closed form optimized transmission conditions for complex diffusion with many subdomains

Optimized transmission conditions in domain decomposition methods have been the focus of intensive research efforts over the past decade. Traditionally, transmission conditions are optimized for two subdomain model configurations, and then used in practice for many subdomains. We optimize here transmission conditions for the first time directly for many subdomains for a class of complex diffusion problems. Our asymptotic analysis leads to closed form optimized transmission conditions for many subdomains, and shows that the asymptotic best choice in the mesh size only differs from the two subdomain best choice in the constants, for which we derive the dependence on the number of subdomains explicitly, including the limiting case of an infinite number of subdomains, leading to new insight into scalability. Our results include both Robin and Ventcell transmission conditions, and we also optimize for the first time a two-sided Ventcell condition. We illustrate our results with numerical experiments, both for situations covered by our analysis and situations that go beyond.

preprint2022arXiv

Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation

We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the subdomain problems satisfy first-order absorbing (impedance) transmission conditions, and exchange of information between subdomains is achieved using a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1,2,3) that it is well-defined in a tensor product of appropriate local function spaces, each with $L^2$ impedance boundary data. Given this, we then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps which ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. These results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.

preprint2022arXiv

Low-rank Parareal: a low-rank parallel-in-time integrator

In this work, the Parareal algorithm is applied to evolution problems that admit good low-rank approximations and for which the dynamical low-rank approximation (DLRA) can be used as time stepper. Many discrete integrators for DLRA have recently been proposed, based on splitting the projected vector field or by applying projected Runge--Kutta methods. The cost and accuracy of these methods are mostly governed by the rank chosen for the approximation. These properties are used in a new method, called low-rank Parareal, in order to obtain a time-parallel DLRA solver for evolution problems. The algorithm is analyzed on affine linear problems and the results are illustrated numerically.

preprint2022arXiv

On the nonlinear Dirichlet-Neumann method and preconditioner for Newton's method

The Dirichlet-Neumann (DN) method has been extensively studied for linear partial differential equations, while little attention has been devoted to the nonlinear case. In this paper, we analyze the DN method both as a nonlinear iterative method and as a preconditioner for Newton's method. We discuss the nilpotent property and prove that under special conditions, there exists a relaxation parameter such that the DN method converges quadratically. We further prove that the convergence of Newton's method preconditioned by the DN method is independent of the relaxation parameter. Our numerical experiments further illustrate the mesh independent convergence of the DN method and compare it with other standard nonlinear preconditioners.

preprint2022arXiv

Schwarz methods by domain truncation

Schwarz methods use a decomposition of the computational domain into subdomains and need to put boundary conditions on the subdomain boundaries. In domain truncation one restricts the unbounded domain to a bounded computational domain and also needs to put boundary conditions on the computational domain boundaries. It turns out to be fruitful to think of the domain decomposition in Schwarz methods as truncation of the domain onto subdomains. The first truly optimal Schwarz method that converges in a finite number of steps was proposed in 1994 and used precisely transparent boundary conditions as transmission conditions between subdomains. Approximating these transparent boundary conditions for fast convergence of Schwarz methods led to the development of optimized Schwarz methods -- a name that has become common for Schwarz methods based on domain truncation. Compared to classical Schwarz methods which use simple Dirichlet transmission conditions and have been successfully used in a wide range of applications, optimized Schwarz methods are much less well understood, mainly due to their more sophisticated transmission conditions. This present situation is the motivation for our survey: to give a comprehensive review and precise exploration of convergence behaviors of optimized Schwarz methods based on Fourier analysis taking into account the original boundary conditions, many subdomain decompositions and layered media. The transmission conditions we study include the lowest order absorbing conditions (Robin), and also more advanced perfectly matched layers (PML), both developed first for domain truncation.

preprint2020arXiv

Asymptotic Analysis for Overlap in Waveform Relaxation Methods for RC Type Circuits

Waveform relaxation (WR) methods are based on partitioning large circuits into sub-circuits which then are solved separately for multiple time steps in so-called time windows, and an iteration is used to converge to the global circuit solution in each time window. Classical WR converges quite slowly, especially when long time windows are used. To overcome this issue, optimized WR (OWR) was introduced which is based on optimized transmission conditions that transfer information between the sub-circuits more efficiently than classical WR. We study here for the first time the influence of overlapping sub-circuits in both WR and OWR applied to RC circuits. We give a circuit interpretation of the new transmission conditions in OWR and derive closed-form asymptotic expressions for the circuit elements representing the optimization parameter in OWR. Our analysis shows that the parameter is quite different in the overlapping case, compared to the nonoverlapping one. We then show numerically that our optimized choice performs well, also for cases not covered by our analysis. This paper provides a general methodology to derive optimized parameters and can be extended to other circuits or system of differential equations or space-time PDEs.