Dear all
I have a problem about caculating the fill factor in the book of "molding
the flow of light"(table 1.page 60). I used the open soures codes(MPB) to
caculate it and i got the same results in the TM case but the different results
in the TE case. My code is attached as follow:
###########################################################
(define-param r 0.2)
(define-param eps 8.9)
(define-param k-interp 8)
(define A (make dielectric (epsilon eps)))
(set! geometry-lattice (make lattice (size 1 1 no-size)))
(set! geometry
(list
(make cylinder
(material A)
(center 0 0) (radius r) (height infinity))))
(define Gamma (vector3 0 0 0))
(define X (vector3 0.5 0 0))
(define M (vector3 0.5 0.5 0))
(set! k-points (interpolate k-interp (list Gamma X)))
(set-param! resolution 64)
(set-param! num-bands 1)
(run-te)
(get-dfield 1)
(compute-field-energy)
(print "energy in dielectric band 1="
(compute-energy-in-objects (make cylinder (center 0)
(radius 0.2) (height infinity)
(material A))) "\n")
(print "energy in dielectric band 1="
(compute-energy-in-dielectric 1.01 8.9) "\n")
##########################################
Could you tell me whether the code is correct?
In addition, compute-energy-in-objects and compute-energy-in-dielectric got
different results.
Thank you!
E_Chui
___________________________________________________________
好玩贺卡等你发，邮箱贺卡全新上线！
http://card.mail.cn.yahoo.com/
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
