Computational study of the macroscopic stability of plasmas is a challenging multi-scale problem. Implicit time integration can be used to relieve stability constraints due to fast Alfvén waves, and adaptive mesh refinement (AMR) can be used to resolve highly localized solution features. The strong nonlinearities and numerical stiffness of magnetohydrodynamics (MHD) models present further chall...

Domain decomposition (DD) methods are important techniques for designing parallel algorithms for solving partial differential equations. Since the decomposition is often performed using automatic mesh partitioning tools, one can in general not make any assumptions on the shape or physical size of the subdomains, especially if local mesh refinement is used. In many of the popular domain decompos...