Print Email Facebook Twitter Reduction of computing time for least-squares migration based on the Helmholtz equation by graphics processing units Title Reduction of computing time for least-squares migration based on the Helmholtz equation by graphics processing units Author Knibbe, H.P. (TU Delft Numerical Analysis) Vuik, Cornelis (TU Delft Numerical Analysis) Oosterlee, C.W. (TU Delft Numerical Analysis; Centrum Wiskunde & Informatica (CWI)) Date 2015-12-30 Abstract In geophysical applications, the interest in least-squares migration (LSM) as an imaging algorithm is increasing due to the demand for more accurate solutions and the development of high-performance computing. The computational engine of LSM in this work is the numerical solution of the 3D Helmholtz equation in the frequency domain. The Helmholtz solver is Bi-CGSTAB preconditioned with the shifted Laplace matrix-dependent multigrid method. In this paper, an efficient LSM algorithm is presented using several enhancements. First of all, a frequency decimation approach is introduced that makes use of redundant information present in the data. It leads to a speedup of LSM, whereas the impact on accuracy is kept minimal. Secondly, a new matrix storage format Very Compressed Row Storage (VCRS) is presented. It not only reduces the size of the stored matrix by a certain factor but also increases the efficiency of the matrix-vector computations. The effects of lossless and lossy compression with a proper choice of the compression parameters are positive. Thirdly, we accelerate the LSM engine by graphics cards (GPUs). A GPU is used as an accelerator, where the data is partially transferred to a GPU to execute a set of operations or as a replacement, where the complete data is stored in the GPU memory. We demonstrate that using the GPU as a replacement leads to higher speedups and allows us to solve larger problem sizes. Summarizing the effects of each improvement, the resulting speedup can be at least an order of magnitude compared to the original LSM method. Subject s Least-squares migrationHelmholtz equationWave equationFrequency domainMultigrid methodGPU accelerationMatrix storage formatFrequency decimation To reference this document use: http://resolver.tudelft.nl/uuid:0a2bcdfd-6b4d-433e-8aae-d876234b9703 DOI https://doi.org/10.1007/s10596-015-9546-z ISSN 1420-0597 Source Computational Geosciences: modeling, simulation and data analysis, 20 (2), 297-315 Part of collection Institutional Repository Document type journal article Rights © 2015 H.P. Knibbe, Cornelis Vuik, C.W. Oosterlee Files PDF 7008244.pdf 1.96 MB Close viewer /islandora/object/uuid:0a2bcdfd-6b4d-433e-8aae-d876234b9703/datastream/OBJ/view