public Geometry <AltVertex> GenerateBorderGeometry <AltVertex>()
where AltVertex : struct, IVertex
{
if (borders == null || borders.Count == 0)
{
return(null);
}
else
{
List <AltVertex> vertices = new List <AltVertex>();
List <uint> newIndices = new List <uint>();
const float h = 0.1f;
foreach (var iter in borders)
{
Int64 borderKey = iter.Key;
int v1index = (int)(borderKey & 0xffffffff);
int v2index = (int)((borderKey >> 32) & 0xffffffff);
Border border = iter.Value;
// The border lies along an edge in the dual geometry.
int plateIndex1 = border.plate0;
int plateIndex2 = border.plate1;
Vector3 v1Pos = geometry.Mesh.GetPosition(v1index);
Vector3 v2Pos = geometry.Mesh.GetPosition(v2index);
Vector3 centroid1 = geometry.Topology.Centroids[border.c1Index].position;
Vector3 centroid2 = geometry.Topology.Centroids[border.c2Index].position;
Vector3 v1g1prime = Math2.BaseProjection(v1Pos, centroid1, centroid2, h);
Vector3 v1g2prime = Math2.BaseProjection(v1Pos, centroid2, centroid1, h);
Vector3 v2g1prime = Math2.BaseProjection(v2Pos, centroid1, centroid2, h);
Vector3 v2g2prime = Math2.BaseProjection(v2Pos, centroid2, centroid1, h);
centroid1 *= 1.01f; // Project the centroids out of the sphere slightly
centroid2 *= 1.01f;
bool edgeOrderIsAnticlockwise = false;
for (int i = 0; i < 3; i++)
{
if (geometry.Indices[border.c1Index * 3 + i] == v1index)
{
if (geometry.Indices[border.c1Index * 3 + (i + 1) % 3] == v2index)
{
edgeOrderIsAnticlockwise = true;
}
break;
}
}
int index = vertices.Count;
geometry.AddTriangle(ref vertices, ref newIndices, ref v1g2prime, ref centroid2, ref centroid1, edgeOrderIsAnticlockwise);
geometry.AddTriangle(ref vertices, ref newIndices, ref v1g2prime, ref centroid1, ref v1g1prime, edgeOrderIsAnticlockwise);
geometry.AddTriangle(ref vertices, ref newIndices, ref v2g1prime, ref centroid1, ref centroid2, edgeOrderIsAnticlockwise);
geometry.AddTriangle(ref vertices, ref newIndices, ref v2g1prime, ref centroid2, ref v2g2prime, edgeOrderIsAnticlockwise);
float p1 = Math2.Clamp(Math.Abs(borderCorners[border.c1Index].stress.pressure) * 1000.0f, 0, 1.0f);
float p2 = Math2.Clamp(Math.Abs(borderCorners[border.c2Index].stress.pressure) * 1000.0f, 0, 1.0f);
float hue1 = borderCorners[border.c1Index].stress.pressure > 0 ? 0 : 0.5f;
float hue2 = borderCorners[border.c2Index].stress.pressure > 0 ? 0 : 0.5f;
Vector3 rgb1 = Math2.HSV2RGB(new Vector3(hue1, p1, 1.0f - p1));
Vector3 rgb2 = Math2.HSV2RGB(new Vector3(hue2, p2, 1.0f - p2));
Vector4 c1Color = new Vector4(rgb1.X, rgb1.Y, rgb1.Z, 1);
Vector4 c2Color = new Vector4(rgb2.X, rgb2.Y, rgb2.Z, 1);
AltVertex v;
for (int i = 0; i < 12; ++i)
{
Vector4 color = new Vector4(1.0f, 1.0f, 1.0f, 1.0f);
if (i == 0 || i == 1 || i == 3 || i == 8 || i == 10 || i == 11)
{
color = c2Color;
}
else if (i == 2 || i == 4 || i == 5 || i == 6 || i == 7 || i == 9)
{
color = c1Color;
}
v = vertices[index + i]; MeshAttr.SetColor(ref v, ref color); vertices[index + i] = v;
}
}
Mesh <AltVertex> newMesh = new Mesh <AltVertex>(vertices.ToArray());
return(new Geometry <AltVertex>(newMesh, newIndices.ToArray()));
}
}
