/// <summary> /// This is the method that actually does the work. /// </summary> /// <param name="DA">The DA object is used to retrieve from inputs and store in outputs.</param> protected override void SolveInstance(IGH_DataAccess DA) { MeshGhData hE_MeshData = new MeshGhData(); if (!DA.GetData(0, ref hE_MeshData)) { return; } Paramdigma.Core.HalfEdgeMesh.Mesh hE_Mesh = hE_MeshData.Value; if (!hE_Mesh.IsTriangularMesh()) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Mesh is not triangular!"); return; } List <Vector3d> normals = new List <Vector3d>(); List <Point3d> centroids = new List <Point3d>(); List <Point3d> circumcenters = new List <Point3d>(); foreach (Paramdigma.Core.HalfEdgeMesh.MeshFace face in hE_Mesh.Faces) { Paramdigma.Core.Geometry.Vector3d v = Paramdigma.Core.Geometry.MeshGeometry.FaceNormal(face); Paramdigma.Core.Geometry.Point3d centroid = Paramdigma.Core.Geometry.MeshGeometry.Centroid(face); normals.Add(new Vector3d(v.X, v.Y, v.Z)); centroids.Add(new Point3d(centroid.X, centroid.Y, centroid.Z)); } DA.SetDataList(0, normals); DA.SetDataList(1, centroids); }
/// <summary> /// This is the method that actually does the work. /// </summary> /// <param name="DA">The DA object can be used to retrieve data from input parameters and /// to store data in output parameters.</param> protected override void SolveInstance(IGH_DataAccess DA) { MeshGhData hE_MeshData = new MeshGhData(); if (!DA.GetData(0, ref hE_MeshData)) { return; } Paramdigma.Core.HalfEdgeMesh.Mesh hE_Mesh = hE_MeshData.Value; if (!hE_Mesh.IsTriangularMesh()) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Mesh is not triangular!"); return; } List <Vector3d> normals = new List <Vector3d>(); foreach (MeshVertex v in hE_Mesh.Vertices) { Paramdigma.Core.Geometry.Vector3d vect = new Paramdigma.Core.Geometry.Vector3d(); switch (Selection) { case NormalSelection.Equal: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalEquallyWeighted(v); break; case NormalSelection.Area: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalAreaWeighted(v); break; case NormalSelection.Angle: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalAngleWeighted(v); break; case NormalSelection.Gauss: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalGaussCurvature(v); break; case NormalSelection.Mean: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalMeanCurvature(v); break; case NormalSelection.Sphere: vect = Paramdigma.Core.Geometry.MeshGeometry.VertexNormalSphereInscribed(v); break; } normals.Add(new Vector3d(vect.X, vect.Y, vect.Z)); } DA.SetDataList(0, normals); }
/// <summary> /// This is the method that actually does the work. /// </summary> /// <param name="DA">The DA object is used to retrieve from inputs and store in outputs.</param> protected override void SolveInstance(IGH_DataAccess DA) { MeshGhData hE_MeshData = new MeshGhData(); if (!DA.GetData(0, ref hE_MeshData)) { return; } Paramdigma.Core.HalfEdgeMesh.Mesh hE_Mesh = hE_MeshData.Value; if (!hE_Mesh.IsTriangularMesh()) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Mesh is not triangular!"); return; } List <double> k1 = new List <double>(); List <double> k2 = new List <double>(); List <double> K = new List <double>(); List <double> km = new List <double>(); foreach (MeshVertex v in hE_Mesh.Vertices) { double[] k = Paramdigma.Core.Geometry.MeshGeometry.PrincipalCurvatures(v); k1.Add(k[0]); k2.Add(k[1]); K.Add(Paramdigma.Core.Geometry.MeshGeometry.ScalarGaussCurvature(v)); km.Add(Paramdigma.Core.Geometry.MeshGeometry.ScalarMeanCurvature(v)); } DA.SetDataList(0, k1); DA.SetDataList(1, k2); DA.SetDataList(2, K); DA.SetDataList(3, km); }
/// <summary> /// This is the method that actually does the work. /// </summary> /// <param name="DA">The DA object is used to retrieve from inputs and store in outputs.</param> protected override void SolveInstance(IGH_DataAccess DA) { MeshGhData hE_MeshData = new MeshGhData(); List <double> scalarValues = new List <double>(); List <double> levels = new List <double>(); string key = "sets1"; if (!DA.GetData(0, ref hE_MeshData)) { return; } if (!DA.GetDataList(1, scalarValues)) { return; } if (!DA.GetDataList(2, levels)) { return; } Paramdigma.Core.HalfEdgeMesh.Mesh hE_Mesh = hE_MeshData.Value; // Check for invalid inputs if (!hE_Mesh.IsTriangularMesh()) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Mesh is not triangular!"); return; } if (levels.Count == 0) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Level count can't be 0!"); return; } if (scalarValues.Count == 0) { AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Value count can't be 0!"); return; } // Assign values to mesh vertices foreach (MeshVertex v in hE_Mesh.Vertices) { if (v.UserValues.ContainsKey(key)) { v.UserValues[key] = scalarValues[v.Index]; } else { v.UserValues.Add(key, scalarValues[v.Index]); } } // Compute level sets Paramdigma.Core.Curves.LevelSets.ComputeLevels(key, levels, hE_Mesh, out List <List <Paramdigma.Core.Geometry.Line> > levelLines); //Convert AR_Lib.Curve.Line to Rhino.Geometry.Line DataTree <Line> resultLevel = new DataTree <Line>(); int count = 0; foreach (List <Paramdigma.Core.Geometry.Line> tempList in levelLines) { List <Line> tempResult = new List <Line>(); foreach (Paramdigma.Core.Geometry.Line l in tempList) { tempResult.Add(new Line(new Point3d(l.StartPoint.X, l.StartPoint.Y, l.StartPoint.Z), new Point3d(l.EndPoint.X, l.EndPoint.Y, l.EndPoint.Z))); } resultLevel.AddRange(tempResult, new GH.Kernel.Data.GH_Path(count)); count++; } // Return list of lines DA.SetDataTree(0, resultLevel); }
/// <summary> /// This is the method that actually does the work. /// </summary> /// <param name="DA">The DA object is used to retrieve from inputs and store in outputs.</param> protected override void SolveInstance(IGH_DataAccess DA) { MeshGhData hE_MeshData = new MeshGhData(); if (!DA.GetData(0, ref hE_MeshData)) { return; } Paramdigma.Core.HalfEdgeMesh.Mesh hE_Mesh = hE_MeshData.Value; List <Point3d> vertices = new List <Point3d>(); List <Line> edges = new List <Line>(); List <Rhino.Geometry.Mesh> faces = new List <Rhino.Geometry.Mesh>(); foreach (MeshVertex v in hE_Mesh.Vertices) { vertices.Add(new Point3d(v.X, v.Y, v.Z)); } foreach (MeshEdge e in hE_Mesh.Edges) { MeshVertex v1 = e.HalfEdge.Vertex; MeshVertex v2 = e.HalfEdge.Twin.Vertex; edges.Add(new Line(new Point3d(v1.X, v1.Y, v1.Z), new Point3d(v2.X, v2.Y, v2.Z))); } foreach (Paramdigma.Core.HalfEdgeMesh.MeshFace f in hE_Mesh.Faces) { List <MeshVertex> vs = f.AdjacentVertices(); List <int> faceVs = new List <int>(); List <Point3d> facePoints = new List <Point3d>(); int vi = 0; foreach (MeshVertex v in vs) { facePoints.Add(new Point3d(v.X, v.Y, v.Z)); faceVs.Add(vi); vi++; } Rhino.Geometry.Mesh m = new Rhino.Geometry.Mesh(); m.Vertices.AddVertices(facePoints); if (vs.Count == 3) { m.Faces.AddFace(0, 1, 2); } else if (vs.Count == 4) { m.Faces.AddFace(0, 1, 2, 3); } faces.Add(m); } DA.SetDataList(0, vertices); DA.SetDataList(1, edges); DA.SetDataList(2, faces); }