Half precision wave simulation

TL;DR

In this work, the negative impact caused by the accumulation of roundoff errors in wave simulations is illustrated with carefully designed numerical experiments and the corresponding remedy in the form of compensated sum is provided.

Abstract

In recent years, half precision floating-point arithmetic has gained wide support in hardware and software stack thanks to the advance of artificial intelligence and machine learning applications. Operating at half precision can significantly reduce the memory footprint comparing to operating at single or double precision. For memory bound applications such as time domain wave simulations, this is an attractive feature. However, the narrower width of the half precision data format can lead to degradation of the solution quality due to larger roundoff errors. In this work, we illustrate with carefully designed numerical experiments the negative impact caused by the accumulation of roundoff errors in wave simulations. Specifically, the energy-conserving property of the wave equations is employed as a convenient diagnosis tool. The corresponding remedy in the form of compensated sum is then provided, with its efficacy demonstrated using numerical examples with both acoustic and elastic wave equations on hardware that support half precision arithmetic natively.

Figures (27)

• Figure 1: Illustration of the grid layout in space.
• Figure 2: Illustration of the grid layout in time.
• Figure 3: Time history of $P$ at the receiver location.
• Figure 4: Time history of $V_x$ at the receiver location.
• Figure 5: Time history of $V_y$ at the receiver location.