Optimization of two-level methods for DG discretizations of reaction-diffusion equations

by   José Pablo Lucero Lorca, et al.

We analyze and optimize two-level methods applied to a symmetric interior penalty discontinuous Galerkin finite element discretization of a singularly perturbed reaction-diffusion equation. Previous analyses of such methods have been performed by Hemker et. al. for the Poisson problem focusing on optimizing the smoother. Our main innovation is that we optimize the complete two-level process, and we obtain explicit formulas for the optimal relaxation parameter of the two-level method for the Poisson problem in 1D, and closed form approximation formulas for the optimal choice in the reaction-diffusion case in all regimes. Our analysis shows that for DG penalization parameter values used in practice, it is better to use cell block-Jacobi smoothers of Schwarz type, in contrast to earlier results suggesting that point block-Jacobi smoothers are preferable, based on a smoothing analysis alone. Our analysis also reveals how the performance of the iterative solver depends on the DG penalization parameter, and what value should be chosen to get the fastest iterative solver, providing a new, direct link between DG discretization and iterative solver performance. We illustrate our analysis with numerical experiments.



There are no comments yet.


page 1

page 2

page 3

page 4


Multilevel Schwarz preconditioners for singularly perturbed symmetric reaction-diffusion systems

We present robust and highly parallel multilevel non-overlapping Schwarz...

Efficient discretization and preconditioning of the singularly perturbed Reaction Diffusion problem

We consider the reaction diffusion problem and present efficient ways to...

An Assessment of Solvers for Algebraically Stabilized Discretizations of Convection-Diffusion-Reaction Equations

We consider flux-corrected finite element discretizations of 3D convecti...

First-order system least squares finite-elements for singularly perturbed reaction-diffusion equations

We propose a new first-order-system least squares (FOSLS) finite-element...

Monolithic multigrid for a reduced-quadrature discretization of poroelasticity

Advanced finite-element discretizations and preconditioners for models o...

Analysis of a new implicit solver for a semiconductor model

We present and analyze a new iterative solver for implicit discretizatio...

Closed form optimized transmission conditions for complex diffusion with many subdomains

Optimized transmission conditions in domain decomposition methods have b...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.