Hi Jian, The problem is this line: (set! default-material (make material-function (epsilon-func feps)))
Since you are using epsilon-func and not material-func, you must specify an actual material. Changing the line to this makes it work: (set! default-material (make dielectric (epsilon-func feps))) Cheers, Matt Miskieiwcz On Mon, Aug 2, 2010 at 11:21 AM, zhang jian <[email protected]> wrote: > Hi, dear Meep users, I am a beginner of Meep. I write a simple meep > code, and run it at the server of my department. However it can not > calculate a final result, though it is just a simple code. Could you > please tell me how to simply this code? Thank you very much.~ > > Initializing structure... > Working in 2D dimensions. > block, center = (-5,5,0) > size (10,1,1e+20) > axes (1,0,0), (0,1,0), (0,0,1) > dielectric constant epsilon = 11.9 > block, center = (2,6.5,0) > size (10,1,1e+20) > axes (1,0,0), (0,1,0), (0,0,1) > dielectric constant epsilon = 11.9 > subpixel-averaging is 19.964% done, 18.708 s remaining > subpixel-averaging is 39.948% done, 7.03004 s remaining > subpixel-averaging is 59.9321% done, 3.10042 s remaining > subpixel-averaging is 79.9161% done, 1.16184 s remaining > subpixel-averaging is 99.9001% done, 0.0046283 s remaining > subpixel-averaging is 19.964% done, 18.5258 s remaining > subpixel-averaging is 39.948% done, 6.83587 s remaining > subpixel-averaging is 59.9321% done, 3.094 s remaining > subpixel-averaging is 79.9161% done, 1.16121 s remaining > subpixel-averaging is 99.9001% done, 0.00463442 s remaining > subpixel-averaging is 13.9447% done, 28.105 s remaining > subpixel-averaging is 27.9034% done, 10.931 s remaining > subpixel-averaging is 41.8621% done, 6.41271 s remaining > subpixel-averaging is 55.8208% done, 3.39532 s remaining > subpixel-averaging is 69.7795% done, 1.92881 s remaining > subpixel-averaging is 83.7381% done, 0.868258 s remaining > subpixel-averaging is 97.6968% done, 0.102055 s remaining > subpixel-averaging is 13.9447% done, 28.8167 s remaining > subpixel-averaging is 27.9034% done, 12.0302 s remaining > subpixel-averaging is 41.8621% done, 6.41883 s remaining > subpixel-averaging is 55.8208% done, 3.66837 s remaining > subpixel-averaging is 69.7795% done, 2.00258 s remaining > subpixel-averaging is 83.7381% done, 0.89728 s remaining > subpixel-averaging is 97.6968% done, 0.109034 s remaining > subpixel-averaging is 1.20386% done, 359.114 s remaining > subpixel-averaging is 2.40893% done, 177.411 s remaining > subpixel-averaging is 3.61399% done, 115.726 s remaining > subpixel-averaging is 4.81906% done, 85.4739 s remaining > subpixel-averaging is 6.02413% done, 66.8641 s remaining > subpixel-averaging is 7.22919% done, 55.6286 s remaining > .......................................... > > > My code: > (reset-meep) > (define eps-si 11.9) > (define eps-sio2 (expt 1.46 2)) > (define wr 1) > (define bendr 5) > (define bend-xcen 0) > (define bend-ycen 0) > (define gap 0.5) > (define (circ p) ( > + (expt (- (vector3-x p) bend-xcen) 2) (expt (- > (vector3-y p) bend-ycen) 2) > ) > ) > (define (fline1 p) (- (vector3-x p) bend-xcen)) > (define (fline2 p) (- (vector3-y p) bend-ycen)) > (define r1 (+ bendr (/ wr 2))) > (define r2 (- bendr (/ wr 2))) > > (define (feps p) > (if > (and (<= (circ p) (expt r1 2)) > (>= (circ p) (expt r2 2)) > (>= (fline1 p) 0) > (>= (fline2 p) 0) > ) > eps-si eps-sio2 > ) > ) > > (set! geometry-lattice (make lattice (size 22 16 no-size))) > > (set! default-material (make material-function (epsilon-func feps))) > > (set! geometry > (append > (list > (make block (center -5 bendr) (size 10 wr infinity) > (material (make dielectric (epsilon eps-si)))) > (make block (center 2 (+ bendr (+ wr gap))) (size 10 wr > infinity) > (material (make dielectric (epsilon eps-si)))) > )) > ) > > (set! pml-layers (list (make pml (thickness 1.0)))) > (set! resolution 10) > > (set! sources (list > (make source > (src (make continuous-src > (wavelength (* 2 (sqrt 12))) (width 20))) > (component Ez) > (center -10 bendr) (size 0 1)))) > (run-until 200 > (at-beginning output-epsilon) > (to-appended "pbs" (at-every 0.6 output-efield-z))) > > > Jian > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss >
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
