I have a problem about a triangular lattice which primitive cell has oblique geometric elements. The .ctl program runs fine, no erros in the terminal. But when I use h5topng to see the epsilon.png, I get wrong png. On tutorial online, we have some light - a mpb-data program is showed. But I still have the wrong png. The right png show a 'X' and a '+' one above the other. The png attached the shows the picture inclined. Below, my .ctl code and the png picture.
define-param dpml 1.0)
(define-param ind 1)
(define-param epdiag (* 3.5 3.5))
(set! num-bands 8);
(set! k-points (list (vector3 0 0 0) ; Gamma
(vector3 0 0.5 0) ; M
(vector3 (/ -3) (/ 3) 0) ;K
(vector3 0 0 0))) ;Gamma
(set! default-material (make dielectric-anisotropic (epsilon-diag epdiag
epdiag epdiag)))
(set! k-points (interpolate 4 k-points))
(set! geometry
(list
(make block
(center 0 0)
(size 0.1 0.9 infinity)
(material (make dielectric (index ind)))
(e1 1 0 0)
(e2 -1 1 0)
(e3 0 0 1)
)
(make block
(center 0 0)
(size 0.1 0.9 infinity)
(material (make dielectric (index ind)))
(e1 1 0 0)
(e2 1 1 0)
(e3 0 0 1)
)
(make block
(center 0 0)
(size 0.1 0.9 infinity)
(material (make dielectric (index ind)))
)
(make block
(center 0 0)
(size 0.9 0.1 infinity)
(material (make dielectric (index ind)))
)
)
)
(set! geometry-lattice (make lattice (size 1 1 no-size)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)
))
(set! resolution 64)
(run-tm
)
(run-te
)
--
Cícero Julião da Silva Junior
____________________________________
Institute of Physics - UFAL
www.ufal.br
<<attachment: epsilon.png>>
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
