On Sep 6, 2009, at 11:53 AM, Jakob Meineke wrote:
Code (I am happy to provide more, but did not want to clutter the mail
unnecessarily):

(define (atoms-func r)
   (make dielectric (epsilon 1)
      (E-polarizations
(make polarizability (omega omegaresonance) (gamma linewidth) (sigma (atomic-density r Rx Ry))))))

(set! extra-materials (list atoms-func))

extra-materials should be a list of materials. In particular, you want something like

(set! extra-materials (list (atoms-func 0)))

This sets the extra-materials list to contain one polarizable material of the same omega and gamma as in your materials-function, which is sufficient to tell Meep to add a polarizability of this form and look for it at each point in space (by evaluating your material function).

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

Reply via email to