[Gmsh] diffrence in mesh with inline command on the same code

[Benoit Pouliot]

Hi,

I have a weird behaviour with the mesh generation.
I do not know if you can duplicate it, but here anyway my problem.


I run gmsh with a 3D mesh 2 times and check the difference between the .msh.
I run it on the same code (2 different named files), but sometimes, I get 2
different .msh and sometime not...

For example, if I send the screen output in the /dev/null, I got 2
different .msh
If I send it in a out1.txt and a out2.txt, then I got the same .msh!

The mesh and the differences are always the same when I retry it.


In attach files, you can see my tests in the bash script "exec".
(I do not have any .opt file in my test folder. I do not have .gmsh-options
file in my $HOME)



Is there any way to force gmsh to give always the same result with inline
command generation?

I can understand that the generating mesh algorithms can give different
results sometimes.
(I can see that when I open the gmsh api and it is normal.)
But with the inline command, it would be useful that this kind of behaviour
did not appear.

I mean, it is weird that the way I send the output change things for gmsh...

Thanks.


-- 
Benoît Pouliot - Développeur locigiel Bodycad
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0001.html>
-------------- next part --------------
Info    : Running '/home/bpouliot/storage/Gmsh_exe/gmsh-3.0.4-Linux/bin/gmsh cylindre.tet-pri-hex.quad.geo -3 -o cylindre.tet-pri-hex.quad.msh -v 99' [Gmsh 3.0.4, 1 node, max. 1 thread]
Info    : Started on Mon Aug  7 08:18:53 2017
Debug   : Destroying model 
Info    : Reading 'cylindre.tet-pri-hex.quad.geo'...
Debug   : Extrude Surface 201
Debug   : Extrude Curve 1
Debug   : Extrude Point 1
Debug   : Extrude Point 2
Debug   : Extrude Curve 5
Debug   : Extrude Point 2
Debug   : Extrude Point 5
Debug   : Extrude Curve 6
Debug   : Extrude Point 5
Debug   : Extrude Point 4
Debug   : Extrude Curve 4
Debug   : Extrude Point 4
Debug   : Extrude Point 1
Debug   : Removed 40 duplicate points
Debug   : Removed 16 duplicate curves
Debug   : Extrude Surface 202
Debug   : Extrude Curve 2
Debug   : Extrude Point 2
Debug   : Extrude Point 3
Debug   : Extrude Curve 3
Debug   : Extrude Point 3
Debug   : Extrude Point 4
Debug   : Extrude Curve -6
Debug   : Extrude Point 4
Debug   : Extrude Point 5
Debug   : Extrude Curve -5
Debug   : Extrude Point 5
Debug   : Extrude Point 2
Debug   : Removed 43 duplicate points
Debug   : Removed 26 duplicate curves
Debug   : Removed 2 duplicate faces
Debug   : Extrude Surface 202
Debug   : Extrude Curve 2
Debug   : Extrude Point 2
Debug   : Extrude Point 3
Debug   : Extrude Curve 3
Debug   : Extrude Point 3
Debug   : Extrude Point 4
Debug   : Extrude Curve -6
Debug   : Extrude Point 4
Debug   : Extrude Point 5
Debug   : Extrude Curve -5
Debug   : Extrude Point 5
Debug   : Extrude Point 2
Debug   : Removed 40 duplicate points
Debug   : Removed 16 duplicate curves
Debug   : Syncing GEO_Internals with GModel
Debug   : Sync is removing 0 model entities
Debug   : GModel imported:
Debug   : 14 vertices
Debug   : 25 edges
Debug   : 15 faces
Debug   : 3 regions
Info    : Done reading 'cylindre.tet-pri-hex.quad.geo'
Info    : Finalized high order topology of periodic connections
Info    : Meshing 1D...
Info    : Meshing curve 1 (Circle)
Info    : Meshing curve 2 (Circle)
Info    : Meshing curve 3 (Circle)
Info    : Meshing curve 4 (Circle)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 204 (extruded)
Info    : Meshing curve 205 (extruded)
Info    : Meshing curve 206 (extruded)
Info    : Meshing curve 207 (extruded)
Info    : Meshing curve 209 (extruded)
Info    : Meshing curve 210 (extruded)
Info    : Meshing curve 214 (extruded)
Info    : Meshing curve 218 (extruded)
Info    : Meshing curve 226 (extruded)
Info    : Meshing curve 227 (extruded)
Info    : Meshing curve 232 (extruded)
Info    : Meshing curve 248 (extruded)
Info    : Meshing curve 249 (extruded)
Info    : Meshing curve 250 (extruded)
Info    : Meshing curve 251 (extruded)
Info    : Meshing curve 253 (extruded)
Info    : Meshing curve 254 (extruded)
Info    : Meshing curve 258 (extruded)
Info    : Meshing curve 262 (extruded)
Info    : Done meshing 1D (0 s)
Info    : Meshing 2D...
Info    : Meshing surface 201 (Plane, Delaunay)
Debug   : Computing edge loops
Debug   : Generating the mesh
Debug   : Meshing of the convex hull (8 points)
Debug   : Meshing of the convex hull (8 points) done
Debug   : Recovering 4 model Edges
Debug   : Recovering 8 mesh Edges (0 not recovered)
Debug   : Boundary Edges recovered for surface 201
Debug   : Computing mesh size field at mesh vertices 8
Debug   : Delaunizing the initial mesh
Debug   : Starting to add internal points
Debug   : Delaunization of the initial mesh done (0 swaps)
Debug   :       8 points created -- Worst tri radius is    0.900
Info    : Blossom: 11 internal 6 closed
Debug   : Perfect Match Starts 17 edges 10 nodes
Debug   : Perfect Match Succeeded in Quadrangulation (0 sec)
Debug   : 0 two-quadrangles vertices removed
Debug   : 0 diamond quads removed
Info    : Blossom completed (0.00 s):     4 quads     2 triangles 0 invalid quads  0 quads with Q < 0.1 Avg Q = 0.647 Min Q 0.590
Debug   : Type 17 2 triangles generated, 2 internal vertices
Info    : Meshing surface 202 (Plane, Delaunay)
Debug   : Computing edge loops
Debug   : Generating the mesh
Debug   : Meshing of the convex hull (8 points)
Debug   : Meshing of the convex hull (8 points) done
Debug   : Recovering 4 model Edges
Debug   : Recovering 8 mesh Edges (0 not recovered)
Debug   : Boundary Edges recovered for surface 202
Debug   : Computing mesh size field at mesh vertices 8
Debug   : Delaunizing the initial mesh
Debug   : Starting to add internal points
Debug   : Delaunization of the initial mesh done (0 swaps)
Debug   :       8 points created -- Worst tri radius is    0.900
Debug   : Type 17 10 triangles generated, 2 internal vertices
Info    : Meshing surface 211 (extruded)
Info    : Meshing surface 215 (extruded)
Info    : Meshing surface 219 (extruded)
Info    : Meshing surface 223 (extruded)
Info    : Meshing surface 224 (extruded)
Info    : Meshing surface 233 (extruded)
Info    : Meshing surface 237 (extruded)
Info    : Meshing surface 246 (extruded)
Info    : Meshing surface 255 (extruded)
Info    : Meshing surface 259 (extruded)
Info    : Meshing surface 263 (extruded)
Info    : Meshing surface 267 (extruded)
Info    : Meshing surface 268 (extruded)
Info    : Done meshing 2D (0.000316 s)
Info    : Meshing 3D...
Info    : Meshing volume 1 (extruded)
Info    : Meshing volume 2 (extruded)
Info    : Meshing volume 3 (extruded)
Info    : Subdividing extruded mesh
Info    : Swapping 0
Info    : Remeshing surface 255
Info    : Meshing surface 255 (extruded)
Info    : Remeshing surface 259
Info    : Meshing surface 259 (extruded)
Info    : Remeshing surface 263
Info    : Meshing surface 263 (extruded)
Info    : Remeshing surface 267
Info    : Meshing surface 267 (extruded)
Info    : Done meshing 3D (0.000124 s)
Info    : Optimizing 3D mesh...
Info    : Done optimizing 3D mesh (2e-06 s)
Info    : Meshing order 2 (curvilinear on)...
Info    : Meshing curve 1 order 2
Info    : Meshing curve 2 order 2
Info    : Meshing curve 3 order 2
Info    : Meshing curve 4 order 2
Info    : Meshing curve 5 order 2
Info    : Meshing curve 6 order 2
Info    : Meshing curve 204 order 2
Info    : Meshing curve 205 order 2
Info    : Meshing curve 206 order 2
Info    : Meshing curve 207 order 2
Info    : Meshing curve 209 order 2
Info    : Meshing curve 210 order 2
Info    : Meshing curve 214 order 2
Info    : Meshing curve 218 order 2
Info    : Meshing curve 226 order 2
Info    : Meshing curve 227 order 2
Info    : Meshing curve 232 order 2
Info    : Meshing curve 248 order 2
Info    : Meshing curve 249 order 2
Info    : Meshing curve 250 order 2
Info    : Meshing curve 251 order 2
Info    : Meshing curve 253 order 2
Info    : Meshing curve 254 order 2
Info    : Meshing curve 258 order 2
Info    : Meshing curve 262 order 2
Info    : Meshing surface 201 order 2
Info    : Meshing surface 202 order 2
Info    : Meshing surface 211 order 2
Info    : Meshing surface 215 order 2
Info    : Meshing surface 219 order 2
Info    : Meshing surface 223 order 2
Info    : Meshing surface 224 order 2
Info    : Meshing surface 233 order 2
Info    : Meshing surface 237 order 2
Info    : Meshing surface 246 order 2
Info    : Meshing surface 255 order 2
Info    : Meshing surface 259 order 2
Info    : Meshing surface 263 order 2
Info    : Meshing surface 267 order 2
Info    : Meshing surface 268 order 2
Info    : Meshing volume 1 order 2
Info    : Meshing volume 2 order 2
Info    : Meshing volume 3 order 2
Info    : Finalized high order topology of periodic connections
Info    : Surface mesh: worst distortion = 1 (0 elements in ]0, 0.2]); worst gamma = 0.683282
Info    : Volume Mesh: worst distortion = 1 (0 elements in ]0, 0.2])
Info    : Done meshing order 2 (0.00111 s)
Info    : 163 vertices 171 elements
Info    : Writing 'cylindre.tet-pri-hex.quad.msh'...
Info    : Done writing 'cylindre.tet-pri-hex.quad.msh'
Info    : Stopped on Mon Aug  7 08:18:53 2017
Debug   : Destroying model cylindre.tet-pri-hex.quad
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cylindre.tet-pri-hex.quad.geo
Type: application/octet-stream
Size: 1213 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0002.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cylindre.tet-pri-hex.quad.msh
Type: model/mesh
Size: 13285 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0002.msh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: exec
Type: application/octet-stream
Size: 1051 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: file.geo
Type: application/octet-stream
Size: 1213 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0003.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: file.msh
Type: model/mesh
Size: 13287 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0003.msh>
-------------- next part --------------
Info    : Running '/home/bpouliot/storage/Gmsh_exe/gmsh-3.0.4-Linux/bin/gmsh file.geo -3 -o file.msh -v 99' [Gmsh 3.0.4, 1 node, max. 1 thread]
Info    : Started on Mon Aug  7 08:18:53 2017
Debug   : Destroying model 
Info    : Reading 'file.geo'...
Debug   : Extrude Surface 201
Debug   : Extrude Curve 1
Debug   : Extrude Point 1
Debug   : Extrude Point 2
Debug   : Extrude Curve 5
Debug   : Extrude Point 2
Debug   : Extrude Point 5
Debug   : Extrude Curve 6
Debug   : Extrude Point 5
Debug   : Extrude Point 4
Debug   : Extrude Curve 4
Debug   : Extrude Point 4
Debug   : Extrude Point 1
Debug   : Removed 40 duplicate points
Debug   : Removed 16 duplicate curves
Debug   : Extrude Surface 202
Debug   : Extrude Curve 2
Debug   : Extrude Point 2
Debug   : Extrude Point 3
Debug   : Extrude Curve 3
Debug   : Extrude Point 3
Debug   : Extrude Point 4
Debug   : Extrude Curve -6
Debug   : Extrude Point 4
Debug   : Extrude Point 5
Debug   : Extrude Curve -5
Debug   : Extrude Point 5
Debug   : Extrude Point 2
Debug   : Removed 43 duplicate points
Debug   : Removed 26 duplicate curves
Debug   : Removed 2 duplicate faces
Debug   : Extrude Surface 202
Debug   : Extrude Curve 2
Debug   : Extrude Point 2
Debug   : Extrude Point 3
Debug   : Extrude Curve 3
Debug   : Extrude Point 3
Debug   : Extrude Point 4
Debug   : Extrude Curve -6
Debug   : Extrude Point 4
Debug   : Extrude Point 5
Debug   : Extrude Curve -5
Debug   : Extrude Point 5
Debug   : Extrude Point 2
Debug   : Removed 40 duplicate points
Debug   : Removed 16 duplicate curves
Debug   : Syncing GEO_Internals with GModel
Debug   : Sync is removing 0 model entities
Debug   : GModel imported:
Debug   : 14 vertices
Debug   : 25 edges
Debug   : 15 faces
Debug   : 3 regions
Info    : Done reading 'file.geo'
Info    : Finalized high order topology of periodic connections
Info    : Meshing 1D...
Info    : Meshing curve 1 (Circle)
Info    : Meshing curve 2 (Circle)
Info    : Meshing curve 3 (Circle)
Info    : Meshing curve 4 (Circle)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 204 (extruded)
Info    : Meshing curve 205 (extruded)
Info    : Meshing curve 206 (extruded)
Info    : Meshing curve 207 (extruded)
Info    : Meshing curve 209 (extruded)
Info    : Meshing curve 210 (extruded)
Info    : Meshing curve 214 (extruded)
Info    : Meshing curve 218 (extruded)
Info    : Meshing curve 226 (extruded)
Info    : Meshing curve 227 (extruded)
Info    : Meshing curve 232 (extruded)
Info    : Meshing curve 248 (extruded)
Info    : Meshing curve 249 (extruded)
Info    : Meshing curve 250 (extruded)
Info    : Meshing curve 251 (extruded)
Info    : Meshing curve 253 (extruded)
Info    : Meshing curve 254 (extruded)
Info    : Meshing curve 258 (extruded)
Info    : Meshing curve 262 (extruded)
Info    : Done meshing 1D (0.00053 s)
Info    : Meshing 2D...
Info    : Meshing surface 201 (Plane, Delaunay)
Debug   : Computing edge loops
Debug   : Generating the mesh
Debug   : Meshing of the convex hull (8 points)
Debug   : Meshing of the convex hull (8 points) done
Debug   : Recovering 4 model Edges
Debug   : Recovering 8 mesh Edges (0 not recovered)
Debug   : Boundary Edges recovered for surface 201
Debug   : Computing mesh size field at mesh vertices 8
Debug   : Delaunizing the initial mesh
Debug   : Starting to add internal points
Debug   : Delaunization of the initial mesh done (0 swaps)
Debug   :       8 points created -- Worst tri radius is    0.900
Info    : Blossom: 11 internal 6 closed
Debug   : Perfect Match Starts 17 edges 10 nodes
Debug   : Perfect Match Succeeded in Quadrangulation (0 sec)
Debug   : 0 two-quadrangles vertices removed
Debug   : 0 diamond quads removed
Info    : Blossom completed (0.00 s):     4 quads     2 triangles 0 invalid quads  0 quads with Q < 0.1 Avg Q = 0.647 Min Q 0.590
Debug   : Type 17 2 triangles generated, 2 internal vertices
Info    : Meshing surface 202 (Plane, Delaunay)
Debug   : Computing edge loops
Debug   : Generating the mesh
Debug   : Meshing of the convex hull (8 points)
Debug   : Meshing of the convex hull (8 points) done
Debug   : Recovering 4 model Edges
Debug   : Recovering 8 mesh Edges (0 not recovered)
Debug   : Boundary Edges recovered for surface 202
Debug   : Computing mesh size field at mesh vertices 8
Debug   : Delaunizing the initial mesh
Debug   : Starting to add internal points
Debug   : Delaunization of the initial mesh done (0 swaps)
Debug   :       8 points created -- Worst tri radius is    0.900
Debug   : Type 17 10 triangles generated, 2 internal vertices
Info    : Meshing surface 211 (extruded)
Info    : Meshing surface 215 (extruded)
Info    : Meshing surface 219 (extruded)
Info    : Meshing surface 223 (extruded)
Info    : Meshing surface 224 (extruded)
Info    : Meshing surface 233 (extruded)
Info    : Meshing surface 237 (extruded)
Info    : Meshing surface 246 (extruded)
Info    : Meshing surface 255 (extruded)
Info    : Meshing surface 259 (extruded)
Info    : Meshing surface 263 (extruded)
Info    : Meshing surface 267 (extruded)
Info    : Meshing surface 268 (extruded)
Info    : Done meshing 2D (0.001076 s)
Info    : Meshing 3D...
Info    : Meshing volume 1 (extruded)
Info    : Meshing volume 2 (extruded)
Info    : Meshing volume 3 (extruded)
Info    : Subdividing extruded mesh
Info    : Swapping 0
Info    : Remeshing surface 255
Info    : Meshing surface 255 (extruded)
Info    : Remeshing surface 259
Info    : Meshing surface 259 (extruded)
Info    : Remeshing surface 263
Info    : Meshing surface 263 (extruded)
Info    : Remeshing surface 267
Info    : Meshing surface 267 (extruded)
Info    : Done meshing 3D (0.001202 s)
Info    : Optimizing 3D mesh...
Info    : Done optimizing 3D mesh (0 s)
Info    : Meshing order 2 (curvilinear on)...
Info    : Meshing curve 1 order 2
Info    : Meshing curve 2 order 2
Info    : Meshing curve 3 order 2
Info    : Meshing curve 4 order 2
Info    : Meshing curve 5 order 2
Info    : Meshing curve 6 order 2
Info    : Meshing curve 204 order 2
Info    : Meshing curve 205 order 2
Info    : Meshing curve 206 order 2
Info    : Meshing curve 207 order 2
Info    : Meshing curve 209 order 2
Info    : Meshing curve 210 order 2
Info    : Meshing curve 214 order 2
Info    : Meshing curve 218 order 2
Info    : Meshing curve 226 order 2
Info    : Meshing curve 227 order 2
Info    : Meshing curve 232 order 2
Info    : Meshing curve 248 order 2
Info    : Meshing curve 249 order 2
Info    : Meshing curve 250 order 2
Info    : Meshing curve 251 order 2
Info    : Meshing curve 253 order 2
Info    : Meshing curve 254 order 2
Info    : Meshing curve 258 order 2
Info    : Meshing curve 262 order 2
Info    : Meshing surface 201 order 2
Info    : Meshing surface 202 order 2
Info    : Meshing surface 211 order 2
Info    : Meshing surface 215 order 2
Info    : Meshing surface 219 order 2
Info    : Meshing surface 223 order 2
Info    : Meshing surface 224 order 2
Info    : Meshing surface 233 order 2
Info    : Meshing surface 237 order 2
Info    : Meshing surface 246 order 2
Info    : Meshing surface 255 order 2
Info    : Meshing surface 259 order 2
Info    : Meshing surface 263 order 2
Info    : Meshing surface 267 order 2
Info    : Meshing surface 268 order 2
Info    : Meshing volume 1 order 2
Info    : Meshing volume 2 order 2
Info    : Meshing volume 3 order 2
Info    : Finalized high order topology of periodic connections
Info    : Surface mesh: worst distortion = 1 (0 elements in ]0, 0.2]); worst gamma = 0.683282
Info    : Volume Mesh: worst distortion = 1 (0 elements in ]0, 0.2])
Info    : Done meshing order 2 (0.001448 s)
Info    : 163 vertices 171 elements
Info    : Writing 'file.msh'...
Info    : Done writing 'file.msh'
Info    : Stopped on Mon Aug  7 08:18:53 2017
Debug   : Destroying model file
-------------- next part --------------
A non-text attachment was scrubbed...
Name: .gmshrc
Type: application/octet-stream
Size: 3301 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170807/d63aa844/attachment-0003.obj>