<!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">
The support of your basis function shouldn't be 'All'. It has to be a
region, even if this region includes everything.<br>
<br>
FunctionSpace {<br>
{ Name Hgrad; Type Form0;<br>
BasisFunction {<br>
{ Name sn; NameOfCoef un; Function BF_Node; <br>
Support Region[Omega]; Entity NodesOf[All]; }<br>
//{ Name sn2; NameOfCoef un2; Function BF_Node_2E; <br>
// Support Region[Omega]; Entity EdgesOf[All]; }<br>
}<br>
Constraint {<br>
{ NameOfCoef un; EntityType NodesOf ; NameOfConstraint Dirichlet;
}<br>
{ NameOfCoef un; EntityType NodesOf ; NameOfConstraint Bloch; }<br>
}<br>
}<br>
}<br>
<br>
There is still a problem with the periodicity BCs. <br>
If you do not use that constraint, Bloch, you get a result.<br>
<br>
Regards,<br>
Ruth<br>
<br>
<br>
<br>
On 09/02/10 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>
<pre class="moz-signature" cols="72">--
Dr. Ir. Ruth V. Sabariego
University of Liege, Dept. of Electrical Engineering & Computer Science,
Applied & Computational Electromagnetics (ACE),
phone: +32-4-3663737 - fax: +32-4-3662910 - <a class="moz-txt-link-freetext" href="http://ace.montefiore.ulg.ac.be/">http://ace.montefiore.ulg.ac.be/</a>
</pre>
</body>
</html>