Short implementation of acoustic wave propagation using finite-differences in time domain and CUDA. The code is solving second order wave equation in pressure formulation, O(2,8). Computational domain is surrounded by reflective boundaries.
The code is intentionally written in a single file to simplify understanding and startup.
CUDA implementation details:
- 2D shared memory tiles
- all kernels are executed in a single stream as it is convenient for a toy application.
Example:
nvcc -O3 -use_fast_math cuwaveprop2d.cu -o cuwaveprop2d && ./cuwaveprop2d
The output snapshots are stored as binary files in snap/ folder. You might visualize these by running show_snaps.m.
Credit:
We use CHECK() function to decode CUDA error codes and overall guidance from this book
@book{cheng2014professional,
title={Professional Cuda C Programming},
author={Cheng, John and Grossman, Max and McKercher, Ty},
year={2014},
publisher={John Wiley \& Sons}
}
We also thank @armandoespindola from KAUST for sharing his expertise in CUDA programming