rand mouradi wrote: > Dear Dr. Steven and meep users, > I'm a new meep user, and I'm trying to run my program after I defined > a function and trying to get the integral of it. But I'm having the > following error: > ERROR: In procedure /: > ERROR: Wrong type argument in position 1: #(-0.05 -0.45 -0.05) > I'm not sure if I need to initial anything, I tried to look for > solution but nothing worked. > also I'm not sure if the integral can be defined over an area (2-D) > instead of volume (3-D). > Please let me know what is wrong if you have any idea. I really > appreciate your help. > > here is my program: > (reset-meep) > ;The dielectric parameters > > (define-param feps 5.5) ;fat permittivity,epsilon,at 900 MHz > > > (define-param fseg 3.493) ;fat D-conductivity, segma, at 900 MHz > > (define-param ro 916e-6) ;Tissue density in Kg/cm3 > > ;the cell dimentions > > (define-param fw 3) ; fat thickness in cm > > (define-param dpml 1) ; PML layer thickness > (define-param sx 8) ; size of cell in x direction > (define-param sy 16) ; size of cell in y direction > > (init-fields) > (define (sar fseg ro ez) > (* (/ fseg (* 0.5 ro))(+(sqr (ez))(* hz 0)))) > > (set! geometry-lattice (make lattice (size sx sy no-size))) > > > (set! geometry > (list > (make block (center 0 0) (size infinity fw infinity) > (material (make dielectric (epsilon feps) (D-conductivity > fseg)))))) > > > (set! sources (list > (make source(src (make continuous-src > (wavelength 25.5))) > (component Ez) > (center 0 0) (size 0 3)))) > > (set! pml-layers (list (make pml (thickness dpml)))) > > (set! resolution 10) > > > (print "The integral of our sar function from (0,-1.5,0) to (0,1.5,0) > is: " > (integrate-field-function (list Ez Hz) sar (volume (size 0 3 0) > (center 0 0 0))) "\n") > > (output-field-function "sar-function" (list Ez Hz) sar) > > > (define (my-sar-output) (output-field-function "sar-function" (list Ez > Hz) sar)) > > (run-until 255 (at-every 1 my-sar-output)) > ------------------------------------------------------------------------ > > _______________________________________________ > meep-discuss mailing list > meep-discuss@ab-initio.mit.edu > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss Hi, rand mouradi (define (sar fseg ro ez) (* (/ fseg (* 0.5 ro))(+(sqr (ez))(* hz 0)))) ============================================= This function looks weired. What do you mean? hz*0 is always equal to 0 and ez is not a function.
Zheng 2008-11-20
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss