# Task

How to show dynamics inside parabolic basin ?
Parameter c is a parabolic point
Attractor (limit set, attracting cycle / fixed point )
* belongs to Julia set ( not as in other cases )
* is weakly attracting ( lazy dynamics)

# Algorithm

Color depends on:
* sign of imaginary part of numerically approximated Fatou coordinate (numerical explanation )
* position of $z_{\infty }$ under $f^{p}$ in relation to one arm of critical orbit star: above or below ( geometric explanation)
* critical orbit is on the boundary between boxes and on the level curves of attraction time

## Name of the algorithm

* Parabolic chessboard = parabolic checkerboard
* zeros of qn - use only interior of Filled Julia sets on real slice of Mandelbrot set ( c is a real number, iaginary part is zero)
* Binary decomposition ( method = BDM ) in parabolic case
* Tessellation of the Interior of Filled Julia Sets ( Tomoki Kawahira in Semiconjugacies in Complex Dynamics with Parabolics)
* internal tiling

## Steps

First choose:
* choose target set, which is a circle with :
** center in parabolic fixed point
** radius such small that width of exterior between components is smaller then pixel width
* Target set consist of fragments of p components ( sectors). Divide each part of target set into 2 subsectors ( above and below critical orbit ) = binary decomposition
Steps:
* take $z_{0}$ = initial point of the orbit ( pixel)
* make forward iterations
** if points escapes = exterior
** if point do not escapes then check if point is near fixed point ( in the target set)
*** if no then make some extra iterations
*** if is then check in what half of target set it ( binary decomposition)

## How to find trap ?

Steps for parabolic basin
* choose component with critical point inside
* choose trap
* decompose trap disk for binary decomposition = divide into 2 parts : above and beyond the critical orbit
Trap features:
* is inside component with critical point inside
* trap has parabolic point on it's boundary
* center of the trap is midpoint between last point of critical orbit and fixed point
* radius of the trap is half of distance between fixed point and last point of critical orbit
Target set
* for periods 1 and 2 in case of complex quadratic polynomial : circle
* for periods >2: curvilinear sectors of circle around fixed point, 2 triangles described by :
Nota that one trap for components is bigger then two traps for BDM
* trap for componets : triangle $(z_{1},z_{f},z_{2})$
* trap for BDM : triangle $(f^{-p}(z_{cr}),z_{f},z_{cr})$
where
* z1, z2 are the last point of external rays bounding component ( usually the biggest one, with critical point inside )
* $z_{cr}$ is the critical point, or it's forward image in case thet critical orbit arm is not straight line
* $z_{f}$ is parabolic fixed point
Because of lazy dynamics one stops the computation before critical orbit and external rays land on the parabolic fixed point. The remaining gap is filled (approximated) with line from the last point to the landing point.
* computed and approximated parts

### How to compute preimages of critical point ?

Here a and b are 2 inverse ( or backward ) iterations ( multivalued with argument adjusted ) used in the program Mandel by Wolf Jung:
* 1-st inverse function = key a = $-f^{-1}(z)$
* 2-nd inverse function = key b = $+f^{-1}(z)$

# dictionary

* The chessboard is the name of this decomposition of A into a graph and boxes
* the chessboard graph
* the chessboard boxes :" The connected components of its complement in A are called the chessboard boxes (in an actual chessboard they are called squares but here they have infinitely many corners and not just four). "
* the two principal or main chessboard boxes
* trap = target set = attracting petal
Visualizing Structures with the Chessboard Graph

# description

"A nice way to visualize the extended Fatou coordinates is to make use of the parabolic graph and chessboard."
Color points according to :
* the integer part of Fatou coordinate
* the sign of imaginary part
Corners of the chessboard ( where four tiles meet ) are precritical points
$\bigcup _{n=0}^{n\geq 0}f^{-n*p}(z_{cr})$
or
$\{z:f^{n*p}(z)=z_{cr}\}$

## 1/1
# Images

Click on the images to see the code and descriptions on the Commons !
* basins
* BDM only in the interior
* BDM in both exterior ( superattracting) and interior ( parabolic)
* only one basin: exterior ( superattracting)
* on the real slice of Mandelbrot set
* on the cusp of main cardioid
* from main cardioid to period 2
* between period 2 and 4 along 1/2 ray
* on the boundary of main cardioid of Mandelbrot set
* 0/1 = 1/1
* 1/2
* 1/3
* 1/4
* 5/11
evolution on the escape route 0
* start with nucleus of component with period 1 (c = 0)
* go along internal ray for angle = 0
* parabolic point (c = 1/4)
* on the external ray for angle 0
* evolution on the internal ray of the main cardioid of Mandelbrot set
* superattracting
* attracting
* parabolic
* only one basin: exterior ( superattracting), Julia set discontinous
* other polynomials
* various number of star branches
* orbits inside sepals
* 0/1 = 1/1
* invariant curves inside main boxes
* DLD
* only boundaries of parabolic BDM
* d=2
* d=3
* d=4
* d=5
* structural chessboard of first renormalization of the degree 2 parabolic Blaschke product, cylinders highted with a normal map.
Examples :
* Tiles: Tessellation of the Interior of Filled Julia Sets by T Kawahira
* coloured cauliflower by A Cheritat

# code

For the internal angle 0/1 and 1/2 critical orbit is on the real line ( Im(z) = 0). It is easy to compute parabolic chessboard because one have to check only imaginary part of z. For other cases it is not so easy

## 0/1

Cpp code by Wolf Jung see function parabolic from file mndlbrot.cpp ( program mandel )
To see effect :
* run Mandel
* (on parameter plane ) find parabolic point for angle 0, which is c=0.25. To do it use key c, in window input 0 and return.
C code :
Gnuplot code :

## 1/2 or fat basilica

Cpp code by Wolf Jung see function parabolic from file mndlbrot.cpp ( program mandel ) To see effect :
* run Mandel
* (on parameter plane ) find parabolic point for angle 1/2, which is c=-0.75. To do it use key c, in window input 0 and return.
C code :
* trap for components

## 1/3

* chessboard
* target set ( precritical points)

## 1/4

* BDM for both exterior and interior

## 5/11

* extarnal rays and critical orbit (11-arm star)

# See also

* parabolic perturbation
* Checkerboard in Hyperbolic tilings by User:Tamfang : images and Python code
* https://plus.google.com/110803890168343196795/posts/Eun6pZVkkmA
* shadertoy: Orbit trapped julia Created by maeln in 2016-Jan-19
* Holomorphic checkerboard by etale_cohomology
* Sepals of cauliflower
* wikipedia :Zebra striping in computer graphics
* aproximation cauliflower by the true tree: tree which converge to the cauliflower julia set. the set of vertices is in one to one correspondence with the grand orbit of the critical point of z squared plus one quarter.

# references

* ↑ Near parabolic renormalization for unisingular holomorphic maps by Arnaud Cheritat
* ↑ About Inou and Shishikura’s near parabolic renormalization by Arnaud Cheritat
* ↑ Applications of near-parabolic renormalization by Mitsuhiro Shishikura
* ↑ Complex Dynamical Systems by Robert L. Devaney, page
* ↑ Antiholomorphic Dynamics: Topology of Parameter Spaces and Discontinuity of Straightening by Sabyasachi Mukherjee
* ↑ tiles by T Kawahira
* ↑ checkerboards by A Cheritat
* ↑ commons:Category:Fractals created with Mandel
* ↑ Program Mandel by Wolf Jung
* ↑ Program Mandel by Wolf Jung
* ↑ Oleg Ivrii (Tel Aviv University), "Shapes of trees"