Nested Domain Decomposition with Polarized Traces for the 2D Helmholtz Equation

From MaRDI portal
Publication:4569314

DOI10.1137/15M104582XzbMATH Open1394.65136arXiv1510.01831OpenAlexW2962789715WikidataQ129641303 ScholiaQ129641303MaRDI QIDQ4569314FDOQ4569314


Authors: Leonardo Zepeda-Núñez, Laurent Demanet Edit this on Wikidata


Publication date: 28 June 2018

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: We present a solver for the 2D high-frequency Helmholtz equation in heterogeneous, constant density, acoustic media, with online parallel complexity that scales empirically as mathcalO(fracNP), where N is the number of volume unknowns, and P is the number of processors, as long as P=mathcalO(N1/5). This sublinear scaling is achieved by domain decomposition, not distributed linear algebra, and improves on the P=mathcalO(N1/8) scaling reported earlier in [L. Zepeda-N'u~nez and L. Demanet, J. Comput. Phys., 308 (2016), pp. 347-388 ]. The solver relies on a two-level nested domain decomposition: a layered partition on the outer level, and a further decomposition of each layer in cells at the inner level. The Helmholtz equation is reduced to a surface integral equation (SIE) posed at the interfaces between layers, efficiently solved via a nested version of the polarized traces preconditioner [L. Zepeda-N'u~nez and L. Demanet, J. Comput. Phys., 308 (2016), pp. 347-388.]. The favorable complexity is achieved via an efficient application of the integral operators involved in the SIE.


Full work available at URL: https://arxiv.org/abs/1510.01831




Recommendations




Cites Work


Cited In (11)

Uses Software





This page was built for publication: Nested Domain Decomposition with Polarized Traces for the 2D Helmholtz Equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4569314)