I still have questions
1) I want to evaluated "ezsum" at every point. So I can make a picture.
However, (print "Ezsum = " Ezsum "\n") only output a number.
Can any method help me??
2) In meep-discuss, I find a function (meep-fields-dt-get fields). I can not
understand it. What can I consult to understand such kinds function??
>> I'm doing a 3D simulation. I want to calculate the integral of
>> |Ez|^2/(|E|^2) over 200 time steps at the end of the sumulation time. I
>This doesn't seem well-defined. What if |E| = 0?
>> need to calculate the whole field. I have looked through the
>>meep-discuss but have no idea. I write the function like:
>>(define Ezsum 0)
>>(define (do-Ez-sum)
>>(set! Ezsum (+ Ezsum
>>(integrate-field-function (list Ez Ex Ey)
>>(lambda (ez ex ey) (/ (sqrt ez) (+ (sqrt ex) (sqrt ey) (sqrt ez))))))))
>>(after-time 300 (run-until 500 do-Ez-sum))
>>Is there anything wrong ???
>I can spot at least two errors.
>First, you should have (lambda r ez ex ey) ...), since
>integrate-field-function alwasy passes the position as the first argument.
>See:
>http://ab-initio.mit.edu/wiki/index.php/Meep_field-function_examples
>Second, after-time should go inside the run statement, i.e.
>(run-until 500 (after-time 300 do-Ez-sum))
>>How can I output the results?? It seems that "output-field-function" is
>>not applicable. Can anyone give me a sample code???
>(print "Ezsum = " Ezsum "\n")
---------------------------------
雅虎免费邮箱，全球第一邮箱品牌！_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
