[Gmsh] Forcing all jacobians to be positive

Christophe Geuzaine cgeuzaine at uliege.be
Sat May 30 14:48:11 CEST 2020

> On 29 May 2020, at 20:48, Carole-Anne Daunais <carole-anne.daunais at polymtl.ca> wrote:
> Hi,
> I am using GMSH 4.5.6 to generate a 3d mesh. I used the following 2 lines in the .geo file to ensure that no element would have a jac.<0.
> Mesh.ElementOrder = 2;
> Mesh.HighOrderOptimize = 1;
> It unfortunately didn't work for me. The mesh is generated, but the software failed (froze). Is there another way I can force all jacobians to be positive when creating my mesh?

It works fine here. However you should remove "Mesh.SubdivisionAlgorithm = 2;" (and set the sizes you want directly), as this makes it harder for the optimizer.

You can also use "Mesh.HighOrderOptimize = 2;", which will apply a global elastic smoothing step before the brute-force optimization.


> Thank you,
> Carole-Anne
> <test2.geo>_______________________________________________
> gmsh mailing list
> gmsh at onelab.info
> http://onelab.info/mailman/listinfo/gmsh

Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 

More information about the gmsh mailing list