Structure
Mini Mandelbrot sets
Name:
- midget
- mini mandelbrot set
- baby Mandelbrots
- island = island mu-molecule = island mu-unit [1]
- primitive component
Parameter rays of mini Mandelbrot sets [2]
" a saddle-node (parabolic) periodic point in a complex dynamical system often admits homoclinic points, and in the case that these homoclinic points are nondegenerate, this is accompanied by the existence of infinitely many baby Mandelbrot sets converging to the saddle-node parameter value in the corresponding parameter plane." Devaney [3]
" Warped midgets in the Mandelbrot set have been measured, using an algorithm that allows the positions of the head, and cardioid atoms (north and south) of any midget to be found, once one has placed the cursor on the computer terminal somewhere inside any midget. We describe two distortions of midgets: linear distortions and angular distortions. When the north and south angles are plotted in the north/south angle plane, families of points are formed. The angle and distance measures of warped midgets from the Sea Horse Valley of the Mandelbrot set and from other sea horse valleys of midgets, whether on the Spike or on tendrils above atoms, all fall closely together in one part of the north/south plane. Measures of warped midgets from tendrils above the major atoms on the surface of the Cardioid fall closely together in another part of the north/south plane. This different way of looking at the Mandelbrot set offers an interesting way of studying the distortions of midgets." A. G. Davis Philip, Michael Frame, Adam Robucci: Warped midgets in the Mandelbrot set. Computers & Graphics 18(2): 239-248 (1994)
Distortion:
- Baby Mandelbrot Sets are Not Exact Self Similar. A baby Mandelbrot sets are homeomorphic to the original one
- fractalforums : how-distorted-can-a-minibrot-be/
- math.stackexchange question: mini-mandelbrots-are-they-exact-copies ?
- math.stackexchange question: why-does-the-mandelbrot-set-contain-slightly-deformed-copies-of-itself
Julia island = Embedded Julia Set = virtual Julia Sets
- "A structure comprised of filaments, resembling a Julia set in appearance, which has a higher delta Hausdorff dimension than filaments in the immediately surrounding region. They are sometimes also called Julia islands or virtual Julia Sets." Robert Munafo[4]
- "This location is called Julia island since it looks like a Julia set but it's actually inside Mandelbrot." [5]
order of hyperbolic componenets
- Myrberg 1963
- Sharkovsky 1964
- Metropolis 1973
- internal adress by Lau and Schleicher 1994
- rotation number by Devaney 1997
- R2 naming system by R Munafo 1999[6]
- orbit portrait by Milnor 2000
- shrubs by M Romero et al 2004 [7]
Topological models of Mandelbrot set
Speed improvements - optimisation
- general optimization
- Speed improvements by Robert Munafo
- fractalforums : help-optimising-mandelbrot
- The Computer Language Benchmarks Game
- Koebe 1/4 theorem by claudiusmaximus
Symmetry
The Mandelbrot set is symmetric with respect to the x-axis in the plane :
colour(x,y) = colour(x,-y)
its intersection with the x-axis ( real slice of Mandelbrot set ) is an interval :
<-2 ; 1/4>
It can be used to speed up computations ( up to 2-times )[13]
Bailout test
Instead of checking if magnitude ( radius = abs(z) ) is greater then escape radius ( ER):
compute square of ER:
ER2 = ER*ER
once and check :
It gives the same result and is faster.
interior detection
Interior detection method[14]
- time with detection versus without detection is : 0.383s versus 8.692s so it is 23 times faster !!!!
// gives last iterate = escape time
// output 0< i < iMax
int iterate(double complex C , int iMax)
{
int i=0;
double complex Z= C; // initial value for iteration Z0
complex double D = 1.0; // derivative with respect to z
for(i=0;i<iMax;i++)
{ if(cabs(Z)>EscapeRadius) break; // exterior
if(cabs(D)<eps) break; // interior
D = 2.0*D*Z;
Z=Z*Z+C; // complex quadratic polynomial
}
return i;
}
Period detection
"When rendering a Mandelbrot or Julia set, the most time-consuming parts of the image are the “black areas”. In these areas, the iterations never escape to infinity, so every pixel must be iterated to the maximum limit. Therefore, much time can be saved by using an algorithm to detect these areas in advance, so that they can be immediately coloured black, rather than rendering them in the normal way, pixel by pixel. FractalNet uses a recursive algorithm to split the image up into blocks, and tests each block to see whether it lies inside a “black area”. In this way, large areas of the image can be quickly coloured black, often saving a lot of rendering time. ... (some) blocks were detected as “black areas” and coloured black immediately, without having to be rendered. (Other) blocks were rendered in the normal way, pixel by pixel." Michael Hogg [15]
// cpp code by Geek3
// http://commons.wikimedia.org/wiki/File:Mandelbrot_set_rainbow_colors.png
bool outcircle(double center_x, double center_y, double r, double x, double y)
{ // checks if (x,y) is outside the circle around (center_x,center_y) with radius r
x -= center_x;
y -= center_y;
if (x * x + y * y > r * r)
return(true);
return(false);
// skip values we know they are inside
if ((outcircle(-0.11, 0.0, 0.63, x0, y0) || x0 > 0.1)
&& outcircle(-1.0, 0.0, 0.25, x0, y0)
&& outcircle(-0.125, 0.744, 0.092, x0, y0)
&& outcircle(-1.308, 0.0, 0.058, x0, y0)
&& outcircle(0.0, 0.25, 0.35, x0, y0))
{
// code for iteration
}
Cardioid and period-2 checking
One way to improve calculations is to find out beforehand whether the given point lies within the cardioid or in the period-2 bulb. Before passing the complex value through the escape time algorithm, first check if:
to determine if the point lies within the period-2 bulb and
to determine if the point lies inside the main cardioid.
another description [16]
// http://www.fractalforums.com/new-theories-and-research/quick-(non-iterative)-rejection-filter-for-mandelbrotbuddhabrot-with-benchmark/
public static void quickRejectionFilter(BlockingCollection<Complex> input, BlockingCollection<Complex> output)
{
foreach(var item in input.GetConsumingEnumerable())
{
if ((Complex.Abs(1.0 - Complex.Sqrt(Complex.One - (4 * item))) < 1.0)) continue;
if (((Complex.Abs(item - new Complex(-1, 0))) < 0.25)) continue;
if ((((item.Real + 1.309) * (item.Real + 1.309)) + item.Imaginary * item.Imaginary) < 0.00345) continue;
if ((((item.Real + 0.125) * (item.Real + 0.125)) + (item.Imaginary - 0.744) * (item.Imaginary - 0.744)) < 0.0088) continue;
if ((((item.Real + 0.125) * (item.Real + 0.125)) + (item.Imaginary + 0.744) * (item.Imaginary + 0.744)) < 0.0088) continue;
//We tried every known quick filter and didn't reject the item, adding it to next queue.
output.Add(item);
}
}
Periodicity checking
Most points inside the Mandelbrot set oscillate within a fixed orbit. There could be anything from ten to tens of thousands of points in between, but if an orbit ever reaches a point where it has been before then it will continually follow this path, will never tend towards infinity and hence is in the Mandelbrot set. This Mandelbrot processor includes periodicity checking (and p-2 bulb/cardioid checking) for a great speed up during deep zoom animations with a high maximum iteration value.
Perturbation theory
"The thing we call perturbation consist of 2 things: Calculate one pixel with high precision and use it as a reference for all other pixels.
This method will fail though, however thanks to Pauldelbrot we have a trustable method of detecting the pixels where the reference fails to compute the pixel with hardware precision. These pixels can be rendered with a reference closer to these pixels, so a typical perturbation render use several references. This method gives a speedup at about 10 times on depths of 10^100 Use a truncated Taylor series to approximate a starting value for a specific iteration, which make you able to skip all previous iterations. This method gives a speedup of typically Another 10 times on depths of 10^100, and together the speed up is typically 100 times. This, which we call Series Approximation, is where we have issues since we do not have any solid theoretically way of finding when too many iterations are skipped - for all pixels in the view. The more terms you include in the Taylor series, the more iterations you are able to skip. So if you stay below say 50 terms, it is not likely that you ever encounter any issues. Because some views can be correctly rendered 1000 or 100,000 times faster than full precision for each pixel with many terms - can you imagine a month turned into seconds! shocked K.I.Martin originally used only 3 terms " - Kalles Fraktaler [17]
Description
Very highly magnified images require more than the standard 64-128 or so bits of precision most hardware floating-point units provide, requiring renderers use slow "bignum" or "arbitrary precision"[18] math libraries to calculate. However, this can be sped up by the exploitation of perturbation theory.[19] Given
as the iteration, and a small epsilon, it is the case that
or
so if one defines
one can calculate a single point (e.g. the center of an image) using normal, high-precision arithmetic (z), giving a reference orbit, and then compute many points around it in terms of various initial offsets epsilon-zero plus the above iteration for epsilon. For most iterations, epsilon does not need more than 16 significant figures, and consequently hardware floating-point may be used to get a mostly accurate image.[20] There will often be some areas where the orbits of points diverge enough from the reference orbit that extra precision is needed on those points, or else additional local high-precision-calculated reference orbits are needed. This rendering method, and particularly the automated detection of the need for additional reference orbits and automated optimal selection of same, is an area of ongoing, active research. Renderers implementing the technique are publicly available and offer speedups for highly magnified images in the multiple orders of magnitude range.[21] [22] [23]
Newton-Raphson zooming
One can "use newton's method to find and progressively refine the precision of the location of the minibrot at the center of a pattern. This allows them to arbitrarily select a magnification between the location they started at and the final minibrot they calculate to be at the center of that location."[24][25]
Glitches
history
Programs
- SuperFractalThing in Java by K.I. Martin[33]
- mightymandel by Claude Heiland-Allen[34]
- mandelbrot-perturbator by Claude Heiland-Allen : http://code.mathr.co.uk/mandelbrot-perturbator/
- Kalles Fraktaler for windows
- perturbation_algebra and et
References
- ↑ the Mandelbrot Set Glossary and Encyclopedia, by Robert Munafo, (c) 1987-2018
- ↑ Parameter rays of mini mandelbrot sets
- ↑ Devaney In Global Analysis of Dynamical Systems, ed.: H. Broer, B. Krauskopf, G. Vegter. IOP Publishing (2001), 329-338 or Homoclinic Points in Complex Dynamical Systems
- ↑ embedded julia set from the Mandelbrot Set Glossary and Encyclopedia, by Robert Munafo, (c) 1987-2017.
- ↑ https://www.flickr.com/photos/nonnameavailable/28654921940/
- ↑ https://www.mrob.com/pub/muency/r2namingsystem.html
- ↑ M. Romera et al, Int. J. Bifurcation Chaos 13, 2279 (2003). https://doi.org/10.1142/S0218127403007941 Shrubs in the Mandelbrot Set Ordering
- ↑ Combinatorics in the Mandelbrot Set - Lavaurs Algorithm by
- ↑ Lavaurs algorithm by Michael Frame, Benoit Mandelbrot, and Nial Neger in lisp by Ruben Berenguel
- ↑ Abstract Mandelbrot tree by claudius maximus
- ↑ The Mandelbrot cactus by Robert L. Devaney
- ↑ Burns A M : : Plotting the escape- An animation of parabolic bifurrcation in the mandelbrot set. Mathematics Magazine: Volume 75, Number 2, Pages: 104-116 page 104
- ↑ How to use symetry of set
- ↑ A Cheritat wiki: Mandelbrot_set#Interior_detection_methods
- ↑ FractalNet by Michael Hogg
- ↑ Fractal forums: quick-(non-iterative)-rejection-filter-for-mandelbrotbuddhabrot-with-benchmark/
- ↑ Fractal Forums > Fractal Software > Help & Support > (C++) How to deep zoom in mandelbrot set?
- ↑ arbitrary precision at wikipedia
- ↑ perturbation theory in wikipedia
- ↑ "Superfractalthing - Arbitrary Precision Mandelbrot Set Rendering in Java by K.I. Martin 2013-05-18". http://www.fractalforums.com/announcements-and-news/superfractalthing-arbitrary-precision-mandelbrot-set-rendering-in-java/.
- ↑ "Kalles Fraktaler 2". http://www.chillheimer.de/kallesfraktaler/.
- ↑ Fast Mandelbrot set with infinite resolution, ver. 2 by Sergey Khashin, July 9, 2011
- ↑ Perturbation techniques applied to the Mandelbrot set by Claude Heiland-Allen October 21, 2013
- ↑ newton-raphson-zooming by quaz0r
- ↑ Newton-Raphson zooming and Evolution zoom method by Dinkydau
- ↑ pertubation-theory-glitches-improvement - fractal forum
- ↑ fractalNotes by Gerrit
- ↑ Fast Mandelbrot set with infinite resolution, ver. 2
- ↑ Fast calculation of the Mandelbrot set with infinite resolution by Sergey Khashin, October 12, 2016
- ↑ fractalforums : superfractalthing-arbitrary-precision-mandelbrot-set-rendering-in-java
- ↑ fractalforums - pertubation-theory-glitches-improvement
- ↑ Perturbation glitches by Claude Heiland-Allen
- ↑ superfractalthing By K.I. Martin
- ↑ mightymandel by Claude Heiland-Allen