/// <summary>
/// </summary>
/// <param name="threads"></param>
/// <param name="start"></param>
/// <param name="end"></param>
/// <param name="region"></param>
/// <param name="maxIterations"></param>
/// <param name="palette"></param>
/// <param name="gradient"></param>
/// <param name="bailout"></param>
/// <returns></returns>
public static byte[] DrawMandelbrot(Size threads, Size start, Size end, Region region, int maxIterations, RgbValue[] palette,
Gradient gradient, double bailout)
{
region = region.NormalizeRegion();
//var cartesianRegion = new Region(
// min: new Complex(region.Min.Real, -region.Min.Imaginary),
// max: new Complex(region.Max.Real, -region.Max.Imaginary));
int globalStartX = start.Width, globalEndX = end.Width;
int globalStartY = start.Height, globalEndY = end.Height;
int globalWidth = globalEndX - globalStartX,
globalHeight = globalEndY - globalStartY;
double globalRealStart = region.Min.Real,
globalImaginaryStart = region.Min.Imaginary;
var image = new byte[globalHeight * globalWidth * BitDepthFor24BppRgb];
int xThreads = threads.Width, yThreads = threads.Height;
var globalScale = MathUtilities.Scale(
globalWidth, globalHeight,
region.Min.Real, region.Max.Real,
region.Min.Imaginary, region.Max.Imaginary);
#region General Configuration For All Threads
if (maxIterations < 1)
{
throw new ArgumentException($"Max iterations must be >= 1, is {maxIterations}");
}
var colors = palette.Length;
var root = gradient.Exponent;
double indexScale = gradient.IndexScale, indexWeight = gradient.Weight;
var rootMinIterations = gradient.RootIndex ? Math.Pow(gradient.MinIterations, root) : 0.0;
var logBase = gradient.LogIndex ? Math.Log((double)maxIterations / gradient.MinIterations) : 0.0;
var logMinIterations = gradient.LogIndex ? Math.Log(gradient.MinIterations, logBase) : 0.0;
var logPaletteBailout = Math.Log(gradient.PaletteBailout);
var halfOverLogPaletteBailout = 0.5 / logPaletteBailout;
var bailoutSquared = bailout * bailout;
var useSqrt = gradient.RootIndex && Math.Abs(gradient.Root - 2) < Epsilon;
var maxIterationColor = gradient.MaxIterationColor;
#endregion
// TODO consider Parallel.For instead of threads
var tasks = new Thread[xThreads * yThreads];
for (var iy = 0; iy < yThreads; ++iy)
{
for (var ix = 0; ix < xThreads; ++ix)
{
var localRegion = MathUtilities.StartEndCoordinates(
globalStartX, globalEndX, globalStartY, globalEndY,
xThreads, ix, yThreads, iy);
int localStartX = localRegion.Item1,
localStartY = localRegion.Item3,
localEndX = localRegion.Item2,
localEndY = localRegion.Item4;
int localWidth = localEndX - localStartX,
localHeight = localEndY - localStartY;
var localStart =
MathUtilities.PixelToArgandCoordinates(
localStartX, localStartY,
globalScale.Real, globalScale.Imaginary,
globalRealStart, globalImaginaryStart);
var taskIndex = iy * xThreads + ix;
tasks[taskIndex] =
new Thread(() =>
DrawMandelbrotDC(
localStartX, localStartY,
localWidth, localHeight,
localStart.Real, localStart.Imaginary, maxIterations,
globalScale.Real, globalScale.Imaginary, globalWidth,
palette, gradient, image,
colors, bailoutSquared, halfOverLogPaletteBailout,
logBase, logMinIterations, root, rootMinIterations,
indexScale, indexWeight, useSqrt, maxIterationColor));
tasks[taskIndex].Start();
}
}
// Wait for completion
foreach (var task in tasks)
{
task.Join();
}
return(image);
}
/// <summary>
/// The following code was decompiled by JetBrains dotTrace.
/// I have subsequently made a couple of micro optimizations,
/// but it otherwise is equivalent to the DrawMandelbrot code
/// above.
/// </summary>
/// <param name="startX"></param>
/// <param name="startY"></param>
/// <param name="width"></param>
/// <param name="height"></param>
/// <param name="realStart"></param>
/// <param name="imaginaryStart"></param>
/// <param name="maxIterations"></param>
/// <param name="realScale"></param>
/// <param name="imaginaryScale"></param>
/// <param name="scan"></param>
/// <param name="palette"></param>
/// <param name="gradient"></param>
/// <param name="image"></param>
/// <param name="colors"></param>
/// <param name="bailoutSquared"></param>
/// <param name="halfOverLogBailout"></param>
/// <param name="logBase"></param>
/// <param name="logMinIterations"></param>
/// <param name="root"></param>
/// <param name="rootMinIterations"></param>
/// <param name="indexScale"></param>
/// <param name="indexWeight"></param>
/// <param name="useSqrt"></param>
/// <param name="maxIterationColor"></param>
private static void DrawMandelbrotDC(
int startX,
int startY,
int width,
int height,
double realStart,
double imaginaryStart,
int maxIterations,
double realScale,
double imaginaryScale,
int scan,
RgbValue[] palette,
Gradient gradient,
byte[] image,
int colors,
double bailoutSquared,
double halfOverLogBailout,
double logBase,
double logMinIterations,
double root,
double rootMinIterations,
double indexScale,
double indexWeight,
bool useSqrt,
Color maxIterationColor)
{
for (int index1 = 0; index1 < height; ++index1)
{
for (int index2 = 0; index2 < width; ++index2)
{
double num1 = 0.0;
double num2 = 0.0;
double num3 = 0.0;
double num4 = 0.0;
double d = 0.0;
double num5 = (double)index2 * realScale + realStart;
double num6 = (double)index1 * imaginaryScale + imaginaryStart;
int num7;
double num1num1 = 0.0;
double num3num3 = 0.0;
for (num7 = 0; d < bailoutSquared && num7 < maxIterations; ++num7)
{
double num8 = num1num1 - num3num3 + num5;
double num9 = 2.0 * num1 * num3 + num6;
double num10 = num8 - num1;
double num11 = num9 - num3;
if (num10 * num10 < 1E-20 && num11 * num11 < 1E-20)
{
num7 = maxIterations;
break;
}
double num12 = num8 - num2;
double num13 = num9 - num4;
if (num12 * num12 < 1E-20 && num13 * num13 < 1E-20)
{
num7 = maxIterations;
break;
}
num2 = num1;
num4 = num3;
num1 = num8;
num3 = num9;
num1num1 = num1 * num1;
num3num3 = num3 * num3;
d = num1num1 + num3num3; // save num1*num1, num3*num3, used above
}
double num14 = Math.Log(Math.Log(d) * halfOverLogBailout) * Mandelbrot.OneOverLog2;
double num15 = indexScale * ((double)(num7 + 1) - indexWeight * num14);
if (useSqrt)
{
num15 = Math.Sqrt(num15) - rootMinIterations;
}
else if (gradient.RootIndex)
{
num15 = Math.Pow(num15, root) - rootMinIterations;
}
if (gradient.LogIndex)
{
num15 = Math.Log(num15, logBase) - logMinIterations;
}
double index3 = MathUtilities.NormalizeIndex(num15, colors);
int index4 = MathUtilities.PreparePaletteIndex(index3, colors, gradient);
byte num16;
byte num17;
byte num18;
if (num7 >= maxIterations)
{
num16 = maxIterationColor.R;
num17 = maxIterationColor.G;
num18 = maxIterationColor.B;
}
else
{
RgbValue rgbValue = RgbValue.LerpColors(palette[index4], palette[(index4 + 1) % colors], index3 - (double)(long)index3);
num16 = (byte)rgbValue.red;
num17 = (byte)rgbValue.green;
num18 = (byte)rgbValue.blue;
}
int index5 = 3 * ((startY + index1) * scan + (startX + index2));
image[index5] = num18;
image[index5 + 1] = num17;
image[index5 + 2] = num16;
}
}
}
public static Image MakeImage(Region region, RgbValue[] palette, Gradient gradient)
{
return(MakeImage(region, palette, gradient, HIGH, WIDE));
}
/// <summary>
/// </summary>
/// <param name="threads"></param>
/// <param name="size"></param>
/// <param name="region"></param>
/// <param name="maxIteration"></param>
/// <param name="palette"></param>
/// <param name="gradient"></param>
/// <param name="bailout"></param>
/// <returns></returns>
public static byte[] DrawMandelbrot(Size threads, Size size, Region region, int maxIteration, RgbValue[] palette,
Gradient gradient, double bailout)
=> DrawMandelbrot(threads, new Size(0, 0), size, region, maxIteration, palette, gradient, bailout);
/// <summary>
/// </summary>
/// <param name="startX"></param>
/// <param name="startY"></param>
/// <param name="width"></param>
/// <param name="height"></param>
/// <param name="realStart"></param>
/// <param name="imaginaryStart"></param>
/// <param name="maxIterations"></param>
/// <param name="realScale"></param>
/// <param name="imaginaryScale"></param>
/// <param name="scan"></param>
/// <param name="palette"></param>
/// <param name="gradient"></param>
/// <param name="image"></param>
/// <param name="colors"></param>
/// <param name="bailoutSquared"></param>
/// <param name="halfOverLogBailout"></param>
/// <param name="logBase"></param>
/// <param name="logMinIterations"></param>
/// <param name="root"></param>
/// <param name="rootMinIterations"></param>
/// <param name="indexScale"></param>
/// <param name="indexWeight"></param>
/// <param name="useSqrt"></param>
/// <param name="maxIterationColor"></param>
private static void DrawMandelbrot(
int startX, int startY,
int width, int height,
double realStart, double imaginaryStart,
int maxIterations,
double realScale, double imaginaryScale,
int scan, RgbValue[] palette, Gradient gradient, byte[] image, int colors,
double bailoutSquared, double halfOverLogBailout,
double logBase, double logMinIterations,
double root, double rootMinIterations,
double indexScale, double indexWeight,
bool useSqrt, Color maxIterationColor)
{
for (var py = 0; py < height; ++py)
{
for (var px = 0; px < width; ++px)
{
double x = 0.0, xp = 0.0; // x;
double y = 0.0, yp = 0.0; // y;
var modulusSquared = 0.0; // x * x + y * y;
var x0 = px * realScale + realStart;
var y0 = py * imaginaryScale + imaginaryStart;
var iterations = 0;
while (modulusSquared < bailoutSquared && iterations < maxIterations)
{
var xtemp = x * x - y * y + x0;
var ytemp = 2 * x * y + y0;
double dx = xtemp - x,
dy = ytemp - y,
dxp = xtemp - xp,
dyp = ytemp - yp;
if ((dx * dx < Epsilon && dy * dy < Epsilon) ||
(dxp * dxp < Epsilon && dyp * dyp < Epsilon))
{
iterations = maxIterations;
break;
}
xp = x;
yp = y;
x = xtemp;
y = ytemp;
modulusSquared = x * x + y * y;
++iterations;
}
var smoothed = Math.Log(Math.Log(modulusSquared) * halfOverLogBailout) * OneOverLog2;
var index = indexScale * (iterations + 1 - indexWeight * smoothed);
if (useSqrt)
{
index = Math.Sqrt(index) - rootMinIterations;
}
else if (gradient.RootIndex)
{
index = Math.Pow(index, root) - rootMinIterations;
}
if (gradient.LogIndex)
{
index = Math.Log(index, logBase) - logMinIterations;
}
index = MathUtilities.NormalizeIndex(index, colors);
var actualIndex = MathUtilities.PreparePaletteIndex(index, colors, gradient);
byte red, green, blue;
if (iterations >= maxIterations)
{
red = maxIterationColor.R;
green = maxIterationColor.G;
blue = maxIterationColor.B;
}
else
{
var outval = RgbValue.LerpColors(
palette[actualIndex],
palette[(actualIndex + 1) % colors],
index - (long)index);
//Palette.Lerp(
// palette[actualIndex],
// palette[(actualIndex + 1) % colors],
// index - (long)index,
// out red, out green, out blue);
red = (byte)outval.red;
green = (byte)outval.green;
blue = (byte)outval.blue;
}
var offset = BitDepthFor24BppRgb * ((startY + py) * scan + (startX + px));
image[offset] = blue;
image[offset + 1] = green;
image[offset + 2] = red;
}
}
}
