Hi,
considering bug.ctl:
(set! eps-averaging? false)
(set-param! resolution 4)
(define-param sx 8)
(define-param sy 8)
(define-param pmlth 4)
(define sz (+ pmlth 15 pmlth))
(set! geometry-lattice (make lattice (size sx sy sz)))
(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(set! pml-layers (list (make pml (direction Z) (thickness pmlth))))
(set! symmetries (list
(make mirror-sym
(direction X)
(phase -1))
(make mirror-sym
(direction Y)
(phase 1))
))
(define fcen 0.05)
(define df 0.06)
(set! sources (list
(make source
(src (make gaussian-src
(frequency fcen)
(fwidth df)))
(component Ex)
(center 0 0 0)
(size sx sy 0))))
(run-until 310
(to-appended "ex" (at-every 10 output-efield-x)))
;-------------------------------------
If you run it like this, the field blow up infinitely. I can't make it in 2
dimensional simulation. If you comment the symmetry for the Y direction or X
direction, everything is fine.
The blowing up appear at the surface of the pml layer on the side of the
simulation.
The field created by the source looks good in the y=0 and x=0 plane, in any
case (X+Y symmetry or only one).
h5topng -R -0 -x0 -t 0:30 -Zc dkbluered bug-ex.h5
only bug-ex.t30.png is not only white
h5topng -0 -x0 -t 0:30 -Zc dkbluered bug-ex.h5
the pulse from the source looks good
The same for -y0
Well, that's all I can say so far.
Cheers.
--
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
