Floating point exception

Lukas Oliva <[email protected]> Thu, 14 Sep 2006 00:05:54 +0200

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

 I attach te source file.

  Lukas Oliva

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

kr.ctl (text/plain, 2.2 KB)

(define-param r1 0.2)
(define-param r2 0.1)
(define-param eps1 5)
(define-param eps2 1)

(set! geometry-lattice (make lattice (size 1 1 no-size)
                         (basis1 0.5 0)
                         (basis2 0 0.5)))

(set! geometry
      (list (make cylinder
	      (center 0 0 0) (radius r1) (height infinity)
	      (material (make dielectric (epsilon eps1))))
	    (make cylinder
	      (center 0 0 0) (radius r2) (height infinity)
	      (material (make dielectric (epsilon eps2))))))

; The k-points list, for the Brillouin zone of a triangular lattice:
(set! k-points (list (vector3 0 0 0)          ; Gamma
                     (vector3 0 0.5 0)        ; M
                     (vector3 0.5 0.5 0) ; K
                     (vector3 0 0 0)))        ; Gamma
(define-param k-interp 4) ; number of k-points to interpolate
(set! k-points (interpolate k-interp k-points))

(set-param! resolution 32)
(set-param! num-bands 8)

(define (tm_gap1 r1 r2 eps1 eps2)
    (set! geometry
      (list (make cylinder
	      (center 0 0 0) (radius r1) (height infinity)
	      (material (make dielectric (epsilon eps1))))
	    (make cylinder
	      (center 0 0 0) (radius r2) (height infinity)
	      (material (make dielectric (epsilon eps2))))))

    (print "Conf: r1 = "r1" r2 = "r2" eps1 = " eps1"eps2 = "eps2"\n")
  (if ( (r1 <= r2)||( eps1<=eps2))
    (print "Chyba: r1 = "r1" r2 = "r2" eps1 = " eps1"eps2 = "eps2"\n"))

    (run-tm)
    (retrieve-gap 1))

(define (te_gap1 r1 r2 eps1 eps2)
    (set! geometry
      (list (make cylinder
	      (center 0 0 0) (radius r1) (height infinity)
	      (material (make dielectric (epsilon eps1))))
	    (make cylinder
	      (center 0 0 0) (radius r2) (height infinity)
	      (material (make dielectric (epsilon eps2))))))

    (print "Conf: r1 = "r1" r2 = "r2" eps1 = " eps1"eps2 = "eps2"\n")
  (if ( (r1 <= r2)||( eps1<=eps2))
    (print "Chyba: r1 = "r1" r2 = "r2" eps1 = " eps1"eps2 = "eps2"\n"))

    (run-te)
    (retrieve-gap 1))

 (run-tm)
 (define-param first_tm (retrieve-gap 1))
 (run-te)
 (define-param first_te (retrieve-gap 1))
(print "First TM maximum: " first_tm"\nFirst TE maximum :" first_te"\n")

(maximize-multiple tm_gap1 0.001)
(maximize-multiple te_gap1 0.001)