Dear MEEP users, Whilst giving the material-function class a whirl using MEEP version 0.20.3, I noticed a possible bug. After running the following script (mf.ctl) and using h5topng on the *-eps*.h5 and *-mu*.h5, epsilon and mu files, the epsilon plot looked at expected, but the mu plot was completely black (see attached pngs). When I set mu without using material-func there is no problem. Any ideas on where I have gone wrong?
mf.ctl
------------------
(define (mat-func p) (make medium
(epsilon (+ 1 (* (vector3-x p) (vector3-x p)) (* (vector3-y p) (vector3-y
p))))
(mu (+ 1 (* (vector3-x p) (vector3-x p)) (* (vector3-y p) (vector3-y p))))
))
(set! geometry-lattice (make lattice (size 10 10 no-size)))
(set! geometry (list
(make cylinder (center 0 0 0) (radius 2.5) (height infinity)
(material (make material-function (material-func mat-func))))
))
(init-fields)
(output-epsilon)
(output-mu)
(quit)
------------------------------
Also, if both epsilon and mu are less than zero (negative index material
- NIM), is MEEP currently able to accomodate this? I know that to model
realistic NIMs one should use a dispersive model for both epsilon and
mu, but at a single frequency, frequency-independent epsilon and mu
should suffice.
Kind regards,
Jonathan Breeze
Research Fellow
Department of Materials
Imperial College London
Exhibition Road
London SW7 2AZ
Phone: + 44 (0)20 7594 6800
Fax : + 44 (0)20 7594 6757
Email: [EMAIL PROTECTED]
Web : www.imperial.ac.uk/materials
<<attachment: mf-eps-000000.00.png>>
<<attachment: mf-mu-000000.00.png>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
