On Wed, 23 Aug 2006, Matt wrote:
2)
As soon as I tell meep a structure is periodic with this
(set! geometry-lattice (make lattice (size a a no-size))) ; 2d periodic
the h5 epsilon data becomes 2d and not 3d. This isn't really what you'd want
if you're trying to visualize a unit cell, right? I'm not sure, but i think
this has something to do with (make lattice (size a a no-size))
having an undefined third dimension.
The third dimension isn't undefined. You've explicitly told it that it
has "no size". i.e. you've specified a 2d computational cell, so it's not
surprising that this is what you get.
3) If the structure is periodic, and you define in the geometry a (size
inf inf finite) slab , does (output-epsilon) truncate the h5 output to
the size defined by the lattice?
You definitely don't want to define the cell size to be infinite. Unless
you have a computer lying around with infinite memory and speed. (If you
do, can I have an account?)
4) I can see where metal is (or where dielectric isn't) with h5topng,
but not with mayavi. Are there any rules for how (output-epsilon)
behaves when dielectric overlaps with dielectric or metal?
output-epsilon should just spit out a large negative epsilon value for
metal. How you visualize this is between you and mayavi.
Cordially,
Steven G. Johnson
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss