<div dir="ltr">Thanks Christophe for speedy explanation.<div><br></div><div>Cheers,</div><div>Nico</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Mar 8, 2017 at 8:17 PM Christophe Geuzaine <<a href="mailto:cgeuzaine@ulg.ac.be">cgeuzaine@ulg.ac.be</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br class="gmail_msg">
Hi Nico,<br class="gmail_msg">
<br class="gmail_msg">
It's not a bug, it's a feature (for quad mesh generation). See <a href="http://link.springer.com/chapter/10.1007%2F978-3-642-24734-7_26" rel="noreferrer" class="gmail_msg" target="_blank">http://link.springer.com/chapter/10.1007%2F978-3-642-24734-7_26</a> for some context.<br class="gmail_msg">
<br class="gmail_msg">
Christophe<br class="gmail_msg">
<br class="gmail_msg">
> On 8 Mar 2017, at 09:46, Nico Schlömer <<a href="mailto:nico.schloemer@gmail.com" class="gmail_msg" target="_blank">nico.schloemer@gmail.com</a>> wrote:<br class="gmail_msg">
><br class="gmail_msg">
> Hi everyone,<br class="gmail_msg">
><br class="gmail_msg">
> When using `-optimize_lloyd`, I noticed that gmsh doesn't produce centroidal tesselations but for some reason favors right angles. For the simple circle geometry<br class="gmail_msg">
> ```<br class="gmail_msg">
> // Points<br class="gmail_msg">
> p1 = newp;<br class="gmail_msg">
> Point(p1) = {0, 0, 0, 0.1};<br class="gmail_msg">
> p2 = newp;<br class="gmail_msg">
> Point(p2) = {1, 0, 0, 0.1};<br class="gmail_msg">
> p3 = newp;<br class="gmail_msg">
> Point(p3) = {0, 1, 0, 0.1};<br class="gmail_msg">
> p4 = newp;<br class="gmail_msg">
> Point(p4) = {-1, 0, 0, 0.1};<br class="gmail_msg">
> p5 = newp;<br class="gmail_msg">
> Point(p5) = {0, -1, 0, 0.1};<br class="gmail_msg">
> // Circle arcs<br class="gmail_msg">
> c1 = newl;<br class="gmail_msg">
> Circle(c1) = {p2, p1, p3};<br class="gmail_msg">
> c2 = newl;<br class="gmail_msg">
> Circle(c2) = {p3, p1, p4};<br class="gmail_msg">
> c3 = newl;<br class="gmail_msg">
> Circle(c3) = {p4, p1, p5};<br class="gmail_msg">
> c4 = newl;<br class="gmail_msg">
> Circle(c4) = {p5, p1, p2};<br class="gmail_msg">
> l1 = newl;<br class="gmail_msg">
> Compound Line(l1) = {c1,c2,c3,c4};<br class="gmail_msg">
> ll1 = newll;<br class="gmail_msg">
> Line Loop(ll1) = {l1};<br class="gmail_msg">
> surf1 = news;<br class="gmail_msg">
> Plane Surface(surf1) = {ll1};<br class="gmail_msg">
> ```<br class="gmail_msg">
> meshed with<br class="gmail_msg">
> ```<br class="gmail_msg">
> gmsh -2 circle.geo -o circle.msh -optimize_lloyd 1000<br class="gmail_msg">
> ```<br class="gmail_msg">
> I'm getting [1] when it should really look like [2].<br class="gmail_msg">
><br class="gmail_msg">
> Any hints on what might be going wrong?<br class="gmail_msg">
><br class="gmail_msg">
> Cheers,<br class="gmail_msg">
> Nico<br class="gmail_msg">
><br class="gmail_msg">
> [1] <a href="http://chunk.io/f/edd6327ec586485b9c8fe7a96ba1df9a" rel="noreferrer" class="gmail_msg" target="_blank">http://chunk.io/f/edd6327ec586485b9c8fe7a96ba1df9a</a><br class="gmail_msg">
> [2] <a href="http://chunk.io/f/5ba441e4700a4fad8be754c4f19dd9c6" rel="noreferrer" class="gmail_msg" target="_blank">http://chunk.io/f/5ba441e4700a4fad8be754c4f19dd9c6</a><br class="gmail_msg">
> _______________________________________________<br class="gmail_msg">
> gmsh mailing list<br class="gmail_msg">
> <a href="mailto:gmsh@onelab.info" class="gmail_msg" target="_blank">gmsh@onelab.info</a><br class="gmail_msg">
> <a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" class="gmail_msg" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a><br class="gmail_msg">
<br class="gmail_msg">
--<br class="gmail_msg">
Prof. Christophe Geuzaine<br class="gmail_msg">
University of Liege, Electrical Engineering and Computer Science<br class="gmail_msg">
<a href="http://www.montefiore.ulg.ac.be/~geuzaine" rel="noreferrer" class="gmail_msg" target="_blank">http://www.montefiore.ulg.ac.be/~geuzaine</a><br class="gmail_msg">
<br class="gmail_msg">
Free software: <a href="http://gmsh.info" rel="noreferrer" class="gmail_msg" target="_blank">http://gmsh.info</a> | <a href="http://getdp.info" rel="noreferrer" class="gmail_msg" target="_blank">http://getdp.info</a> | <a href="http://onelab.info" rel="noreferrer" class="gmail_msg" target="_blank">http://onelab.info</a><br class="gmail_msg">
<br class="gmail_msg">
</blockquote></div>