Dear Steven and meep users, I have a question regarding the MEEP field-function examples as given in
http://ab-initio.mit.edu/wiki/index.php/Meep_field-function_examples On my personal computer I have meep 0.10 and on a cluster computer I have access to meep-mpi 0.20-3. When I run the example functions from the above website, they work just fine on meep 0.10, namely > (define (f r ex hz eps) > (- (+ (* (vector3-x r) (vector3-norm r)) ex) (* eps hz))) > (print "The integral of our weird function is: " > (integrate-field-function (list Ex Hz Dielectric) f) "\n") gives a number as a result. However, when I run the same code on the cluster computer, I get "##,###,##i" instead of a number. If I remove the dependence of the function on "Dielectric", then everything works fine on the cluster computer as well: > (define (f r ex hz) > (- (+ (* (vector3-x r) (vector3-norm r)) ex) h))) > (print "The integral of our weird function is: " > (integrate-field-function (list Ex Hz) f) "\n") Could it be that meep 0.20-3 has problems with "Dielectric"? (I'm still in the process of trying to install 0.20 on my own computer to do do the testing) Or is it the fact that on the cluster I am using meep-mpi? I checked the Meep release notes (http://ab-initio.mit.edu/wiki/index.php/Meep_release_notes) and found out that in meep 0.20, the following change was made: <Renamed "|dielectric|" to "|medium|" in libctl interface,> That is why in the above code, I tried using "Medium" instead of "Dielectric", however in that case Meep exits with an error saying that "Medium" is an unbound variable. Any suggestions will be greatly appreciated. As a side question, on the same field-function example website above the following example is given: > (print "The integral of our weird function from (-1,0,0) to (1,0,0) is: " > (integrate-field-function (list Ex Hz Dielectric) f (volume (size 1 0 > 0) (center 0 0 0))) "\n") > My question is if (volume (size 1 0 0) (center 0 0 0)), doesn't that give the integral from (-0.5, 0, 0) to (0.5, 0, 0) as opposed to (-1, 0, 0) to (1, 0, 0)? I was just wondering if that is a typo, or a misunderstanding on my part. Thanks a lot, George _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
