Hi Chris,
The error you are having is due to the flux region calculated in the
holey-wvg example.
When you define a flux region in 3D you have to define also the normal
direction of your flux-plane (if it is ambiguous) as you can read in Meep
reference.
http://ab-initio.mit.edu/wiki/index.php/Meep_Reference#flux-region
You just have to add to the flux definition '(direction Y)' for instance and
that's all.
Regards.
David Domenech Gomez
Optical & Quantum Communications Group
Universidad Politecnica de Valencia
www.gco.upv.es
----------------------------------------------------------------------
Message: 1
Date: Fri, 18 May 2007 11:08:42 -0700
From: Chris <[EMAIL PROTECTED]>
Subject: Re: [Meep-discuss] meep: Could not determine normal direction
for given volume.
Greetings,
Does anyone have an example of a libctl file for a meep calcualtion on
a 3d ( preferably laminar ) structure ?
Thanks,
Chris
On May 15, 2007, at 9:35 AM, Chris wrote:
> Greetings,
>
> I am trying to extend the "holey waveguide" example to 3D ( as part of
> an incremental effort to model more complex structures. )
>
> When I try to run meep on the new 3D structure, I get the message
> "meep: Could not determine normal direction for given volume."
>
> The code for defining the structure is as follows :
>
> (define sx (+ (* 2 (+ pad dpml N)) d -1)) ; size of cell in x direction
> (define lx (+ (* 2 (+ pad dpml N)) d -1)) ; size of cell in x direction
> (set! geometry-lattice (make lattice (size sx sy sz)))
>
> (set! geometry
> (append ; combine lists of objects:
> (list (make block (center 0 0) (size lx w 0.3)
> (material (make dielectric (epsilon eps)))))
> (geometric-object-duplicates (vector3 1 0) 0 (- N 1)
> (make cylinder (center (/ d 2) 0) (radius r) (height 0.3)
> (material air)))
> (geometric-object-duplicates (vector3 -1 0) 0 (- N 1)
> (make cylinder (center (/ d -2) 0) (radius r) (height 0.3)
> (material air)))))
>
> (set! pml-layers (list (make pml (thickness dpml))))
>
> I am quite new to meep ( and to libctl and guile for that matter ) and
> would appreciate any help with resolving this.
>
> Regards,
>
> Chris
>
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
