Logo by Trifox - Contribute your own Logo!


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: Visit the official fractalforums.com Youtube Channel
Welcome, Guest. Please login or register. January 29, 2023, 05:16:53 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
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: extremly high iteration buddhabrot  (Read 2214 times)
0 Members and 1 Guest are viewing this topic.
Fractal Molossus
Posts: 795

« on: March 10, 2010, 11:57:10 PM »

I'm still working on a stable, parallel, optimized version of my fortan buddhabrot generator.
while i work on it i'm playing with my PureBasic version.

i started explring buddhabroth with buddhabrot iteration > 1 million.
The problem is probably in the precision of the random number generator, but, for now, it seem that it's very hard to find plots that :
1) are not in the mset
2) do not escape before 1 million iteration.

i'm not sure about what i'm expecting to find yet ...

i currently run 2 version of the generator :
- one with 1 million maximum iteration but start plotting after 10 iteration
- one with 2 millions maximum iteration but start plotting after 1k iteration

the 1st one start to show some interesting result but i should wait for many more plot to really explore the result.
the 2nd run only show plot with less than 1 million iteration.

anyone tried something similar ? any tought about this experiment ?

often times... there are other approaches which are kinda crappy until you put them in the context of parallel machines
(en) http://www.blog-gpgpu.com/ , (fr) http://www.keru.org/ ,
Sysadmin & DBA @ http://www.over-blog.com/
Fractal Molossus
Posts: 795

« Reply #1 on: March 11, 2010, 12:42:26 AM »

i tought about something :

I first tought about a metropolis-hastings algorithm, but then tought about something else.

1) randomly explore the space
2) find a plot with high iteration
3) when a good candidate is found
3a) plot the orbit
3b) start a brownian motion with very tiny step
3c) plot the high iteration plot found
3d) stop the motion after <find a good stop condition>
4) goto 1

the stop condition could be a simple "do only 100 step"
something more complex like : stop the motion if it couldn't find another high iteration plot after 3 consecutive step

or something much more complex but probably too much cpu intensive like euristic, statistic, prediction, and so on... why not something inspired by http://en.wikipedia.org/wiki/Simulated_annealing (something i'd like to code sooner or later, sound fun)

often times... there are other approaches which are kinda crappy until you put them in the context of parallel machines
(en) http://www.blog-gpgpu.com/ , (fr) http://www.keru.org/ ,
Sysadmin & DBA @ http://www.over-blog.com/
Fractal Fertilizer
Posts: 362

The pwnge.

« Reply #2 on: March 11, 2010, 02:00:32 AM »

Rather than doing brownian, use directional estimation.

Hey, I just had an idea of how to actually do that! I'd been wondering about that for ages! YAAAY!

Anyway, you can derive a formula, similar to that for distance estimation, which will tell you the direction in which to head for higher iteration. You'll want to use a step size that's related to the current iteration value... or, more naturally, distance estimation. So something like Newz=OldZ+1/2*distEst*direst should lead you towards the Mset's deepest value at a (hopefully) exponential decrease in distance - here dirEst is a root of unity, btw. The only potential problem is that it might not "find" it's way into embedded Julia's very easily, due to their small gaps. Maybe doing some studies with the DE normals in that vicinity could lead to clues.

I recommend you try this. In order to do that direction estimation, you'll need to find the value of t such that d(t*DE(z))/dz is at it's minimum. Since DE is computed analytically using the chain rule, I think you should be able to do the same for direction estimation.

Wow, I wonder if anyone has ever thought of that before... they probably have.  crazy

Someday, man will understand primary theory; how every aspect of our universe has come about. Then we will describe all of physics, build a complete understanding of genetic engineering, catalog all planets, and find intelligent life. And then we'll just puzzle over fractals for eternity.
Fractal Senior
Posts: 1863

« Reply #3 on: March 11, 2010, 06:37:41 PM »

for those harder cases there could be a random treshold where, say, 10% of the points are tried out randomly, rather than with the direction estimation smiley
Pages: [1]   Go Down
Jump to:  

Related Topics
Subject Started by Replies Views Last post
First! (iteration) Philosophy oftakofta 1 2972 Last post December 19, 2009, 11:01:43 AM
by jehovajah
Low iteration anti-buddhabrot Mandelbrot & Julia Set ker2x 1 1558 Last post December 30, 2009, 03:51:32 AM
by ker2x
max iteration count as a function of zoom Mandelbrot & Julia Set makc 2 1391 Last post March 31, 2010, 10:08:48 PM
by Timeroot
Mandelbox smooth iteration count 3D Fractal Generation subblue 5 3570 Last post April 16, 2010, 05:14:25 PM
by knighty
High detail Buddhabrot - "Buddhabrot Echos" Images Showcase (Rate My Fractal) aluminumstudios 4 2262 Last post December 13, 2012, 04:42:29 PM
by tit_toinou

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.158 seconds with 24 queries. (Pretty URLs adds 0.006s, 2q)