1887

Abstract

Summary

Half precision floating point numbers is becoming increasingly supported by new processors, often with a significant throughput gain over single precision operations. In this article, we investigate whether half precision is suitable for finite difference based seismic modeling in the context of imaging and inversion. By scaling the finite difference expression of the isotropic elastic wave equation, we manage to obtain a stable solution despite the very narrow dynamic range of the half-precision format. We present a CUDA implementation of this code, which, on most recent GPUs, is nearly twice as fast and uses half the memory of the equivalent single precision version. The error on seismograms caused by the reduced precision is shown to correspond to a fraction of a percent of the total seismic energy, and is mostly incoherent with seismic phases. Thus, half precision modeling could accelerate full waveform inversion or migration with a negligible impact on the quality of their output.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201801351
2018-06-11
2024-04-18
Loading full text...

Full text loading...

References

  1. Billette, F. and Brandsberg-Dahl, S.
    [2005] The 2004 BP velocity benchmark. In: 67th EAGE Conference & Exhibition.
    [Google Scholar]
  2. Fabien-Ouellet, G., Gloaguen, E. and Giroux, B.
    [2017] Time-domain seismic modeling in viscoelastic media for full waveform inversion on heterogeneous computing platforms with OpenCL. Computers & Geosciences, 100, 142–155.
    [Google Scholar]
  3. Foley, D. and Danskin, J.
    [2017] Ultra-Performance Pascal GPU and NVLink Interconnect. IEEE Micro, 37(2), 7–17.
    [Google Scholar]
  4. Graves, R.W.
    [1996] Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences. Bulletin of the Seismological Society of America, 86(4), 1091–1106.
    [Google Scholar]
  5. Komatitsch, D., Erlebacher, G., Göddeke, D. and Michéa, D.
    [2010] High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster. Journal of Computational Physics, 229(20), 7692–7714.
    [Google Scholar]
  6. Micikevicius, P.
    [2009] 3D finite difference computation on GPUs using CUDA. In: Proceedings of 2nd workshop on general purpose processing on graphics processing units. 79–84.
    [Google Scholar]
  7. Virieux, J.
    [1986] P-SVwave propagation in heterogeneous media: Velocity-stress finite-difference method. Geophysics, 51(4), 889–901.
    [Google Scholar]
  8. Zuras, D., Cowlishaw, M., Aiken, A., Applegate, M., Bailey, D., Bass, S., Bhandarkar, D., Bhat, M., Bindel, D. and Boldo, S.
    [2008] IEEE standard for floating-point arithmetic. IEEE Std 754-2008, 1–70.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201801351
Loading
/content/papers/10.3997/2214-4609.201801351
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error