Hello,
For the simple controlfile below, I have to disable eps-averaging, otherwise
NaNs "spread out" from the corners of the metal piece on top of the dielectric
slab and eventually fill up the whole space, destroying any field. It only
happens if the metal and the dielectric piece are <= 0.1 apart (parameter d),
though. If eps-averaging is off, no problems.
I am using meep 0.10.1 / libctl 3.0.2 / Guile 1.8.4
Marcel
---------------------------------------------------------------
(define-param d 0)
(set! geometry-lattice (make lattice (size 50 20 no-size)))
(define semicon
(make block (center 0 -5 0) (size infinity 10 0)
(material (make dielectric (epsilon 12.5) ))) )
(define metalm
(make block (center 0 (+ d 0.5) 0) (size 15 1 0)
(material (make perfect-metal))) )
(set! geometry (list semicon metalm))
(set! pml-layers (list (make pml (thickness 1))))
(set! sources (list (make source
(src (make continuous-src (frequency 0.2)))
(component Hz) (end-time infinity)
(center 0 -5 0) (size 0 10 0)
)))
(run-until 1
(at-beginning output-epsilon)
(at-every 0.1 (output-png+h5 Hz "-Zc bluered -m -1 -M 1 -C $EPS"))
)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
