Re: Effective epsilon used for each polarization

On Fri, 13 Oct 2006, Minghui Han wrote:
> H field by doing the curl operation in the FFT manner. The discrepancy
> occurs at the media interfaces. So what does MPB really do with the epsilon
> tensor epsilon_inverse? Isn't it right to first multiply each Fourier
> coefficient of H with (k+G) and then do the inverse FFT to transform curl H
> to the spatial domain, and then multiply it with tensor epsilon_inverse to
> get the E field?

That's correct.  If you're trying to do this manually, probably you have a 
bug.  e.g. you could be using G vectors with all one sign instead of equal 
amounts of positive and negative, which would cause high-frequency errors. 
Or some other problem.

(Unless I'm forgetting something tricky about MPB's data output.)

Steven