<div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi Thomas, <br></div><div><br></div><div>I first tried to set the random factors using: <br></div><div>Mesh.RandomFactor = 1e-9;<br>Mesh.RandomFactor3D = 1e-9;<br></div><div><br></div><div>I think the bug is elsewhere however...<br></div><div>With your file, I get a warning the first time and none the subsequent times.</div><div>I think it is due to internal state with the mesher remembering the choice of MeshAdapt?</div><div>Switching from Delauney to MeshAdapt halfway might cause different point spacing etc. <br></div><div><br></div><div>Here are some pictures showing the two consecutive mesh logs <br></div><div><br></div><div>1st time: <br></div><div><div><img src="cid:ii_jnqkbbgn0_166b26678a9a53e3" width="475" height="158"><br></div></div><div><br></div><div>2nd time (same window): <br></div><div><div><img src="cid:ii_jnqkd66w1_166b267cd5d8aa9c" width="407" height="179"><br></div></div><div><br></div><div>
On closing Gmsh and trying again, this repeats.</div><div><br></div><div>Instead, to fix your problem, you can set the meshing algo to MeshAdapt ahead of time by changing the line here: <br></div><div>Mesh.Algorithm = 1; // changed to MeshAdapt<br></div><div><br></div><div>Three repeated trials gave identical output files. <br></div><div><br></div><div>p.s. just in case, all of the current / modified options can be looked through and searched using Ctrl-Shift-H or Help -> Current Options and Workspace</div><div><br></div><div>Sincerely, <br></div><div>Max<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 26, 2018 at 12:44 PM, Thomas Kiel <span dir="ltr"><<a href="mailto:Thomas.Kiel@physik.hu-berlin.de" target="_blank">Thomas.Kiel@physik.hu-berlin.<wbr>de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Gmshers,<br>
<br>
as I already pointed out in an earlier email (s. <a href="http://onelab.info/pipermail/gmsh/2018/012045.html" rel="noreferrer" target="_blank">http://onelab.info/pipermail/g<wbr>msh/2018/012045.html</a>), I had problems with reproducing the exact same mesh twice. Since this email was unfortunately html encoded I think it was overlooked. Yet I like to apologize to ask again, this time in text formatted email, how to set the random number generator seed with an outside option?<br>
<br>
I know, that a grep to "srand" within the gmsh git repository gives several times hard-coded initialized seeds, e.g.<br>
<br>
Geo/SOrientedBoundingBox.cpp:3<wbr>26:  srand((unsigned)time(0));<br>
Mesh/Levy3D.cpp:624:  srand(time(NULL));<br>
Mesh/Levy3D.cpp:1359:    srand(time(NULL));<br>
Mesh/meshGFaceRecombine.cpp:24<wbr>5:  srand(time(NULL));<br>
Mesh/meshGFaceLloyd.cpp:394:  srand(time(NULL));<br>
<br>
yet I am not sure whether these influence my mesh generation. A minimalistic test file which causes different meshes at consecutive runs in my case, is<br>
<br>
<br>
//+<br>
SetFactory("OpenCASCADE");<br>
<br>
Mesh.CharacteristicLengthMin = 0.001;<br>
Mesh.CharacteristicLengthMax = 2.0;<br>
<br>
Mesh.CharacteristicLengthExten<wbr>dFromBoundary = 1;<br>
<br>
Mesh.Algorithm = 5;<br>
Mesh.Smoothing = 3;<br>
Mesh.Algorithm3D = 1;<br>
<br>
Box(1) = {0, -0.5, -0.5, 1, 1, 1};<br>
Box(2) = {-1, -0.5, -0.5, 1, 1, 1};<br>
Sphere(3) = {0.5, 0.0, 0.0, 0.3, -Pi/2, Pi/2, 2*Pi};<br>
<br>
Coherence;<br>
<br>
surfs_box1() = Unique(Abs(Boundary{ Volume{4}; } ));<br>
lines_box1() = Unique(Abs(Boundary{ Surface{surfs_box1()}; } ));<br>
pts_box1() = Unique(Abs(Boundary{ Line{lines_box1()}; } ));<br>
<br>
surfs_sphere() = Unique(Abs(Boundary{ Volume{3}; } ));<br>
lines_sphere() = Unique(Abs(Boundary{ Surface{surfs_sphere()}; } ));<br>
pts_sphere() = Unique(Abs(Boundary{ Line{lines_sphere()}; } ));<br>
<br>
surfs_box2() = Unique(Abs(Boundary{ Volume{2}; } ));<br>
lines_box2() = Unique(Abs(Boundary{ Surface{surfs_box2()}; } ));<br>
pts_box2() = Unique(Abs(Boundary{ Line{lines_box2()}; } ));<br>
<br>
// sorted from largest to smallest:<br>
Characteristic Length{pts_box2()} = 0.5;<br>
Characteristic Length{pts_sphere()} = 0.2;<br>
Characteristic Length{pts_box1()} = 0.1;<br>
<br>
<br>
<br>
I tested this on a machine running Ubuntu 16.04 and installed the out-of-the-box gmsh 4.0.4 version from the webpage. If any further detail on the running machine is required please tell me.<br>
I hope you can help me out.<br>
<br>
Thanks a lot in advance and best regards,<br>
<br>
Thomas Kiel<br>
<br>
-- <br>
Thomas Kiel<br>
Theoretische Optik und Photonik<br>
<br>
Humboldt-Universität zu Berlin, Institut für Physik<br>
<a href="https://maps.google.com/?q=Newtonstra%C3%9Fe+15,+12489+Berlin&entry=gmail&source=g" target="_blank">Newtonstraße 15, 12489 Berlin</a><br>
<br>
Tel.:   +49 30 2093 7999<br>
E-Mail: <a href="mailto:thomas.kiel@physik.hu-berlin.de" target="_blank">thomas.kiel@physik.hu-berlin.d<wbr>e</a><br>
Web:    <a href="https://top.physik.hu-berlin.de/people/thomas-kiel" rel="noreferrer" target="_blank">https://top.physik.hu-berlin.d<wbr>e/people/thomas-kiel</a><br>
<br>
______________________________<wbr>_________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/lis<wbr>tinfo/gmsh</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div class="m_7126871681863280953gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Max Orok</div><div>Contractor<br></div></div><div dir="ltr"><a href="http://www.mevex.com" target="_blank">www.mevex.com</a><br><br><img src="https://docs.google.com/uc?export=download&id=1fHTIiW4OMUjQr1iOkspQ7wiEsxunoOs0&revid=0B6x5w-5zVaEjSkpwbm5oY29jbG1XMzJoYldXTmJpNGFtb3dVPQ" width="164" height="42"><br></div></div></div></div></div></div>
</div></div>