Logo by Dinkydau - Contribute your own Logo!

END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

it was a great time but no longer maintainable by c.Kleinhuis contact him for any data retrieval,
thanks and see you perhaps in 10 years again

this forum will stay online for reference
News: Did you know ? you can use LaTex inside Postings on fractalforums.com!
 
*
Welcome, Guest. Please login or register. March 28, 2024, 02:40:39 PM


Login with username, password and session length


The All New FractalForums is now in Public Beta Testing! Visit FractalForums.org and check it out!


Pages: [1]   Go Down
  Print  
Share this topic on DiggShare this topic on FacebookShare this topic on GoogleShare this topic on RedditShare this topic on StumbleUponShare this topic on Twitter
Author Topic: Bugs  (Read 2204 times)
0 Members and 1 Guest are viewing this topic.
Vega
Global Moderator
Conqueror
******
Posts: 143


From Russia With Love


« on: September 20, 2011, 03:46:10 PM »

Bugs



http://www.fractalforums.com/index.php?action=gallery;sa=view;id=8655

The Mandelbrot Set embarrass
Logged

cKleinhuis
Administrator
Fractal Senior
*******
Posts: 7044


formerly known as 'Trifox'


WWW
« Reply #1 on: September 20, 2011, 04:06:23 PM »

ehrm ... does this belong to our adult section !?

Logged

---

divide and conquer - iterate and rule - chaos is No random!
KRAFTWERK
Global Moderator
Fractal Senior
******
Posts: 1439


Virtual Surreality


WWW
« Reply #2 on: September 20, 2011, 04:38:51 PM »

Must be a bug!  afro
Logged

Pauldelbrot
Fractal Senior
******
Posts: 2592



pderbyshire2
« Reply #3 on: September 20, 2011, 07:29:05 PM »

I notice you have branch cuts in the spirals inside some of the buds (particularly the period-3 buds) inside the M-set.

I can get rid of those. Use this algorithm:

  • Iterate an interior, or suspected-interior, point w while keeping a saved earlier iteration t, initially equal to the starting point (for the normal quadratic M-set, this is the critical point 0). After every iteration compare the current iterate w with the saved one q. If they are extremely close, closer than some very small distance o, the point is periodic. At iterations 1, 2, 4, 8, 16, etc. update the saved iteration q to the current one w (after the comparison, of course). That's Brent's periodicity detection algorithm.
  • If the point escapes, do whatever. If it doesn't, but Brent doesn't fire, color it black.
  • Save the iterate w as t -- a point on the attractor. Iterate t forward until it's close to w again, counting these iterations, and save this count as p, the period of the attractor. Iterate t forward some more until the total number of iterations is congruent to zero mod p. Now t is very close to the specific attractor point that's in the same attracting basin component as the critical point we're studying.
  • Reiterate w from the starting point and check its distance from t whenever the iteration count i is congruent to zero mod p. Stop when w gets within a small disk (but much larger than o) around t. Save w, and save the iteration number i.
  • Iterate from w forward p iterations and save this point as y. Iterate from y forward p iterations and save this point as z. It should now be the case that w, y, and z are three successive iterates to land within the primary basin component, and that t is very close to the attractor point in the same component, in particular, very much closer than w is.
  • Compute a = (arg (t - y) - arg (t - w))/2π. This is the internal angle of the starting point.
  • Compute b = |t - w|, c = |t - y|, and d = |t - z| -- the distances of x, y, and z from the attractor.
  • Compute the ratio r = (b - c)/(a - b). This is the acceleration towards the attractor.
  • Compute log (((a - r)/(1 - r))(1/r - 1) + 1)/log(1/r). This is the fractional smoothed iteration.
  • Add i to that to get the smoothed iteration j.
  • The best spirals seem to result from setting colorindex = a + sf(j), where s determines the steepness of the spiral and f is a function. (It's assumed that colorindex 0 maps to one end of the gradient, 1 to the other, a gradient that starts and ends on the same color is chosen, and colorindex values out of the range [0,1] are taken as their fractional parts only.)
  • s = 10 and f(x) = 1/(px) (where p is the period from before) gives very similar results to this image, except that there is exactly one full repetition of the gradient around a circle in each bud or cardioid. In particular, if you used a green gradient that repeats four times you'd get a main cardioid very much like in this image, but there'd be four spiral arms in every bud, not two in the front bud and 4/3 in the side buds (which non-integral value may be the source of your branch cuts).
  • Note that all of the above can be done for any map of the complex plane, not just the traditional quadratic Mandelbrot.
Logged

Vega
Global Moderator
Conqueror
******
Posts: 143


From Russia With Love


« Reply #4 on: September 24, 2011, 06:58:57 PM »

I use other method. It is similar to yours. My method is intended for a 2D texture. Is your method intended for a 1D texture?


* F1051.jpg (53.83 KB, 500x500 - viewed 386 times.)
Logged

Pauldelbrot
Fractal Senior
******
Posts: 2592



pderbyshire2
« Reply #5 on: September 24, 2011, 08:41:33 PM »

For that, you want the internal angle of the point within the component. You seem to be currently getting the internal angle times the period, for whatever reason. My method for getting internal angles is to let an iterate get extremely close to a previous iterate, thus finding a point P of the attractor; then make a target disk, large in comparison to "extremely close", about this point and look at the first two iterates X, Y that lie inside that disk. The angle XPY is then a good approximation to the internal angle through the start point of the iteration, and can be computed as arg (X - P) - arg (Y - P) where arg is your complex number argument-determining function (likely based ultimately on atan, or atan2 if your programming language has that function) and subtraction of points is (equivalently) vector subtraction and complex number subtraction.

Note also that the number of iterations between X and Y gives the component's period.

The accuracy improves the smaller the target disk is compared to the component's size, and thus in general the pixel size, and also the smaller the "extremely close" distance is compared to the target disk size.
Logged

Vega
Global Moderator
Conqueror
******
Posts: 143


From Russia With Love


« Reply #6 on: September 28, 2011, 08:19:17 AM »

For that, you want the internal angle of the point within the component. You seem to be currently getting the internal angle times the period, for whatever reason. My method for getting internal angles is to let an iterate get extremely close to a previous iterate, thus finding a point P of the attractor; then make a target disk, large in comparison to "extremely close", about this point and look at the first two iterates X, Y that lie inside that disk. The angle XPY is then a good approximation to the internal angle through the start point of the iteration, and can be computed as arg (X - P) - arg (Y - P) where arg is your complex number argument-determining function (likely based ultimately on atan, or atan2 if your programming language has that function) and subtraction of points is (equivalently) vector subtraction and complex number subtraction.

Note also that the number of iterations between X and Y gives the component's period.

The accuracy improves the smaller the target disk is compared to the component's size, and thus in general the pixel size, and also the smaller the "extremely close" distance is compared to the target disk size.


I wait for a fraclal which you will draw by your technique.



* F1061.jpg (62.59 KB, 500x500 - viewed 341 times.)
Logged

Pauldelbrot
Fractal Senior
******
Posts: 2592



pderbyshire2
« Reply #7 on: December 14, 2011, 08:04:12 AM »

How's this grab you?



smiley
Logged

Pauldelbrot
Fractal Senior
******
Posts: 2592



pderbyshire2
« Reply #8 on: December 14, 2011, 09:47:23 AM »

Here's another variation on the theme, this time mapping the radial coordinate such that the images don't crowd toward the border of the M-set. It's very much like the one with the maple leaves or the one with the stars in terms of how the images are mapped ... save that the buds don't have fractional images, and indeed all rings of images in this rendering, in all buds and cardioids, have seven repetitions.

Logged

Vega
Global Moderator
Conqueror
******
Posts: 143


From Russia With Love


« Reply #9 on: December 14, 2011, 11:41:33 AM »

Pauldelbrot
Wow! Congratulate!


* 100.gif (41.53 KB, 100x100 - viewed 239 times.)
Logged

Pauldelbrot
Fractal Senior
******
Posts: 2592



pderbyshire2
« Reply #10 on: December 14, 2011, 08:39:28 PM »

Thanks!
Logged

Alef
Fractal Supremo
*****
Posts: 1174



WWW
« Reply #11 on: December 18, 2011, 06:00:33 PM »

Nice ladybirds. Could there be smaller images made of the large final images shocked
Logged

fractal catalisator
Pages: [1]   Go Down
  Print  
 
Jump to:  

Related Topics
Subject Started by Replies Views Last post
Bugs and Bugs II Images Showcase (Rate My Fractal) Pauldelbrot 0 743 Last post June 03, 2009, 06:50:58 AM
by Pauldelbrot
Bat Bugs Mandelbulb3D Gallery Weber 0 578 Last post March 17, 2012, 11:30:28 PM
by Weber
About bugs and glitches Movies Showcase (Rate My Movie) SeryZone 2 1129 Last post June 02, 2014, 11:00:15 PM
by SeryZone
Bugs Ultra Deep Mandelbrot CFJH 0 717 Last post May 16, 2016, 05:48:13 PM
by CFJH
v 2.12 alpha 6 bugs and missing features Bug Reporting joe 1 1250 Last post September 27, 2017, 06:17:59 AM
by mclarekin

Powered by MySQL Powered by PHP Powered by SMF 1.1.21 | SMF © 2015, Simple Machines

Valid XHTML 1.0! Valid CSS! Dilber MC Theme by HarzeM
Page created in 0.25 seconds with 24 queries. (Pretty URLs adds 0.011s, 2q)