Hi Christopher,
I have played with pml for Helmholtz, and it works great.
But I got some trouble turning to Maxwell.
I tried to use the same setup for the scatering.py problem in pde_tutorial/pml/, by only changing finite element to Hcurl and bilinear form to curl*curl
It gave me the error:
PML not supported for diffop curl
Also, I noticed the old way:
a += BFI('PML_laplace',coef=1)
stop working for pml, and the error is:
RuntimeError: please set global variable constant_table_for_FEM
Any idea how to fix this?
Thanks in advance,
Guosheng