Hello,
When I run the scheme code below, I'm receiving the following error:

ERROR: In procedure meep-fields-get-field:
ERROR: No matching method for generic function
`meep_fields_get_field'

I can't seem to debug it.  Does anyone know a quick fix?

Chad



; Calculate band Projected Band Structure (PBS) for 2D PC
; ---------------------------------------------------

(define-param eps 7.45)  ; dielectric constant of n_AlGaAs = 3.3
(effective index = 2.77)
(define-param a 1)  ; lattice parameter 410 nm
(define-param r 0.275)   ; radius of air holes
(define-param dpml 1)  ; PML thickness
(define-param sx 10)     ; size of cell in x direction
(define-param sy 1)      ; size of cell in y direction  ;linear case
run with sy=28
(define-param c (* (/ (sqrt 3) 2) a))    ; define c cell size

(set-param! resolution 12)

; source parameters
(define-param fcen 0.25)
(define-param df 1.5)      ; pulse width (in frequency) - .01 runs
for 1000 units

(set! geometry-lattice (make lattice (size sx sy no-size)))

(set! pml-layers (list (make pml (direction X) (thickness dpml))))

(set! geometry
        (append
                (list (make block (center 0 0 0) (size sx sy no-size) (material
(make dielectric (epsilon eps)))))

                (geometric-objects-lattice-duplicates
                (list
                   (make cylinder (center 0 0 0) (radius r) (height
infinity) (material air))
                   (make cylinder (center c (* -1 .5 a) 0) (radius r) (height
infinity) (material air))
                ) (* a (sqrt 3)) a )

                ;MAKE PCWG
                (list (make block (center 0 0 0) (size a sy infinity) (material
(make dielectric (epsilon eps)))))
        ))

(set! sources (list
                (make source    ;transverse to propagation direction y
                        (src (make gaussian-src (frequency fcen) (fwidth df)))
                        (component Ex) (center 0 0.1234))

        ))

        (set! symmetries (list (make mirror-sym (direction X) (phase +1))))

        (set-param! k-point (vector3 0 0.4))

        (run-sources+ 300 (after-sources (harminv Ey (vector 0 0.1234) fcen
df)))

        (define-param k-interp 10)

        (run-k-points 300 (interpolate k-interp (list (vector3 0 0)
(vector3 0 0.5))))

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

Reply via email to