/// <summary>
///
/// </summary>
public static void Run()
{
// import halfedge mesh
var mesh = HeMesh3d.Factory.CreateFromOBJ(Paths.Resources + _fileIn);
// create particles
var particles = mesh.Vertices.Select(v => new Particle(v.Position)).ToArray();
// create constraints
var constraints = mesh.Faces.Where(f => !f.IsDegree3)
.Select(f => new PlanarNgon(f.Vertices.Select(v => v.Index)))
.ToArray();
// create solver
var solver = new ConstraintSolver();
// wait for keypress to start the solver
Console.WriteLine("Press return to start the solver.");
Console.ReadLine();
// step the solver until converged
while (!solver.IsConverged)
{
solver.StepParallel(particles, constraints);
Console.WriteLine($" step {solver.StepCount}");
}
Console.WriteLine("\nSolver converged! Press return to exit.");
// update mesh vertices
mesh.Vertices.Action(v => v.Position = particles[v.Index].Position);
// mesh.Vertices.Action(v => v.Position = particles[v].Position); // mesh elements (vertices, halfedges, faces) are converted to their indices implicitly so this also works
// compute vertex normals & write to file
mesh.GetVertexNormals(v => v.Position, (v, n) => v.Normal = n, true);
mesh.WriteToOBJ(Paths.Resources + _fileOut);
Console.ReadLine();
}
