Hi, all
I have a question related to use of mpi meep.
Our system administrator installed Meep 1.0 in our server with single
processor version and MPI. When I run the following code, the error
showed "unbound variable synchronized-magnetic" in both version. If I remove
synchronized-magnetic, the error showed "too many field components
". The system is debian in the server. No problem to run these codes in my
local computer with Meep 0.20.4 and fedora system. I like to know what
happened with them. Your suggestion to solve it is thanked in advance.
My code is here:
(define-param fiber 450) ;set the diameter of the fiber
(define-param l 20) ;set the length of the fiber
(define-param guide 250) ;set the normalized value
(define-param dfiber (/ fiber guide)) ;normalize the diameter of fiber
(define-param MFD (* (/ dfiber 2) 0.9267)) ;set the MFD value
(define-param wlength (/ 1000 guide)) ;set the normalized wavelength of light
(define-param eps (* 2.0355 2.0355)) ;define the index of fiber
(define-param sx (* 1 l))
(define-param sy 4)
(define-param sz 4)
(define-param sc 4)
(set! geometry-lattice (make lattice (size sx sy sz)))
(define-param hole 10) ;set the diameter of the fiber
(define-param n_hole 1.0)
(define-param rhole (/ hole guide)) ;normalize the diameter of central hole
(set-param! force-complex-fields? true)
(set-param! eps-averaging? false)
(define-param no-inclined? false)
(set! geometry
(list
(make cylinder
(center (+ (* 0.5 l) (* -0.5 sx)) 0 0)
(radius (/ dfiber 2))
(height l )
(axis 1 0 0)
(material (make dielectric (epsilon eps))))
)
)
(define-param gauss? false)
(define-param fcen (/ 1 wlength))
(define-param df 0.1)
(define ((circular-amp w) r)
(if (> (vector3-norm r) (/ dfiber 2))
(* (exp (* -1 (/ (vector3-dot r r) (* w w)))) eps)
(exp (* -1 (/ (vector3-dot r r) (* w w))))))
(set! sources
(if (not gauss?)
(list
(make source
(src (make continuous-src (frequency fcen)(end-time 150) ))
(component Ez)
(center (+ 1 (* -0.5 sx)) 0 0)
(size 0 sy sz)
(amp-func (circular-amp MFD))))
(list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez)
(center (+ 1 (* -0.5 sx)) 0 0)
(size 0 (/ dfiber 2) (/ dfiber 2))
;(amp-func (circular-amp MFD))
))))
(set! symmetries (list (make mirror-sym (direction Y)) (make mirror-sym
(direction Z) (phase -1))))
(set! pml-layers (list (make pml (thickness 1))))
(set-param! resolution 14)
(define-param nfreq 1) ; number of frequencies at which to compute
flux
(run-until 150
(at-beginning output-epsilon)
(at-end
(to-appended "left1" (in-volume (volume (center (- 0.048 (- l (* 0.5 sx))) 0
0) (size 0 sc sc))
(synchronized-magnetic output-tot-pwr)))
(at-end output-dpwr)
))
Yinlan
University Research Fellow
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
