Full text loading...
-
Multifrontal Hierarchically Semi-separable Solver for 3D Helmholtz Problem Using 27-point Finite-difference Scheme
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 77th EAGE Conference and Exhibition 2015, Jun 2015, Volume 2015, p.1 - 5
Abstract
We present a multi-frontal hierarchically semi-separable solver to perform forward modeling of the 3D Helmholtz acoustic problem. Our frequency-domain solver combines two efficient approaches. First, it uses an optimal 27-point finite-difference scheme to decrease numerical dispersion and reduce required discretization of the model in terms of points per wavelength from 15 to about 4. Second, it uses a supernodal multi-frontal method based on low-rank approximation and hierarchically semi-separable (HSS) structure to improve performance, decrease memory usage and make it practical for realistic size 3D models required by full-waveform inversion. We validate our new solver using a 3D synthetic model and benchmark it against other codes.