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 _______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
