Simulating wave propagation in 3D configurations is becoming a very active area of research. Themain advantage of the BEM is that only the domain boundaries are discretized. As a result, thismethod is well suited to dealing with unbounded domains. However, the standard BEM leadsto fully-populated matrices, which results in high computational costs in CPU time and mem-ory require...