<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Dear Guillaume,<br>
<br>
The 'Region[All]' is not supported as a support of a BasisFunction.
Sorry for that...<br>
<br>
You should define it explicitly, with<br>
<br>
Support Region[Omega];<br>
<br>
instead of<br>
<br>
Support Region[All];<br>
<br>
Thank you for reporting this bug. It will be corrected.<br>
<br>
Also, the constraint of type 'LinkCplx' asks for an exact
correspondence of nodes between the surfaces to be linked. You can
specify this in your .geo file with:<br>
<br>
Transfinite Line {1,4} = (14e-006-10.8e-006)/lc_air + 1;<br>
Transfinite Line {2,8} = (10.8e-006-3.1e-006)/lc_air + 1;<br>
Transfinite Line {3,12} = (3.1e-006-1e-007)/lc_air + 1;<br>
<br>
... and do not forget to correct the Function that applies Region
SurfBlochRight on SurfBlochLeft, i.e.<br>
<br>
Function Vector[$X-period,$Y,$Z] ;<br>
<br>
instead of<br>
<br>
Function Vector[$X+period,$Y,$Z] ;<br>
<br>
<br>
Regards,<br>
Patrick<br>
<br>
<br>
<br>
On 09/02/2010 15:35, <a class="moz-txt-link-abbreviated" href="mailto:gdemesy@physics.utoronto.ca">gdemesy@physics.utoronto.ca</a> wrote:
<blockquote
cite="mid:20100209093522.r8bayi16ucosgg0g@webmail.physics.utoronto.ca"
type="cite">Dear Getdp List,
<br>
<br>
I am trying to calculate the field diffracted by a grating enlighten by
a plane wave in the visible range (for 1 polarization case = 2D pb, in
harmonic regime). In short, it amounts to solve a scalar Helmholtz
equation together with Bloch conditions, PMLs and an appropriate source
term located inside the diffractive element (just like we did in
<a class="moz-txt-link-abbreviated" href="http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-15-26-18089">www.opticsinfobase.org/oe/abstract.cfm?uri=oe-15-26-18089</a>).
<br>
<br>
I have been struggling for a while on the rather simple attached
problem. This gmsh/getdp code is a "benchmark problem" for which I
already have the solution via another FEM code. I am stuck with this
"Segmentation violation (invalid memory reference)" at Pre-processing
step.
<br>
<br>
I have unsuccessfully tried to simplify each step of the model as much
as possible:
<br>
- a geometry generated whithout using "Extrude"
<br>
- a rough mesh
<br>
- Neumann instead of Bloch conditions
<br>
- source = 0
<br>
...
<br>
My conclusion: I'm missing something elementary somewhere.
<br>
Any hint would be greatly appreciated...
<br>
<br>
Thanks a lot for your time.
<br>
<br>
Best,
<br>
Guillaume Demésy,
<br>
Postdoc at University of Toronto,
<br>
Sajeev John's Group
<br>
<br>
----------------------------------------------------------------
<br>
This message was sent using IMP, the Internet Messaging Program.
<br>
<br>
<pre wrap="">
<fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
getdp mailing list
<a class="moz-txt-link-abbreviated" href="mailto:getdp@geuz.org">getdp@geuz.org</a>
<a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/getdp">http://www.geuz.org/mailman/listinfo/getdp</a>
</pre>
</blockquote>
<br>
</body>
</html>