/// <summary>
/// Creates a prism with the specified base.
/// </summary>
/// <param name="polygonBase2D">A 2D <see cref="GraphicsPath"/> representing the base of the prism.</param>
/// <param name="triangulationResolution">The resolution that will be used to linearise curve segments in the <see cref="GraphicsPath"/>.</param>
/// <param name="bottomOrigin">A <see cref="Point3D"/> that will correspond to the origin of the 2D reference system of the bottom base.</param>
/// <param name="topOrigin">A <see cref="Point3D"/> that will correspond to the origin of the 2D reference system of the top base.</param>
/// <param name="baseXAxis">A <see cref="NormalizedVector3D"/> that will correspond to the x axis of the 2D reference system of the bases. This will be orthonormalised to the <paramref name="baseYAxis"/>.</param>
/// <param name="baseYAxis">A <see cref="NormalizedVector3D"/> that will correspond to the y axis of the 2D reference system of the bases.</param>
/// <param name="fill">A collection of materials that will be applied to the <see cref="Triangle3DElement"/>s returned by this method.</param>
/// <param name="tag">A tag that will be applied to the <see cref="Triangle3DElement"/>s returned by this method.</param>
/// <param name="zIndex">A z-index that will be applied to the <see cref="Triangle3DElement"/>s returned by this method.</param>
/// <returns>A list of <see cref="Triangle3DElement"/>s that constitute the prism.</returns>
public static List <Element3D> CreatePrism(GraphicsPath polygonBase2D, double triangulationResolution, Point3D bottomOrigin, Point3D topOrigin, NormalizedVector3D baseXAxis, NormalizedVector3D baseYAxis, IEnumerable <IMaterial> fill, string tag = null, int zIndex = 0)
{
baseXAxis = (baseXAxis - baseYAxis * (baseXAxis * baseYAxis)).Normalize();
List <Element3D> tbr = new List <Element3D>();
bool orientation = (baseXAxis ^ baseYAxis) * (bottomOrigin - topOrigin) > 0;
double[,] matrix1 = Matrix3D.RotationToAlignAWithB(new NormalizedVector3D(0, 1, 0), baseYAxis);
double[,] matrix2 = Matrix3D.RotationToAlignAWithB(((Vector3D)(matrix1 * new Point3D(1, 0, 0))).Normalize(), baseXAxis);
List <List <NormalizedVector3D> > normals = (from el2 in polygonBase2D.GetLinearisationPointsNormals(triangulationResolution) select(from el in el2 select(el.X * baseXAxis + el.Y * baseYAxis).Normalize()).ToList()).ToList();
polygonBase2D = polygonBase2D.Linearise(triangulationResolution);
tbr.AddRange(CreatePolygon(polygonBase2D, triangulationResolution, bottomOrigin, baseXAxis, baseYAxis, orientation, fill, tag, zIndex));
tbr.AddRange(CreatePolygon(polygonBase2D, triangulationResolution, topOrigin, baseXAxis, baseYAxis, !orientation, fill, tag, zIndex));
List <List <Point3D> > bottomPoints = (from el2 in polygonBase2D.GetPoints() select(from el in el2 select(Point3D)(el.X * baseXAxis + el.Y * baseYAxis + (Vector3D)bottomOrigin)).ToList()).ToList();
List <List <Point3D> > topPoints = (from el2 in polygonBase2D.GetPoints() select(from el in el2 select(Point3D)(el.X * baseXAxis + el.Y * baseYAxis + (Vector3D)topOrigin)).ToList()).ToList();
if (orientation)
{
for (int i = 0; i < bottomPoints.Count; i++)
{
for (int j = 0; j < bottomPoints[i].Count - 1; j++)
{
tbr.AddRange(CreateRectangle(bottomPoints[i][j], bottomPoints[i][j + 1], topPoints[i][j + 1], topPoints[i][j], normals[i][j], normals[i][j + 1], normals[i][j + 1], normals[i][j], fill, tag, zIndex));
}
}
}
else
{
for (int i = 0; i < bottomPoints.Count; i++)
{
for (int j = 0; j < bottomPoints[i].Count - 1; j++)
{
tbr.AddRange(CreateRectangle(bottomPoints[i][j], topPoints[i][j], topPoints[i][j + 1], bottomPoints[i][j + 1], normals[i][j], normals[i][j], normals[i][j + 1], normals[i][j + 1], fill, tag, zIndex));
}
}
}
return(tbr);
}
