Hello
Meep is doing something I don't understand. I have a geometry that I
have computed fields on before, so I don't think there's anything
wrong with it.
I define the following functions:
;txx xx comp of stress tensor
(define (txx ex ey ez hx hy hz)
(- (+ (* ex ex) (* hx hx))
(* 0.5 (+ (* ex ex) (* ey ey) (* ez ez) (* hx hx) (* hy hy) (*
hz hz))) )
)
;tyx yx comp of stress tensor
(define (tyx ex ey ez hx hy hz)
(+ (* ex ey) (* hx hy))
)
; txx+tyx
(define (txxplustyx ex ey ez hx hy hz)
(+
(- (+ (* ex ex) (* hx hx))
(* 0.5 (+ (* ex ex) (* ey ey) (* ez ez) (* hx hx) (* hy hy) (*
hz hz))) );txy
(+ (* ex ey) (* hx hy)) ;tyx
)
)
(define (txxplustyx-output)
(output-field-function "txxplustyx-output" (list Ex Ey Ez Hx Hy Hz)
txxplustyx))
As far as I can tell this all fits the instructions given here
http://ab-initio.mit.edu/wiki/index.php/Meep_field-function_examples.
Then I try to run a simulation
(run-until (/ 100 fcen) (at-beginning output-epsilon)
(synchronized-magnetic (at-every (/ 1 fcen 2) txxplustyx-output)))
Meep outputs the epsilon and then tells me there is a "Segmentation Fault".
I would be grateful if somebody could explain this problem and how to fix it.
Alan Forrester
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
