private static void ReflectSpheresRecursive(Sphere[] simplex, Vector3D[] poles, HashSet <Vector3D> completedPoles)
{
if (0 == poles.Length)
{
return;
}
HashSet <Vector3D> newPoles = new HashSet <Vector3D>(new AntipodeEqualityComparer(stereo: true));
foreach (Vector3D pole in poles)
{
for (int m = 0; m < simplex.Length; m++)
{
Sphere mirror = simplex[m];
Vector3D reflected = mirror.ReflectPoint(pole);
GreatSphere gs = new GreatSphere(reflected);
if (completedPoles.Add(gs.Pole))
{
newPoles.Add(gs.Pole);
}
}
}
ReflectSpheresRecursive(simplex, newPoles.ToArray(), completedPoles);
}
public void Gen(int p, int q, int r)
{
Geometry g = Util.GetGeometry(p, q, r);
if (g != Geometry.Spherical)
{
throw new System.Exception("Point group code only for spherical geometry.");
}
Simplex simplex = new Simplex();
simplex.Facets = SimplexCalcs.Mirrors(p, q, r);
Vector3D cen = new Vector3D();
Vector3D faceCenter = SimplexCalcs.FaceCenterSpherical(p, q, r);
Vector3D edgeMid = SimplexCalcs.EdgeMidpointSpherical(p, q, r);
Vector3D vertex = SimplexCalcs.VertexSpherical(p, q, r);
List <Vector3D> startingPoles = new List <Vector3D>();
startingPoles.Add(Sterographic.S3toR3(GreatSphere.FromSphere(simplex.Facets[0]).Pole));
GreatSphere[] spheres = CalcSpheres(simplex.Facets, startingPoles.ToArray());
string filename = "point_group.pov";
using (StreamWriter sw = File.CreateText(filename))
{
//foreach( var greatSphere in spheres )
//sw.WriteLine( PovRay.Sphere( greatSphere.ToSphere() ) );
}
/*
* List<H3.Cell.Edge> startingEdges = new List<H3.Cell.Edge>();
* startingEdges.Add( new H3.Cell.Edge( cen, faceCenter ) );
* H3.Cell.Edge[] edges = Recurse.CalcEdges( simplex.Facets, startingEdges.ToArray(), new Recurse.Settings() { G = Geometry.Spherical, Threshold = 0.001 } );
* //edges = edges.Where( e => !( Infinity.IsInfinite( e.Start ) || Infinity.IsInfinite( e.End ) ) ).ToArray();
*
* PovRay.WriteEdges( new PovRay.Parameters() { AngularThickness = 0.01 }, Geometry.Spherical, edges, filename, append: true );
*/
double minRad = 0;
double thick = 0.005;
System.Func <Vector3D, Sphere> sizeFunc = v =>
{
Vector3D c;
double rad;
H3Models.Ball.DupinCyclideSphere(v, thick / 2, g, out c, out rad);
return(new Sphere()
{
Center = c, Radius = Math.Max(rad, minRad)
});
};
// All geodesics
List <Circle3D> startingCircles = new List <Circle3D>();
startingCircles.Add(GeodesicFrom2Points(cen, edgeMid));
Circle3D[] geodesics = CalcGeodesics(simplex.Facets, startingCircles.ToArray());
Vector3D color = new Vector3D(0, 0, 1);
using (StreamWriter sw = File.CreateText(filename))
{
Shapeways shapeways = new Shapeways();
foreach (Circle3D geodesic in geodesics)
{
Vector3D[] points;
if (Infinity.IsInfinite(geodesic.Radius))
{
double cutoff = 15;
Segment seg = Segment.Line(geodesic.Normal * cutoff, geodesic.Normal * -cutoff);
points = seg.Subdivide(42);
}
else
{
List <Vector3D> tempPoints = geodesic.Subdivide(150).ToList();
tempPoints.Add(tempPoints[0]);
tempPoints.Add(tempPoints[1]);
points = tempPoints.ToArray();
}
List <Vector3D> ePoints = new List <Vector3D>();
List <double> eRadii = new List <double>();
foreach (Vector3D pNE in points)
{
Sphere sphere = sizeFunc(pNE);
ePoints.Add(sphere.Center);
eRadii.Add(sphere.Radius);
}
shapeways.AddCurve(ePoints.ToArray(), eRadii.ToArray());
sw.WriteLine(PovRay.EdgeSphereSweep(points, sizeFunc, color));
}
STL.SaveMeshToSTL(shapeways.Mesh, "533.stl");
}
}