A two-scale method using a list of active sub-domains for a fully parallelized solution of wave equations

Abstract

tWave form modeling is used in a vast number of applications. Therefore, different methods have beendeveloped that exhibit different strengths and weaknesses in accuracy, stability and computational cost.The latter remains a problem for most applications. Parallel programming has had a large impact onwave field modeling since the solution of the wave equation can be divided into independent steps. Thefinite difference solution of the wave equation is particularly suitable for GPU acceleration; however,one problem is the rather limited global memory current GPUs are equipped with. For this reason, mostlarge-scale applications require multiple GPUs to be employed. This paper proposes a method to optimallydistribute the workload on different GPUs by avoiding devices that are running idle. This is done by usinga list of active sub-domains so that a certain sub-domain is activated only if the amplitude inside the sub-domain exceeds a given threshold. During the computation, every GPU checks if the sub-domain needsto be active. If not, the GPU can be assigned to another sub-domain. The method was applied to syntheticexamples to test the accuracy and the efficiency of the method. The results show that the method offersa more efficient utilization of multi-GPU computer architectures.