<div dir="ltr"><div><div><div>Thanks again Geordie, <br><br></div>I forgot the freefem file, but it was just <br>mesh Th = readmesh("freefem_mesh.mesh");<br><br></div>so It seems the problem commes from my version of gmsh, I am still at 2.5<br></div>I'll update it.<br><div><div><br></div><div>or else Mathias mentioned optimize3D, what is the syntax/use of this option ? <br><br><br></div><div>Julien<br></div><div>************************************************************************<br><p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)">Hi,</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)"> </span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)" lang="EN-US">I
 didn’t run your example, but how about Mesh > Optimize3D? This has 
solved problems with malformed meshes for me in the past.</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)" lang="EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)" lang="EN-US">HTH,</span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)" lang="EN-US"> </span></p>
<span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)" lang="EN-US">Matthias</span><br></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 7, 2015 at 2:02 AM, Geordie McBain <span dir="ltr"><<a href="mailto:gdmcbain@freeshell.org" target="_blank">gdmcbain@freeshell.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">2015-05-06 22:53 GMT+10:00 Julien Derr <<a href="mailto:julien.derr@gmail.com">julien.derr@gmail.com</a>>:<br>
> Hi Georgie,<br>
><br>
> here is a minimum example: if you run the script.sh you would get the error,<br>
<br>
</span>Hello.  No, actually, it seemed to work fine here.  I wasn't sure what<br>
freefem.edp was supposed to contain so wrote a simple one solving a<br>
Dirichlet problem for the oblong island using P1 triangular elements.<br>
I think that that should fail if any elements had nonpositive<br>
Jacobians.<br>
<br>
mesh Th = readmesh("freefem_mesh.mesh");<br>
fespace Vh(Th, P1);<br>
Vh u, v;<br>
solve laplace(u, v) =<br>
  int2d(Th)(dx(v) * dx(u) + dy(v) * dy(u))<br>
  + on(27, u=1) + on(16, 17, 18, 19, u=0);<br>
plot(u, ps="u.eps");<br>
<br>
My versions are Gmsh 2.8.4 and FreeFem++ 3.340000 (10 février 2015),<br>
both compiled locally from SVN/hg sources on 64-bit Debian testing.<br>
</blockquote></div><br></div>