public static ECMesh GetECMesh(Mesh m)
{
ECMesh ecm = new ECMesh();
for (int i = 0; i < m.Vertices.Count; i++)
{
ecm.AddVertex(m.Vertices[i]);
}
for (int i = 0; i < m.Faces.Count; i++)
{
ecm.AddFace(m.Faces[i].P0Index, m.Faces[i].P1Index, m.Faces[i].P2Index);
}
ecm.InitTypes();
return(ecm);
}
static void Main(string[] args)
{
Mesh m = Sample.GetSampleMesh2(100);
PlyManager.Output(m, "D://VTKproj//sample.ply");
ECMesh ecm = ECMesh.GetECMesh(m);
List <EdgeAndWeight> edgeWeights = new List <EdgeAndWeight>(ecm.Edges.Count);
Random r = new Random();
for (int i = 0; i < ecm.Edges.Count; i++)
{
EdgeAndWeight ew = new EdgeAndWeight();
ew.edgeID = i;
ew.weight = r.Next(0, 10000);
edgeWeights.Add(ew);
}
edgeWeights.Sort();
int deciCount = ecm.Edges.Count / 10;
int index = 0;
while (deciCount != 0)
{
if (index == edgeWeights.Count)
{
break;
}
EdgeAndWeight ew = edgeWeights[index];
ECEdge edge = ecm.Edges[ew.edgeID];
if (edge.IsValid() && edge.OrgV().Type == VertexType.INTERIOR && edge.DestV().Type == VertexType.INTERIOR)
{
ECVertex v0 = edge.OrgV();
ECVertex v1 = edge.DestV();
Point3d mid = new Point3d(v0.X / 2 + v1.X / 2, v0.Y / 2 + v1.Y / 2, v0.Z / 2 + v1.Z / 2);
ecm.Contract(v0, v1, mid);
deciCount--;
index++;
}
else
{
index++;
continue;
}
}
m = ecm.GetMesh();
PlyManager.Output(m, "test2.ply");
Console.WriteLine("\nCMP!");
Console.Read();
}
public QEMProcessor(Mesh m)
{
this.mesh = m;
ecm = ECMesh.GetECMesh(m);
heap = new RAHeap(ecm.Edges.Count);
}