We review the derivation of weakly penalized discontinuous Galerkin methods for scattering dominated radiation transport and extend the asymptotic analysis to non-isotropic scattering.
We present robust and highly parallel multilevel non-overlapping Schwarz preconditioners, to solve an interior penalty discontinuous Galerkin finite element discretization of a system of steady state, singularly perturbed reaction-diffusion equations with a singular reaction operator, using a GMRES solver. We provide proofs of convergence for the two-level setting and the multigrid V-cycle as well as numerical results for a wide range of regimes.
We consider a discontinuous interpolation operator with a parameter that controls the discontinuity, and determine the optimal choice for the discontinuity, leading to the fastest solver for a specific 1D symmetric interior penalty DG discretization model problem. We show in addition that our optimization delivers a perfectly clustered spectrum with a high geometric multiplicity, which is very advantageous for a Krylov solver using the method as its preconditioner.
We present an efficient solver for the linear system that arises from the Hierarchical Poincaré-Steklov (HPS) discretization of three dimensional variable coefficient Helmholtz problems. This work is the first efficient iterative solver for the linear system that results from the HPS discretization. Numerical results show the capacity to tackle problems of approximately 100 wavelengths in each direction, requiring more than a billion unknowns to achieve approximately 4 digits of accuracy taking less than 20 minutes to solve.
We apply two-level methods to a symmetric interior penalty discontinuous Galerkin finite element discretization of a singularly perturbed reactions-diffusion equation. The main innovation is that explicit formulas for the optimal relaxation parameter of the two-level method for the Poisson problem in 1D are obtained, as well as very accurate closed form approximation formulas for the optimal choice in the reaction-diffusion case in all regimes. Numerical experiments and comparisons show the applicability of the expressions obtained in higher dimensions and general geometries.
We revisit the Hierarchical Poincaré–Steklov (HPS) method for the Poisson equation using standard Q1 finite elements, building on the original in [1]. While corner degrees of freedom were implicitly handled in that work, subsequent spectral-element implementations have typically avoided them. In Q1-FEM, however, corner coupling cannot be factored out, and we show how the HPS merge procedure naturally accommodates it when corners are enclosed by elements. This clarification bridges a conceptual gap between algebraic Schur-complement methods and operator-based formulations, providing a consistent path for the FEM community to adopt HPS to retain the Poincaré–Steklov interpretation at both continuous and discrete levels.