Time domain PML complex scaling


I am modifying a 2D time-domain PML code (see attached by Joachim) to incorporate the existence of evanescent waves.

In some other problems, there exist evanescent waves. In this case I want to set the pml scaling parameter has a negative imaginary part to accelerate the decay of evanescent waves. Does the pml support this feature?


Attachment: acoustics_cart_pml_2018-04-17-2.py