<div dir="ltr">Thanks for the answer, Lars! That solved the problem.<br><br>That's right, I'm starting to work with 3D photonic crystal structures and I believe GetDP is<br>general enough to allow me to produce the band diagrams.<br>
However I'm not sure how to formulate the problem.<br>Since you are also working around the same subject, could you point me in the right direction?<br><br><span style="font-family: courier new,monospace;">Constraint {</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">   { Name BoundCond; Type Link;</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">     Case {</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">       { Region fxp; RegionRef fxm; Coefficient Complex[Cos[kx*ax], -Sin[kx*ax]]; Function Vector[$X-ax,$Y,$Z]; }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">       { Region fyp; RegionRef fym; Coefficient Complex[Cos[ky*ay], -Sin[ky*ay]]; Function Vector[$X,$Y-ay,$Z]; }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">       { Region fzp; RegionRef fzm; Coefficient Complex[Cos[kz*az], -Sin[kz*az]]; Function Vector[$X,$Y,$Z-az]; }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">     }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">   }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">}</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">FunctionSpace {</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">  { Name Curl_3D; Type Form1;</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">    BasisFunction {</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">      { Name W1; NameOfCoef e1; Function BF_Edge; Support Domain; Entity EdgesOf[All]; }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">    }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">    Constraint {</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">      { NameOfCoef e1; EntityType EdgesOf; NameOfConstraint BoundCond; }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">    }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">  }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">}</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Formulation {</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">  { Name EMW_3D_Eigen; Type FemEquation;</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">    Quantity {</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">      { Name e; Type Local; NameOfSpace Curl_3D; }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">    }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">    Equation {</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">      Galerkin { [ -Dof{Curl e} , {Curl e} ];</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">                 In Domain; Jacobian JVol; Integration I2; }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">      Galerkin { DtDt [ mu[] * eps[] * Dof{e} , {e} ];</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">                 In Domain; Jacobian JVol; Integration I2; }</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">    }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">  }</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">}</span><br style="font-family: courier new,monospace;">
<br><br>Thank you for your help,<br>Lucas<br><br><br><br>On Wed, Jul 30, 2008 at 7:13 AM, Lars Rindorf <<a href="mailto:Lars.Rindorf@teknologisk.dk">Lars.Rindorf@teknologisk.dk</a>> wrote:<br>><br>> Hi Lucas<br>
>  <br>> You need to tell getdp the specific name of your domain. I.e. you should set<br>> Group {<br>>   Domain = Region[4000];<br>> }<br>> instead of<br>> Group {<br>>   Domain = Region[All];<br>
> }<br>>  <br>> I'm also working on 3D electromagnetic bloch simulations, and I can help you if have questions.<br>>  <br>> Are you work on band structure calc of photonic crystals?<br>>  <br>> Kind regards<br>
> Lars<br>>  <br>> ________________________________<br>> Fra: <a href="mailto:getdp-bounces@geuz.org">getdp-bounces@geuz.org</a> [mailto:<a href="mailto:getdp-bounces@geuz.org">getdp-bounces@geuz.org</a>] På vegne af Lucas<br>
> Sendt: 29. juli 2008 21:51<br>> Til: <a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>> Emne: [Getdp] Electromagnetic 3D model - Crash info (bug?)<br>><br>> Hi all,<br>><br>> I've just started working with GetPD and I'm not sure I'm doing the right thing here.<br>
> The attached files should solve an 3D eigenvalue EM problem (no boundary conditions or<br>> sources yet) but I can't seem to make them work:<br>><br>> getdp test -solve EMW_3D_Eigen -pos EMW_3D_Eigen<br>
><br>> Loading   : Problem definition 'em\<a href="http://test.pro">test.pro</a>'<br>> Loading   : Problem definition 'em\<a href="http://em.pro">em.pro</a>'<br>> Info      : Selected Resolution 'EMW_3D_Eigen'<br>
> Loading   : Geometric data 'em\test.msh'<br>> Info      : System 'A' : Complex, Frequency = 0 Hz<br>> P r e - P r o c e s s i n g . . .<br>> Operation : Treatment Formulation 'EMW_3D_Eigen'<br>
> Error     : Segmentation violation (invalid memory reference)<br>>             ------------------------------------------------------<br>>             You have discovered a bug in GetDP! You may report it<br>>             by e-mail (together with any helpful data permitting to<br>
>             reproduce it) to <<a href="mailto:getdp@geuz.org">getdp@geuz.org</a>><br>><br>> The geometry file is transfinite just cause I plan to use Bloch periodic conditions later (any<br>> help in that matter would also be welcome).<br>
><br>> If I can manage to understand how to correctly create a model in GetDP, I would say it's one<br>> of the best sw packages for FEM, but I couldn't do it following only the documentation.<br>><br>
> Thanks for any help,<br>> Lucas<br><br></div>