The method of polarized traces for the 3D Helmholtz equation

Author(s)
Zepeda Nunez, Leonardo Andres; Scheuer, Adrien; Hewett, Russell; Demanet, Laurent
Thumbnail
DownloadAccepted version (4.022Mb)
Terms of use
Creative Commons Attribution-Noncommercial-Share Alike http://creativecommons.org/licenses/by-nc-sa/4.0/
Metadata
Show full item record
Abstract
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.
Date issued
2019-07
URI
https://hdl.handle.net/1721.1/124220
Department
Massachusetts Institute of Technology. Department of Mathematics; Massachusetts Institute of Technology. Earth Resources Laboratory
Journal
Geophysics
Publisher
Society of Exploration Geophysicists
Citation
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)
Version: Author's final manuscript
ISSN
1942-2156
0016-8033
Collections