[Gmsh] Lloyd smoothing broken?

[Nico Schlömer]

Hi everyone,

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
```
// Points
p1 = newp;
Point(p1) = {0, 0, 0, 0.1};
p2 = newp;
Point(p2) = {1, 0, 0, 0.1};
p3 = newp;
Point(p3) = {0, 1, 0, 0.1};
p4 = newp;
Point(p4) = {-1, 0, 0, 0.1};
p5 = newp;
Point(p5) = {0, -1, 0, 0.1};
// Circle arcs
c1 = newl;
Circle(c1) = {p2, p1, p3};
c2 = newl;
Circle(c2) = {p3, p1, p4};
c3 = newl;
Circle(c3) = {p4, p1, p5};
c4 = newl;
Circle(c4) = {p5, p1, p2};
l1 = newl;
Compound Line(l1) = {c1,c2,c3,c4};
ll1 = newll;
Line Loop(ll1) = {l1};
surf1 = news;
Plane Surface(surf1) = {ll1};
```
meshed with
```
gmsh -2 circle.geo -o circle.msh -optimize_lloyd 1000
```
I'm getting [1] when it should really look like [2].

Any hints on what might be going wrong?

Cheers,
Nico

[1] http://chunk.io/f/edd6327ec586485b9c8fe7a96ba1df9a
[2] http://chunk.io/f/5ba441e4700a4fad8be754c4f19dd9c6
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20170308/4a0c3ddf/attachment-0001.html>