------

        Steven-

believe i found a minor error, but it's probably just numerics. in a control file, place:

(define a 1)
(define sx (/ a (sqrt 2)))
(define sy (* sx (sqrt 3)))

(set! geometry-lattice (make lattice (size sx sy 10)))

(set-param! k-point (vector3 0))

(set! pml-layers (list (make pml (thickness 1) (direction Z) (side Low)) (make pml (thickness 1) (direction Z) (side High))))

(set-param! resolution 32)

(set! sources (list (make source (src (make continuous-src (frequency 1))) (component Ex) (size sx sy 0) (center 0 0 -4.5))))

(run-until 20 (at-end output-efield-x))


running meep (either in linux or a mac), i get:

-----------
Initializing structure...
Working in 3D dimensions.
time for set_epsilon = 0.484472 s
-----------
meep: Source width > cell width in y direction!


replacing sy in the size for source, w/ (- sy 6e-3) avoids the problem. of course, 6e-3 depends on the resolution used...

thanks,

        gp

------


_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to