Perhaps work with polar coordinates - points in the complement of the Mandelbrot set have external angle (I use Newton's method to trace rays, billtavis has an alternative using Runge-Kutta integration), and continuous iteration count heading towards infinity at the boundary of the set. Use the continuous iteration count as the hyperbolic distance from the origin. To convert hyperbolic distance to Euclidean distance in the Poincaré disc model, you need to invert the formula on the wiki page, that is, you have d(0,r) and want to find r, turns out
r = tanh(d/2)
(at least according to
http://www.wolframalpha.com/input/?i=solve+d+%3D+acosh%281+%2B+2+x^2+%2F+%281+-+x^2%29%29 which I didn't verify by hand, and I might have made a mistake in the input...).
But tanh(20) is already rounded to 1.0 in double precision, so this is only good to 40 iterations max (and will already be horribly grainy due to insufficient precision before then), and then double precision for external angles runs out at 53 iterations max. High precision maths is needed, no way out that I can see.
March 25, 2014, 01:22:16 AM
Dilber MC Theme by HarzeM