; Copyright (C) 1998, 1999, 2000, 2001, 2002, Steven G. Johnson
;
; This file may be used without restriction.  It is in the public
; domain, and is NOT restricted by the terms of any GNU license.
;
; This library is distributed in the hope that it will be useful,
; but WITHOUT ANY WARRANTY; without even the implied warranty of
; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
; Lesser General Public License for more details. 

; Sample control file (for example.scm specification).
; You should edit this file in Scheme mode if you are using emacs:
; type M-x scheme-mode or modify your .emacs file to automatically
; put .ctl files in Scheme mode.

; Define some dielectric materials:
(define GaAs (make material-type (epsilon 11.56)))
(define AlOx (make material-type (epsilon 2.25)))

; Set the dimensions.  Use set! rather than define to change the value
; of an existing variable.
(set! dimensions 2)

; Set the k-point list:
(set! k-points 
  (list
   (vector3 0 0)
   (vector3 0.5 0)
   (vector3 0.5 0.5)
   (vector3 0 0)))

; Reset the k-point list for fun, in a fancier way:
; Here, we will use the built-in interpolate function to interpolate
; points between the corners of the Brillouin zone.

(define Gamma-point (vector3 0 0))
(define X-point (vector3 0.5 0))
(define M-point (vector3 0.5 0.5))
(set! k-points (interpolate 4 (list Gamma-point X-point M-point Gamma-point)))

; Define a parameter, R, that can be set from the command-line (e.g.
; by "R=0.3", and which assumes a default value of 0.1 if not set:
(define-param R 0.1)  ; (used as the radius of the sphere, below)

; Set the geometry:

(set! geometry
      (list
       (make cylinder (material air) (center 1.5 0)
	(radius 0.2) (height 0.5) (axis 3 4))
       (make sphere (material GaAs) (center 1.1 -1.1) (radius R))
       (make block (material AlOx) (center -1 0.5) (size 0.2 1.3))
       (make block (material air) (center -1.5 -1) (size 1.2 0.3) (e2 1 1))))

; Append a 2d 3x3 lattice of cylinders to the geometry:

(set! geometry-lattice (make lattice (size 5 5 1)))

(set! geometry
      (append geometry
	      (geometric-objects-lattice-duplicates
	       (list (make cylinder (material GaAs) 
			   (center 0 0) (radius 0.1) (height 1.0))))))

(run)

; The (run) command normally turns interactive mode off; turn it back
; on so that we are dropped into a Guile command line:
(set! interactive? true)