Operator Expansions for Linear Waves: Parallel Implementation and Multilayer Inversion

In this contribution we discuss a novel Boundary Perturbation approach to compute general Dirichlet–Neumann Operators arising in elastodynamics in a rapid, high–order, and robust fashion. In addition, we have also implemented a parallel version of the algorithm using OpenMP and FFTW-MPI, and investigated the inversion method for 2D interface reconstruction using the Operator Expansions approach.