Clarification and Guidance on Coupling PML Boundary in Devito Solver

[zhangxiaoshuotttt]

Hi team,

I'm trying to incorporate a PML boundary into the solver to make it easier for new users like myself to apply Devito to field data. My primary reference has been your notebook.

However, I'm encountering some uncertainties:

PML Coupling Approach:

I'm not sure if the method I've used to couple the PML boundary is correct. My implementation produces results nearly identical to the original damping method, which suggests I might have made a mistake. I find it challenging to figure out exactly where the error is occurring and would appreciate your insights.

Interaction with Free-Surface Condition:
Because I defined a subdomain specifically for the PML, I'm also unsure if this setup could inadvertently affect the free-surface boundary conditions.

I've provided additional details and code snippets for your reference.

Thanks in advance for your support and patience.

Best regards,
Shuo