/// <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)
{
int N1 = 0, N2 = 0, N3 = 0, N4 = 0, N5 = 0, N6 = 0;
DA.GetData(0, ref N1);
DA.GetData(1, ref N2);
DA.GetData(2, ref N3);
DA.GetData(3, ref N4);
DA.GetData(4, ref N5);
DA.GetData(5, ref N6);
IPrismElement e = new IPrismElement(N1, N2, N3, N4, N5, N6);
DA.SetData(0, e);
}
public static IMesh ExtrudeTwoDimensionalElements(IMesh mesh, IEnumerable <int> eKeys, double length)
{
IMesh dM = mesh.CleanCopy();
ITopologicVertex v, vv;
IVector3D n, pos;
IElement e, nE;
int next_vKey = mesh.FindNextVertexKey();
foreach (int eK in eKeys)
{
e = mesh.GetElementWithKey(eK);
if (e.TopologicDimension == 2)
{
dM.DeleteElement(eK, false);
nE = null;
mesh.Topology.ComputeTwoDimensionalElementNormal(eK, out n, out pos);
n *= length;
List <int> vertices = new List <int>(e.Vertices);
foreach (int vK in e.Vertices)
{
v = mesh.GetVertexWithKey(vK);
vv = new ITopologicVertex(v.Position + n, next_vKey);
dM.AddVertex(next_vKey, vv);
vertices.Add(next_vKey);
next_vKey++;
}
if (vertices.Count == 8)
{
nE = new IHexahedronElement(vertices.ToArray());
}
else if (vertices.Count == 6)
{
nE = new IPrismElement(vertices.ToArray());
}
if (nE != null)
{
dM.AddElement(nE);
}
}
}
dM.BuildTopology(true);
return(dM);
}
public static HashSet <int> TryParseToIguanaElement(int elementType, long[] nodes, int nodes_per_element, int number_of_elements, ref HashSet <int> parsedNodes, ref IMesh mesh)
{
if (IsElementImplemented(elementType))
{
for (int j = 0; j < number_of_elements; j++)
{
int[] eD = new int[nodes_per_element];
IElement e = null;
for (int k = 0; k < nodes_per_element; k++)
{
eD[k] = (int)nodes[j * nodes_per_element + k];
parsedNodes.Add(eD[k]);
}
switch (elementType)
{
//1st-order Triangle Face
case 2:
e = new ISurfaceElement(eD);
break;
//1st-order Quadrangle Face
case 3:
e = new ISurfaceElement(eD);
break;
//2nd-order 6-node triangle
case 9:
e = new ISurfaceElement.HighOrder.ITriangle6(eD);
break;
//2nd-order 9-node quadrangle
case 10:
e = new ISurfaceElement.HighOrder.IQuadrangle9(eD);
break;
//2nd-order 8-node quadrangle
case 16:
e = new ISurfaceElement.HighOrder.IQuadrangle8(eD);
break;
//3rd-order 9-node incomplete triangle
case 20:
e = new ISurfaceElement.HighOrder.ITriangle9(eD);
break;
//4th-order 12-node incomplete triangle
case 22:
e = new ISurfaceElement.HighOrder.ITriangle12(eD);
break;
//5th-order 15-node incomplete triangle
case 24:
e = new ISurfaceElement.HighOrder.ITriangle15(eD);
break;
//2nd-order 10-node tetrahedron
case 11:
e = new ITetrahedronElement.HighOrder.ITetrahedron10(eD);
break;
//1s-order 4-node tetrahedron element
case 4:
e = new ITetrahedronElement(eD);
break;
//2n-order 13-node pyramid
case 19:
e = new IPyramidElement.HighOrder.IPyramid13(eD);
break;
//1st-order 5-node pyramid element
case 7:
e = new IPyramidElement(eD);
break;
//1st-order 6-node prism element
case 6:
e = new IPrismElement(eD);
break;
//2nd-order 15-node prism
case 18:
e = new IPrismElement.HighOrder.IPrism15(eD);
break;
//1st-order 8-node hexahedron element
case 5:
e = new IHexahedronElement(eD);
break;
//2nd-order 20-node hexahedron
case 17:
e = new IHexahedronElement.HighOrder.IHexahedron20(eD);
break;
}
if (e != null)
{
mesh.AddElement(e);
}
}
}
return(parsedNodes);
}
