A method is presented to simulate the propagation of elastic waves in a 2-D piecewise homogeneous domain composed of several regions of arbitrary shapes. The regions can be isolated within an elastic background or also embedded among themselves. The method is based on the indirect boundary element method (IBEM) in which the diffracted field is obtained as the contributions of distributed source…