On Wed, 27 Jun 2007, matt wrote:
> I'm solving a simple multilayer loss-free dielectric slab scattering problem,
> for the case of normal incidence.
>
> The problem is simple enough that it can be solved analytically. The
> analytic solution matches published results, so I have confidence in its
> correctness.
>
> As you can see in the attached plot, the resonance in the reflectivity is
> frequency shifted by about 10 percent, and the scaling isn't quite right.
Hi Matt,
The Meep solution should definitely converge to the exact result as you
increase resolution etc. (In fact, transmission through a multilayer film
is one of the test cases: see tests/bragg_transmission.cpp.)
Someone posted with a similar problem in March - they were comparing the
transmission through a thin film to the analytical formula, and weren't
getting convergence:
http://thread.gmane.org/gmane.comp.science.electromagnetism.meep.general/746
The problem turned out to be that they were calling (set! resolution ...)
too late in the file. When you call add-flux to add a flux plane the
fields are initialized and the resolution etc. is set; changing the
resolution etc. after that point has no effect (unless you call
reset-meep! or similar). Perhaps your problem is similar?
Regards,
Steven G. Johnson
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
