Question about concentration factor

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