Hello
I am trying to simulate a perfectly matched layers (PML) for an elastic wave propagating in a solid. For this purpose I need to define a specific function which is called "stretching function". It is defined in the below form:
a0 is a constant and "fi" is an arbitrary function (in X and Y) directions. The question is:
How can I define the above function? As you can see it is a complex function (with a term multiplied by "i" ). As far as I know FEniCS does not support complex numbers and functions. Does anybody know if it is possible to implement it in FEniCS or not? If it is, how can I define it?
Thanks in advance for your help.