internal static CoxeterImages.Settings AutoCalcScale(HoneycombDef def, Sphere[] mirrors)
{
// Calculate the color scale.
int size = 200;
CoxeterImages.Settings settings = new CoxeterImages.Settings()
{
Honeycomb = def,
Width = size,
Height = size,
Bounds = 1.1,
Mirrors = mirrors,
FileName = def.FormatFilename(),
};
CoxeterImages imageCalculator = new CoxeterImages();
imageCalculator.AutoCalcScale(settings);
if (settings.ColorScaling < 1)
{
settings.ColorScaling = 15;
}
return(settings);
}
public void AnimationSections(Settings config)
{
HoneycombDef imageData = new HoneycombDef(config.P, config.Q, config.R);
int p = imageData.P, q = imageData.Q, r = imageData.R;
string filename = imageData.FormatFilename();
Sphere[] mirrors = SimplexCalcs.Mirrors(p, q, r);
double bounds = 1.0; //config.UhsBoundary.Bounds;
bounds = 9.0;
// Calculate the color scale.
int size = 200;
CoxeterImages.Settings settings = new CoxeterImages.Settings()
{
Honeycomb = imageData,
Width = size,
Height = size,
Bounds = bounds,
Mirrors = mirrors,
FileName = imageData.FormatFilename(),
};
CoxeterImages imageCalculator = new CoxeterImages();
//imageCalculator.AutoCalcScale( settings );
if (settings.ColorScaling < 1)
{
settings.ColorScaling = 15;
}
settings.ColorScaling = 11;
Program.Log("\nGenerating sections...");
size = 500;
settings.Width = size;
settings.Height = size;
settings.FileName = filename;
double max = Spherical2D.e2sNorm(15);
double min = Spherical2D.e2sNorm(1.0 / 15);
DonHatch.e2hNorm(max);
int numSteps = 1800; // 1 minute
double step = (max - min) / numSteps;
for (int i = 0; i < 1; i++)
{
Program.Log("\nSection " + i);
imageCalculator.m_z = 1.0 / 0.5;
Spherical2D.s2eNorm(min + step * i);
DonHatch.h2eNorm(step * i);
settings.FileName = string.Format("533_{0:D4}.png", i);
imageCalculator.GenImage(settings);
}
}
public static void Create(HoneycombDef def, string filename)
{
int p = def.P;
int q = def.Q;
int r = def.R;
double scale = 5.0;
Vector3D cen = HoneycombPaper.InteriorPointBall;
Sphere[] simplex = SimplexCalcs.Mirrors(p, q, r, moveToBall: false);
// Apply transformations.
simplex = simplex.Select(s =>
{
Sphere.ScaleSphere(s, scale);
return(H3Models.UHSToBall(s));
}).ToArray();
for (int i = 0; i < 4; i++)
{
if (simplex[i].IsPointInside(cen))
{
simplex[i].Invert = true;
}
}
Sphere[] simplexForColorScale = SimplexCalcs.Mirrors(p, q, r, moveToBall: true);
CoxeterImages.Settings temp = HoneycombPaper.AutoCalcScale(def, simplexForColorScale);
int maxDepth = (int)temp.ColorScaling;
bool ball = true;
bool dual = false;
H3.Cell[] simplicesFinal = HoneycombPaper.GenCell(simplex, null, cen, ball, dual);
simplicesFinal = simplicesFinal.Where(s => s.Depths[0] < 1).ToArray();
//simplicesFinal = simplicesFinal.Where( s => s.)
// Output the facets.
using (StreamWriter sw = File.CreateText(filename)) // We need to reuse this StreamWriter (vs. calling AppendSimplex) for performance.
{
sw.WriteLine("#include \"hyper_ball.pov\"");
int[] include = new int[] { 0 };
foreach (H3.Cell cell in simplicesFinal)
{
Sphere[] facets = cell.Facets.Select(f => f.Sphere).ToArray();
int depth = cell.Depths[0] + 1;
Color c = Coloring.ColorAlongHexagon(maxDepth, depth);
PovRay.AddSimplex(sw, facets, cell.Center, include, filename, Coloring.ToVec(c));
}
}
}
public static void OneImage(Settings config, double t = 0.0)
{
HoneycombDef imageData = new HoneycombDef(config.P, config.Q, config.R);
int p = imageData.P, q = imageData.Q, r = imageData.R;
string filename = imageData.FormatFilename();
//if( File.Exists( filename ) )
// return;
Sphere[] mirrors = SimplexCalcs.Mirrors(p, q, r);
double bounds = config.UhsBoundary.Bounds;
// Calculate the color scale.
int size = 200;
CoxeterImages.Settings settings = new CoxeterImages.Settings()
{
Honeycomb = imageData,
Width = size,
Height = size,
Bounds = bounds,
Mirrors = mirrors,
FileName = imageData.FormatFilename(),
};
CoxeterImages imageCalculator = new CoxeterImages();
imageCalculator.AutoCalcScale(settings);
if (settings.ColorScaling < 1)
{
settings.ColorScaling = 15;
}
Program.Log("\nGenerating full image...");
size = 500;
settings.Width = config.UhsBoundary.ImageWidth;
settings.Height = config.UhsBoundary.ImageHeight;
settings.FileName = filename;
imageCalculator.GenImage(settings, t);
}
private static void OneAnimationFrame(HoneycombDef imageData, string filename, double colorScaling, double t = 0.0)
{
int p = imageData.P, q = imageData.Q, r = imageData.R;
Sphere[] mirrors = SimplexCalcs.Mirrors(p, q, r);
int size = 750;
size = 1024;
CoxeterImages.Settings settings = new CoxeterImages.Settings()
{
Honeycomb = imageData,
Width = size * 2,
Height = size,
Bounds = 1.0,
Mirrors = mirrors,
FileName = filename,
};
CoxeterImages imageCalculator = new CoxeterImages();
settings.ColorScaling = colorScaling;
imageCalculator.GenImage(settings, t);
}
private static void CreateCellPovRay(HoneycombDef def, string filename, double t = 0)
{
int p = def.P;
int q = def.Q;
int r = def.R;
//Vector3D trans = new Vector3D( 1.0/3, 0 ) * (2 + 2 * Math.Sin( Math.PI / 6 )) * t;
//double scale = 1.8;
Vector3D trans = new Vector3D();
double scale = 1.0;
Vector3D[] sVerts = null; // SimplexCalcs.VertsBall( p, q, r );
Vector3D vUHS = H3Models.BallToUHS(SimplexCalcs.VertexPointBall(p, q, r));
// Just did this for everything. Non-general position working better and will make all heads consistent.
scale = 2.0;
if (Geometry2D.GetGeometry(q, r) != Geometry.Hyperbolic) // Vertex-centered if possible
{
scale = 1.0 / vUHS.Z;
}
//else if( Geometry2D.GetGeometry( p, q ) == Geometry.Hyperbolic ) // Make the biggest head somewhat smaller.
// scale = 2.0;
Vector3D cen = InteriorPointBall;
/*var kleinVerts = sVerts.Select( v => HyperbolicModels.PoincareToKlein( v ) );
* Vector3D avg = new Vector3D();
* foreach( Vector3D v in kleinVerts )
* avg += v;
* avg /= kleinVerts.Count();
* Vector3D cen = HyperbolicModels.KleinToPoincare( avg );*/
cen = H3Models.BallToUHS(cen);
cen += trans;
//cen *= scale;
cen = H3Models.UHSToBall(cen);
Sphere[] simplex = SimplexCalcs.Mirrors(p, q, r, moveToBall: false);
// Apply transformations.
simplex = simplex.Select(s =>
{
Sphere.TranslateSphere(s, trans);
Sphere.ScaleSphere(s, scale);
return(H3Models.UHSToBall(s));
}).ToArray();
for (int i = 0; i < 4; i++)
{
if (simplex[i].IsPointInside(cen))
{
simplex[i].Invert = true;
}
}
Sphere[] simplexForColorScale = SimplexCalcs.Mirrors(p, q, r, moveToBall: true);
CoxeterImages.Settings temp = AutoCalcScale(def, simplexForColorScale);
int maxDepth = (int)temp.ColorScaling;
//Random rand = new Random( p+q+r );
//int randOffset = rand.Next( maxDepth );
bool ball = true;
bool dual = false;
H3.Cell[] simplicesFinal = GenCell(simplex, null, cen, ball, dual);
using (StreamWriter sw = File.CreateText(filename)) // We need to reuse this StreamWriter (vs. calling AppendSimplex) for performance.
{
sw.WriteLine("#include \"hyper_ball.pov\"");
//int[] include = new int[] { 0, 1, 2, 3 };
int[] include = new int[] { 0 };
if (dual)
{
include = new int[] { 3 }
}
;
// Output the facets.
foreach (H3.Cell cell in simplicesFinal)
{
Sphere[] facets = cell.Facets.Select(f => f.Sphere).ToArray();
if (m_toKlein)
{
facets = facets.Select(s => H3Models.BallToKlein(s)).ToArray();
}
int depth = cell.Depths[0] + 1;
Color c = Coloring.ColorAlongHexagon(maxDepth, depth);
if (cell.Depths.Sum() % 2 == 0)
{
c = Coloring.Inverse(c);
}
PovRay.AddSimplex(sw, facets, cell.Center, include, filename, Coloring.ToVec(c));
}
/*include = new int[] { 1, 2, 3 };
* foreach( H3.Cell cell in simplicesFinal )
* {
* Sphere[] facets = cell.Facets.Select( f => f.Sphere ).ToArray();
* Color c = Color.Red;
* Vector3D cv = Coloring.ToVec( c );
* cv.W = 0.9;
* PovRay.AddSimplex( sw, facets, cell.Center, include, filename, cv );
* }*/
}
// Output the edges/verts.
bool includeEdges = false;
if (includeEdges)
{
sVerts = sVerts.Select(v =>
{
v = H3Models.BallToUHS(v);
v += trans;
v *= scale;
return(H3Models.UHSToBall(v));
}).ToArray();
H3.Cell.Edge[] edges = Recurse.CalcEdges(simplex.Skip(1).ToArray(),
new H3.Cell.Edge[] { new H3.Cell.Edge(sVerts[2], sVerts[3], order: false) },
new Recurse.Settings()
{
Threshold = 0.01
});
PovRay.WriteH3Edges(new PovRay.Parameters {
AngularThickness = 0.01
}, edges, filename, append: true);
HashSet <Vector3D> verts = new HashSet <Vector3D>();
foreach (H3.Cell.Edge e in edges)
{
verts.Add(e.End);
}
PovRay.WriteVerts(new PovRay.Parameters {
AngularThickness = 0.02
}, Geometry.Hyperbolic, verts.ToArray(), filename, append: true);
}
}
