The method of polarized traces for the 2D Helmholtz equation

8 years 1 months ago

Download math.mit.edu

We present a solver for the 2D high-frequency Helmholtz equation in heterogeneous acoustic media, with online parallel complexity that scales optimally as O(N L ), where N is the number of volume unknowns, and L is the number of processors, as long as L grows at most like a small fractional power of N. The solver decomposes the domain into layers, and uses transmission conditions in boundary integral form to explicitly deﬁne “polarized traces”, i.e., up- and down-going waves sampled at interfaces. Local direct solvers are used in each layer to precompute traces of local Green’s functions in an embarrassingly parallel way (the oﬄine part), and incomplete Green’s formulas are used to propagate interface data in a sweeping fashion, as a preconditioner inside a GMRES loop (the online part). Adaptive low-rank partitioning of the integral kernels is used to speed up their application to interface data. The method uses second-order ﬁnite diﬀerences. The complexity scalings ar...

Leonardo Zepeda-Núñez, Laurent Deman

Real-time Traffic