output-dpwr-in-objects

Hi All,

Continuing with the H1 cavity, the additional line below should output the energy density of all the bands, since the condition is "greater than 0" but nothing is generated (no error but no output). What is wrong?

Again, Thanks,

Fred

...
(run-zeven)

(output-dpwr-in-objects output-dpwr 0 (make cylinder (center 0)
			(radius (- 1 r)) (height h)
			(material (make dielectric (epsilon eps)))))
...

reminder of the full prog:

; A triangular lattice of air holes with H1 defect


(define-param sx 5)
(define-param sy 5)
(define-param sz 2)
(define-param resxy 16)
(define-param resz 10)
(define-param bands 50)

(define r 0.36)
(define h 0.75)

(define eps 8.4)
(set! mesh-size 6)
(set! num-bands bands)

(define Gamma (vector3 0))
(define M0 (vector3 0 0.5 0))
(define K0 (vector3 (/ -3) (/ 3) 0))

(set! geometry-lattice (make lattice (size sx sy sz)
                         (basis1 (/ (sqrt 3) 2) 0.5)
                         (basis2 (/ (sqrt 3) 2) -0.5)))

(set! geometry   (list (make block (material air)
			(center 0)
			(size infinity infinity infinity))
		     (make block (material (make dielectric (epsilon eps)))
			(center 0)
			(size infinity infinity h))))

(set! geometry (append geometry (geometric-objects-lattice-duplicates (list (make cylinder
                       (center 0) (radius r) (height infinity)
                       (material air))))))

(set! geometry (append geometry
			(list (make cylinder
				(center 0) (radius r) (height h)
				(material (make dielectric (epsilon eps)))))))
				
(set! k-points (list M0))

(set! resolution (vector3 resxy resxy resz)) 

(set! tolerance 1e-5)

(set! filename-prefix "3DH1-")

(run-zeven)

(output-dpwr-in-objects output-dpwr 0 (make cylinder (center 0)
			(radius (- 1 r)) (height h)
			(material (make dielectric (epsilon eps)))))


-- 
Frederic Brossard
Hitachi Cambridge Laboratory
Cavendish Laboratory
J J Thomson Avenue
Cambridge CB3 0HE