MPB vs MPBI discrepancy

MHerrera <[email protected]>
Newsgroups gmane.comp.science.photonic-bands
Message-ID <[email protected]>
Hello,

When I run the following ctl file using MPBI vs MPB, I get different results
between the two, though I believe I should be getting the same output:

*******************************************************************
(set! default-material (make dielectric (index 1.0)) )  
(define-param eps_air 1.0)
(define-param eps_sub 12)
(define-param thickness_air 10.0)
(define-param thickness_sub 1.0)
(set! geometry-lattice ( make lattice
        (basis1 1 0 0)
        (basis2 0 1 0)
        (basis3 0 0 1)
        (size 1 1 10)))
(set! num-bands 8)
(set! grid-size (vector3 16 16 160))
(set! k-points (list (vector3 0 0 0)     ; Gamma
                     (vector3 0.5 0 0)   ; X
                     (vector3 0.5 0.5 0) ; M
                     (vector3 0 0 0)))   ; Gamma

(set! k-points (interpolate 4 k-points))
(set! geometry
        (list
; PHOTONIC LAYER
         (make cylinder
                (material (make dielectric (epsilon eps_sub)))
                (center 0 0 0)
                (radius 2.0)
                (height thickness_sub))

; AIR HOLES IN SUBSTRATE
         (make cylinder
                (material (make dielectric (epsilon eps_air)))
                (center 0 0 0)
                (radius 0.4)
                (height thickness_air) 
          )
         )
)
(run-even)

*****************************************************************
Any help that could be provided would be greatly appreciated.
Mark Herrera
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.