> I do not know, why p explodes, though. In principle, even if a p is set
> somewhere outside the metal region because of some yee-grid shift, it
> should decay as fast as anywhere else, or not?
>
> The whole problem, btw, disappears if you just put
> if(s[i] != 0){...
> inside the
Dear Meep community,
I sent a problem a couple of days ago and tracked it down a bit. the
problem was, that if you define a non-diagonal sigma in your
lorentzian_susceptibility, the calculation becomes an unstable and
fields and polarizations grow exponentially.
The problem, I think, is a bug in
2 matches
Mail list logo