[Gmsh] Forcing all jacobians to be positive

[Christophe Geuzaine]

> 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.

Christophe




> 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 
http://www.montefiore.ulg.ac.be/~geuzaine