I developed a 2D/3D acoustic wavefield simulator as an addition to the MADAGASCAR seismic package. It is GPU enabled and can operate in both cartesian and spherical coordinate systems.
This program was originally developed to act as an analogue to simulate radar propagation through the ice sheets of Mars. Due to the large processing requirements for simulating over such large regions of the martian polar deposits, an acoustic wavefield model was chosen as a proxy as it does not use nearly as much computational power as a full waveform finite difference based EM simulator would require.
This also comes with the benefits of the package being largely adaptable and having applications for seismic/acoustic applications.
The source code and compiled versions can be found at on GitHub here. Simple examples of constant velocity exist in each model directory.
In order to run the GPU code a CUDA enabled GPU is required as well as sufficient memory for the task at hand. Additionally, the Madagascar software needs to be installed. If the ./bashrc file on the system is configured properly, the python interface for the packages should automatically detect if the system is GPU enabled or not. This can be done by defining "nGPU" as greater or equal to 1 in the ./bashrc.