The method of polarized traces for the 3D Helmholtz equation
Author(s)Zepeda Nunez, Leonardo Andres; Scheuer, Adrien; Hewett, Russell; Demanet, Laurent
MetadataShow full item record
We present a fast solver for the 3D high-frequency Helmholtz equation in heterogeneous, constant density, acoustic media. The solver is based on the method of polarized traces, coupled with distributed linear algebra libraries and pipelining to obtain an empirical online runtime O(max(1, R/n)N log N) where N = n³ is the total number of degrees of freedom and R is the number of right-hand sides. Such a favorable scaling is a prerequisite for large-scale implementations of full waveform inversion (FWI) in frequency domain.
DepartmentMassachusetts Institute of Technology. Department of Mathematics; Massachusetts Institute of Technology. Earth Resources Laboratory
Society of Exploration Geophysicists
Zepeda-Núñez, Leonardo, Adrien Scheuer, Russell J. Hewett, and Laurent Demanet, "The Method of Polarized Traces for the 3D Helmholtz Equation." Geophysics 84, 4 (July 2019): T313-T333. ©2019 Author(s)
Author's final manuscript