[Gmsh] problem with gmsh to Freefem mesh?

Geordie McBain gdmcbain at freeshell.org
Thu May 7 02:02:59 CEST 2015


2015-05-06 22:53 GMT+10:00 Julien Derr <julien.derr at gmail.com>:
> Hi Georgie,
>
> here is a minimum example: if you run the script.sh you would get the error,

Hello.  No, actually, it seemed to work fine here.  I wasn't sure what
freefem.edp was supposed to contain so wrote a simple one solving a
Dirichlet problem for the oblong island using P1 triangular elements.
I think that that should fail if any elements had nonpositive
Jacobians.

mesh Th = readmesh("freefem_mesh.mesh");
fespace Vh(Th, P1);
Vh u, v;
solve laplace(u, v) =
  int2d(Th)(dx(v) * dx(u) + dy(v) * dy(u))
  + on(27, u=1) + on(16, 17, 18, 19, u=0);
plot(u, ps="u.eps");

My versions are Gmsh 2.8.4 and FreeFem++ 3.340000 (10 février 2015),
both compiled locally from SVN/hg sources on 64-bit Debian testing.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: u.png
Type: image/png
Size: 43544 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20150507/66a345bc/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: freefem.edp
Type: application/octet-stream
Size: 199 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20150507/66a345bc/attachment.edp>


More information about the gmsh mailing list