This program is simple two dimensional FDTD program in free space with gaussian
source in position at (jc,kc) cell of the problem space with PML absorbing boundary conditions.
When i run this program i get the same result graph in the book for nmax=100
But when i run the program for nmax=800, reflection occurs. Is this due error in the program or what is the reason? Regards
yes, the theory in sullivan's book is not quit right. the PML sullivan deduced is only applicable for 1D and in the edge part of the PML, at the corner where both sigx, sigy exit, the program is not right . you must using the original berenger's PML theory and deduce the correct formula
and CD of taflove's book contains 2D FDTD example with PML (matlab)
I think it is not easy to get the source file written in c.
So, I recommend you to read Berenger's original paper.
It is easy to understand his idea with some figures.
Engmfarhat !.. I think you mean 2D PML program that is working chek it out again at T = 1860 everything is absorbed for 2D but in 3D there is problem i Think in Sullivan's book or may be I am doing the mistake the program in C codes that i am running is already in the form you can download it the link is below
and i fyou have other working code please could you share them in the link that we give see you soon
Dear kubba , wher is the example of the PML in Taflove s book you mean the chapter 7 i think , other wise there is no written code fro pml in taflove's book as much as i know,
The main website of ToyFDTD is removed the authors personal page. I find the now site but lost it very next time i restarted (Thanks to the virus). Fortunately i manged to save a copy of the code. and here it is