<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.E-MailFormatvorlage17
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi Peter,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">there are some issues in your coil.pro.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">E.g. the crash comes due to this line:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">  Omega = Region[ All ];<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">You have to specify the regions by their names:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">  Omega = Region[ {Omega_air, Omega_coil} ];<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Also your function space is of form 0 (i.e. nodal<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">based scalar basis functions). For the vector<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">potential you need form 1 (edge based vector<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">basis functions). In your case of a 2D simulation<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">you need Form1P.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Maybe the best starting point would be to<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">modify the magnetic example delivered with<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">GetDP. You will find it in your installation directory:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">C:\Path\to\GetDP\demos\magnet.pro<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">There are two resolutions: One calculates the<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">magnetic vector potential. Alternatively the other<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">calculates the magnetic scalar potential.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">You can also look into the examples at:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><a href="http://onelab.info/wiki/Inductor">http://onelab.info/wiki/Inductor</a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Cheers,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Michael<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="DE" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span lang="DE" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> getdp [mailto:getdp-bounces@ace20.montefiore.ulg.ac.be]
<b>On Behalf Of </b>Peter Kis<br>
<b>Sent:</b> Friday, August 07, 2015 12:27 AM<br>
<b>To:</b> getdp@geuz.org<br>
<b>Subject:</b> [Getdp] [solve] Segmentation fault (core dumped)<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Dear List Members,<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I am at the very beginning of learning how to use getdp. In fact I bumped into a problem at my first model, which is about an axisymmetric 2D coil in time domain. I use vector potential formalism with excitation
 current perpendicular to the plane. Therefore this is a nodal FEM problem because the vector potential has only phi component as the current density.
<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I did my best to tackle this problem but something is wrong with my *.pro file. However the example files run well.
<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I use getdp version 2.6.1 64 bit on Linux and get this message:<br>
<br>
<span style="font-family:"Courier New"">getdp <a href="http://coil.pro">coil.pro</a> -solve '#1' -pos '#1' -msh coil.msh<br>
Info    : Running 'getdp <a href="http://coil.pro">coil.pro</a> -solve #1 -pos #1 -msh coil.msh' [GetDP 2.6.1, 1 node]<br>
Info    : Started (Fri Aug  7 00:13:51 2015, CPU = 0.076s, Mem = 255.719Mb)<br>
Info    : Initializing Gmsh<br>
Info    : Loading problem definition '<a href="http://coil.pro">coil.pro</a>'<br>
Info    : Loading problem definition 'param.geo'<br>
Info    : Selected Resolution 'Solution'<br>
Info    : Loading Geometric data 'coil.msh'<br>
Info    : System 'Syst' : Real<br>
P r e - P r o c e s s i n g . . .<br>
Info    : Treatment Formulation 'VectorPotentialFormulation'<br>
make: *** [solve] Segmentation fault (core dumped)</span><br>
<br>
<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Could you take a look at my project file what's wrong with it?  This is really my very first project in getDP. So any clue is welcome.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Thank you for your help in advance.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Peter<br>
<br>
-----------<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="http://coil.pro">coil.pro</a>:<br>
------------<br>
<span style="font-family:"Courier New"">// <a href="http://coil.pro">coil.pro</a>: weak formulation<br>
<br>
Include "param.geo";<br>
<br>
Group{ // OK<br>
    Omega_air = Region[{Omega_air}];<br>
    Omega_coil = Region[{Omega_coil}];<br>
    Omega = Region[ All ];<br>
    Gama_D = Region[{Gama_D}]; // Dirichlet for A=0<br>
}<br>
<br>
Function{ // OK<br>
    InitialState[Omega] = 0.;<br>
    mu0 = 4e-7*Pi;<br>
    Freq = 50;<br>
    t0 = 0;<br>
    tmax = 2.;<br>
    dt = 0.1;<br>
    theta = 1.; // backward Euler<br>
    source_current_t = {0, 1, 2};<br>
    source_current_j = {0, 1, 1};<br>
    Js[] = InterpolationLinear[$1]{ListAlt[source_current_t, source_current_j]};<br>
//    Jext[] = Js[$Time];<br>
    Jext[] = 50*(1-Exp[-0.5*$Time]);<br>
    sigma[Omega_air] = 1;<br>
    sigma[Omega_coil] = 6e7;<br>
}<br>
<br>
Jacobian { // OK<br>
 { Name JVol; // for 2D regions<br>
   Case {<br>
     {Region All; Jacobian VolAxi;}<br>
   }<br>
 }<br>
 { Name JSur; // for 1D regions<br>
   Case {<br>
     {Region All; Jacobian SurAxi;}<br>
   }<br>
 }<br>
}<br>
<br>
Integration {<br>
  { Name I1 ;<br>
    Case {<br>
      { Type Gauss ;<br>
        Case {<br>
          { GeoElement Point ; NumberOfPoints  1 ; }<br>
          { GeoElement Line ; NumberOfPoints  4 ; }<br>
          { GeoElement Triangle ; NumberOfPoints  6 ; }<br>
          { GeoElement Quadrangle ; NumberOfPoints 7 ; }<br>
          { GeoElement Tetrahedron ; NumberOfPoints 15 ; }<br>
          { GeoElement Hexahedron ; NumberOfPoints 34 ; }<br>
        }<br>
      }<br>
    }<br>
  }<br>
}<br>
<br>
Constraint{ // OK<br>
    // Initial State<br>
    {Name InitialData; Type Init;<br>
        Case {<br>
            {Region Omega; Value InitialState[];}<br>
        }<br>
    }<br>
    // Dirichlet in Inf.<br>
    {Name DirichletBC; Type Assign;<br>
        Case {<br>
            {Region Gama_D; Value 0;}<br>
        }<br>
    }<br>
    // No constraint for the axis<br>
}<br>
<br>
FunctionSpace{ // OK<br>
  { Name Afs; Type Form0;<br>
    BasisFunction{<br>
      {Name Ni; NameOfCoef Ai; Function BF_Node;<br>
       Support Omega; Entity NodesOf[All];}<br>
    }<br>
    Constraint {<br>
        {NameOfCoef Ai; EntityType NodesOf;<br>
         NameOfConstraint DirichletBC;}<br>
        {NameOfCoef Ai; EntityType NodesOf;<br>
         NameOfConstraint InitialData;}<br>
    }<br>
  }<br>
}<br>
<br>
Formulation{ // OK<br>
  {Name VectorPotentialFormulation; Type FemEquation;<br>
   Quantity{<br>
     {Name a; Type Local; NameOfSpace Afs;}<br>
   }<br>
   Equation{<br>
    Galerkin{ [1/mu0*Dof{Curl a},{Curl a}];<br>
      In Omega; Jacobian JVol; Integration I1;}<br>
<br>
    Galerkin{ DtDof[sigma[]*Dof{a},{a}];<br>
      In Omega; Jacobian JVol; Integration I1;}<br>
<br>
    Galerkin{ [-Jext[],{a}];<br>
      In Omega_coil; Jacobian JVol; Integration I1;}<br>
   }<br>
  }<br>
}<br>
<br>
Resolution{<br>
  {Name Solution;<br>
    System{<br>
      {Name Syst; NameOfFormulation VectorPotentialFormulation;}<br>
    }<br>
    Operation{<br>
        InitSolution[Syst];SaveSolution[Syst];<br>
        TimeLoopTheta[t0,tmax,dt,theta] {<br>
            Generate[Syst]; Solve[Syst]; SaveSolution[Syst];<br>
        }<br>
    }<br>
  }<br>
}<br>
<br>
PostProcessing{<br>
  {Name VectorPotentialPostProc; NameOfFormulation VectorPotentialFormulation;<br>
    Quantity {<br>
      {Name vec_pot; Value {Local{[{a}];In Omega;Jacobian JVol;}}}<br>
    }<br>
  }<br>
}<br>
<br>
PostOperation{<br>
  {Name Map_Vectorpotential; NameOfPostProcessing VectorPotentialPostProc;<br>
    Operation{<br>
      Print[vec_pot, OnElementsOf Omega, File "A.pos"];<br>
    }<br>
  }<br>
}<br>
</span><br>
---------------------<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">param.geo<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">---------------------<br>
<span style="font-family:"Courier New"">// param.geo: the auxiliary file<br>
<br>
// PhysicalIDs<br>
Omega_coil = 700;<br>
Omega_air = 1000;<br>
Gama_D = 500;<br>
<br>
// Geometry<br>
r_enc = 60e-3;<br>
z_enc = 40e-3;<br>
inner_radius = 59e-3/2;<br>
thick = 100e-6;<br>
width = 4e-3;<br>
<br>
// Mesh<br>
lc_air = 0.008;<br>
lc_coil = 0.5e-3;<br>
</span><br>
<br>
---------------------<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">coil.geo<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">---------------------<br>
<br>
<span style="font-family:"Courier New"">// coil.geo: creation of the geometry with GMSH<br>
// <br>
<br>
Include "param.geo";<br>
<br>
// Enclosure points<br>
p1 = newp; Point(p1) = {0,-z_enc,0,lc_air};<br>
p2 = newp; Point(p2) = {r_enc,-z_enc,0,lc_air};<br>
p3 = newp; Point(p3) = {r_enc,z_enc,0,lc_air};<br>
p4 = newp; Point(p4) = {0,z_enc,0,lc_air};<br>
<br>
// First turn points<br>
p5 = newp; Point(p5) = {inner_radius,-width/2,0,lc_coil};<br>
p6 = newp; Point(p6) = {inner_radius,width/2,0,lc_coil};<br>
p7 = newp; Point(p7) = {inner_radius+thick,width/2,0,lc_coil};<br>
p8 = newp; Point(p8) = {inner_radius+thick,-width/2,0,lc_coil};<br>
<br>
// Enclosure lines<br>
L1 = newreg; Line(L1) = {p1,p2};<br>
L2 = newreg; Line(L2) = {p2,p3};<br>
L3 = newreg; Line(L3) = {p3,p4};<br>
L4 = newreg; Line(L4) = {p4,p1};<br>
<br>
// First turn lines<br>
L5 = newreg; Line(L5) = {p5,</span>p6};<br>
<span style="font-family:"Courier New"">L6 = newreg; Line(L6) = {p6,p7};<br>
L7 = newreg; Line(L7) = {p7,p8};<br>
L8 = newreg; Line(L8) = {p8,p5};<br>
<br>
// Line loops<br>
Bound = newreg; Line Loop(Bound) = {L1,L2,L3,L4};<br>
Turn = newreg; Line Loop(Turn) = {L5,L6,L7,L8};<br>
<br>
SurfaceOmega_air = newreg; Plane Surface(SurfaceOmega_air) = {Bound,Turn};<br>
SurfaceOmega_coil = newreg; Plane Surface(SurfaceOmega_coil) = {Turn};<br>
<br>
Physical Surface(Omega_air) = {SurfaceOmega_air};<br>
Physical Surface(Omega_coil) = {SurfaceOmega_coil};<br>
<br>
// For Dirichlet bnd cnd<br>
Physical Line(Gama_D) = {L1,L2,L3};<br>
<br>
Mesh.Algorithm = 1;<br>
//Mesh.SurfaceFaces = 1;<br>
//Mesh.SurfaceNumbers = 1;<br>
//Mesh.LabelType = 2; // Show Phys ID<br>
//Mesh 2;</span><o:p></o:p></p>
</div>
</div>
</div>
</body>
</html>