Finite Element Time Domain implementation of C-PML

I'm in process of doing a Finite-Element Time-Domain implementation of PML using the recursive convolution technique (CPML).

I suspect that I don't carry out the recursive convolution update correctly at each time step or that my global matrices are not assembled correctly. I really need some help on this issue. It would be very helpful if anybody can provide me a 1D or 2D code in whatever programming language.

Regards René

hi, it seems we hv the same problems, i need helps for matlab, may be we can cooperate, come on my email afsharibawili@gmail.com

Dear Rene,

It seems I have the same problem as yours. Have you solved the problem now? If so, could you send me the related part of your code for reference? The file in fortran language would be preferable. It would be very helpful. My email is yujianxing@126.com. Thank you very much.

Regards, Jeremy