Dear Steven,
I wonder if Harminv is necessary to calculate band structure, for I get band
structure without using Harminv. The code is as follows,
************************************************************************************************
(define-param eps 13) ; dielectric constant of waveguide
(define-param w 1.2) ; width of waveguide
(define-param r 0.36) ; radius of holes
; The cell dimensions
(define-param sy 12) ; size of cell in y direction (perpendicular to wvg.)
(define-param dpml 1) ; PML thickness (y direction only!)
(set! geometry-lattice (make lattice (size 1 sy no-size)))
(set! geometry
(list (make block (center 0 0) (size infinity w infinity)
(material (make dielectric (epsilon eps))))
(make cylinder (center 0 0) (radius r) (height infinity) (material
air))))
(set! pml-layers (list (make pml (direction Y) (thickness dpml))))
(set-param! resolution 20)
(define-param fcen 0.25) ; pulse center frequency
(define-param df 1.5) ; pulse freq. width: large df = short impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center 0.1234 0))))
(define-param k-interp 39) ; # k-points to interpolate, otherwise
(run-k-points 300 (interpolate k-interp (list (vector3 0) (vector3
0.5))))*************************************************************************************************
I can get the same result which is shown on
http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Band_diagram%2C_resonant_modes%2C_and_transmission_in_a_holey_waveguide
My question is if Harminv is not necessary to calculate band structure or
(run-k-points T k-points) call Harminv automatically?
Thank you in advance.
Li
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
