Is there a way to implement perfectly matched layers?

[bleyerj]

Hi,
this has been discussed here

How to treat complex damping functions for PMLs?

Hi everyone, I’m building knowledge on perfectly-matched-layers and how to define them on FEniCS, but I’m stuck now. Consider the weak form of the Helmholtz equation in the bounded domain x = [-a^*, a^*] and y = [-b^*, b^*], where a<x<a^* and b<y<b^* define the PML. \large{\int_{\Omega} (\frac{\gamma_y}{\gamma_x} \frac{\partial p}{\partial x} \frac{\partial q}{\partial x} + \frac{\gamma_x}{\gamma_y} \frac{\partial p}{\partial y} \frac{\partial q}{\partial y}} - \kappa^2 \gamma_x \gamma_y p q)…