Domain decomposition: parallel multilevel methods for elliptic partial differential equations
Domain decomposition: parallel multilevel methods for elliptic partial differential equations
Robust Ordering of Sparse Matrices using Multisection
SIAM Journal on Matrix Analysis and Applications
Computer Solution of Large Sparse Positive Definite
Computer Solution of Large Sparse Positive Definite
A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation
Journal of Computational Physics
Flexible Variants of Block Restarted GMRES Methods with Application to Geophysics
SIAM Journal on Scientific Computing
Hi-index | 0.00 |
This is the second paper in a two-part series that describes a massively parallel code that performs 2D full-waveform inversion of wide-aperture seismic data for imaging complex structures. We present several numerical validation of the full-waveform inversion code with both canonical and realistic synthetic examples. We illustrate how different multiscale strategies can be applied by either successive mono-frequency inversions or simultaneous multifrequency inversions and their impact on the convergence and the robustness of the inversion. We present a scalability analysis using a real marine data set recorded by a dense array of ocean bottom seismometers to image the crustal structure of a subduction zone. We obtained a speedup of 20 when using 50 processes on a PC cluster which allowed us to iteratively invert 13 frequencies of the full data set in less than 2 days. This computational performance will allow in the future more extensive analysis of full-waveform tomography methods when applied to representative case studies or when considering 3D geometries.