private static CadFigure GetTriangleWithCenterPoint(VertexList pointList, int cpIndex)
{
int i1 = cpIndex;
int endi = pointList.Count - 1;
int i0 = i1 - 1;
int i2 = i1 + 1;
if (i0 < 0)
{
i0 = endi;
}
if (i2 > endi)
{
i2 = 0;
}
var triangle = CadFigure.Create(CadFigure.Types.POLY_LINES);
CadVertex tp0 = pointList[i0];
CadVertex tp1 = pointList[i1];
CadVertex tp2 = pointList[i2];
triangle.AddPoint(tp0);
triangle.AddPoint(tp1);
triangle.AddPoint(tp2);
triangle.IsLoop = true;
return(triangle);
}
public static List <CadFigure> Split(CadFigure fig, int curveSplitNum = 32)
{
CadVertex p0 = default(CadVertex);
var triangles = new List <CadFigure>();
int i1 = -1;
int state = 0;
CadFigure triangle;
VertexList pointList = fig.GetPoints(curveSplitNum);
i1 = CadUtil.FindMaxDistantPointIndex(p0, pointList);
if (i1 == -1)
{
return(triangles);
}
triangle = GetTriangleWithCenterPoint(pointList, i1);
Vector3d tp0 = triangle.PointList[0].vector;
Vector3d tp1 = triangle.PointList[1].vector;
Vector3d tp2 = triangle.PointList[2].vector;
Vector3d dir = CadMath.Normal(tp1, tp0, tp2);
Vector3d currentDir = Vector3d.Zero;
while (pointList.Count > 3)
{
if (state == 0)
{
i1 = CadUtil.FindMaxDistantPointIndex(p0, pointList);
if (i1 == -1)
{
return(triangles);
}
}
triangle = GetTriangleWithCenterPoint(pointList, i1);
tp0 = triangle.PointList[0].vector;
tp1 = triangle.PointList[1].vector;
tp2 = triangle.PointList[2].vector;
currentDir = CadMath.Normal(tp1, tp0, tp2);
bool hasIn = ListContainsPointInTriangle(pointList, triangle);
double scala = CadMath.InnerProduct(dir, currentDir);
if (!hasIn && (scala > 0))
{
triangles.Add(triangle);
pointList.RemoveAt(i1);
state = 0;
continue;
}
if (state == 0)
{
state = 1;
i1 = 0;
}
else if (state == 1)
{
i1++;
if (i1 >= pointList.Count)
{
break;
}
}
}
if (pointList.Count == 3)
{
triangle = CadFigure.Create(CadFigure.Types.POLY_LINES);
triangle.AddPoints(pointList, 0, 3);
triangle.IsLoop = true;
triangles.Add(triangle);
}
return(triangles);
}
