Hi<div><br></div><div>I'm using the centreline field to mesh an artery. The stl for the artery and the centrelines along with the .geo file are attached (250 KB altogether).</div><div><br></div><div>Here's my my geometry file:</div>
<div>//////////////////////////////</div><div><div>Mesh.Algorithm = 6;</div><div>Mesh.Algorithm3D = 1;</div><div>Mesh.RemeshAlgorithm=1;</div><div>Mesh.RemeshParametrization = 1;</div><div><br></div><div>Merge "AbdominalAorta.stl";</div>
<div>Field[1] = Centerline;</div><div>Field[1].FileName = "VMTKCenterlinesOut.vtk";</div><div>Field[1].hLayer = 0.2;</div><div>Field[1].nbPoints = 20;</div><div>Field[1].closeVolume = 1;</div><div>Field[1].extrudeWall = 1;</div>
<div>Field[1].nbElemLayer = 4;</div><div>Field[1].hSecondLayer = 0;</div><div>Field[1].nbElemSecondLayer = 0;</div><div>Field[1].reMesh = 1;</div><div>Background Field = 1;</div><div>Field[1].run;</div></div><div>///////////////////////////////</div>
<div><br></div><div>I thought if give a negative value for hLayer then the extrusion will reverse direction, but it doesn't work. So I looked into the relevant source code and found the following in CenterlineField.cpp</div>
<div><div>  /////////////</div><div>  //orient extrude direction outward</div><div>  int dir = 0;</div><div>  MElement *e = current->getFaceByTag(1)->getMeshElement(0);</div><div>  SVector3 ne = e->getFace(0).normal();</div>
<div>  SVector3 ps(e->getVertex(0)->x(), e->getVertex(0)->y(), e->getVertex(0)->z());</div><div>  double xyz[3] = {ps.x(), ps.y(), ps.z()};</div><div>  kdtree->annkSearch(xyz, 1, index, dist);</div><div>
  SVector3 pc(nodes[index[0]][0], nodes[index[0]][1], nodes[index[0]][2]);</div><div>  SVector3 nc = ps-pc;</div><div>  if (dot(ne,nc) < 0) dir = 1;</div><div>  if (dir ==1 && hLayer > 0 ) hLayer *= -1.0;</div>
</div><div>/////////////</div><div><br></div><div>For testing purposes, I changed the value for dir to 1 (int dir = 1;) and commented out the conditional statement in the second last line and recompiled gmsh. Nothing happened, the extrusion is still oriented outwards. Is this a bug?</div>
<div><br></div><div>Thanks</div><div>Abdulhaq</div>