Hi, all,

Sorry I have to correct my previous email for help.

The meep version in our server should be 0.20.3-1 directly from the following
website
http://packages.debian.org/testing/science/meep.

Thus it produced the errors as I described in the following email. I am not sure
if this is the reason of the errors. Do you have any idea how to update the
version? In debina, only old version 0.20.3-1 available. We hope installation of
the latest version can avoid this problem. Thanks a lot.

Cheers

Yinlan




----- Forwarded message from Yinlan Ruan <[email protected]> -----
    Date: Tue, 30 Jun 2009 00:50:22 +0930
    From: Yinlan Ruan <[email protected]>
Reply-To: Yinlan Ruan <[email protected]>
 Subject: [Meep-discuss] power output problem with MPI MEEP 1.0
      To: [email protected]

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

Reply via email to