Logo by Fiery - 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: Visit us on facebook
 
*
Welcome, Guest. Please login or register. January 09, 2026, 02:08: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
  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: 3D distance estimates without using a derivative or Delta DE  (Read 5957 times)
0 Members and 1 Guest are viewing this topic.
fractalrebel
Fractal Lover
**
Posts: 211



WWW
« on: February 12, 2011, 10:06:24 PM »

I have a distance estimate method that doesn't use a delta to estimate the derivative. In fact, it doesn't use a derivative of any kind. Here are the equations:

pot = exp(-siter*log(2)
grad = exp(-(siter-iter)*log(2))
distance = dfactor*sinh(pot)/(2*exp(pot)*grad)

siter is a smoothed iteration. It can be obtained from a number of approaches.

If the Vepstas smoothing method is used
siter = iter + log(log(sqrt(|zri|+|zjk|)))/log(2)

The potential (pot) comes from some of Vepstas' treatments
The gradient (grad) is normally where the derivative would appear, and is approximated with a special derivative of the potential (also fromm Vepstas' writings). The distance equation comes from the book Hypercomplex Iterations.

Here is a redo of an image I posted here a while ago using the the distance estimate method above:



* MandelFlowers.jpg (115.98 KB, 640x640 - viewed 478 times.)
Logged

David Makin
Global Moderator
Fractal Senior
******
Posts: 2286



Makin' Magic Fractals
WWW
« Reply #1 on: February 13, 2011, 01:38:16 PM »

pot = exp(-siter*log(2)
grad = exp(-(siter-iter)*log(2))
distance = dfactor*sinh(pot)/(2*exp(pot)*grad)

Hi Ron, I think I'm correct in that the formula you quoted gives DE for a divergence of degree 2, I know how to change "siter" for other values but what, if anything, needs changing for other divergence values in the above (and how) ?
Also does this work generically e.g. for the Mandelbox and KIFS and all the hybrids ?
Logged

The meaning and purpose of life is to give life purpose and meaning.

http://www.fractalgallery.co.uk/
"Makin' Magic Music" on Jango
fractalrebel
Fractal Lover
**
Posts: 211



WWW
« Reply #2 on: February 13, 2011, 06:19:30 PM »

Hi Dave,

It works for all types and is especially good for mandelbox, conformal mandelbrot and the hybrids.

Other forms for siter include:

more general
Vepstas/Harkonen:
                pow = log(|zri - oldzri| + |zjk-oldzjk|)/log(|oldzri - oldzri2| + |oldzjk-oldzjk2|)
                siter = i + (log(log(bailout^2)) - log(log(1/(|zri - oldzri| + |zjk-oldzjk|))))/log(pow)

most general
General smoothing:
                modhold = sqrt(|oldzri-zri|+|oldzjk-zjk|)
                modh = sqrt(|zri|+|zjk|)
                siter = siter + exp(-modh-0.5/modhold)

The Vepstas/Harkonen method is applied after bailout. General Smoothing must be aplied each iteration cycle.

Logged

fractalrebel
Fractal Lover
**
Posts: 211



WWW
« Reply #3 on: February 13, 2011, 06:26:59 PM »

I forgot to mention that for mandelbox and the hybrids I used Buddhi's method to calculate what he calls distance, and then use that value in place of sqrt(|zri|+|zjk|)
« Last Edit: February 13, 2011, 07:19:25 PM by fractalrebel » Logged

fractalrebel
Fractal Lover
**
Posts: 211



WWW
« Reply #4 on: February 13, 2011, 07:56:31 PM »

Dave,

I must have missed it somewhere along the way, and search of this site didn't provide much. What is KIFS?
Logged

Buddhi
Fractal Iambus
***
Posts: 895



WWW
« Reply #5 on: February 13, 2011, 08:08:16 PM »

Kaleidoscopic Iterated Function System - http://www.fractalforums.com/3d-fractal-generation/kaleidoscopic-(escape-time-ifs)/
Logged

fractalrebel
Fractal Lover
**
Posts: 211



WWW
« Reply #6 on: February 13, 2011, 11:38:18 PM »

I was hospitalized for most of the time from March through June of 2010, so I missed a lot. In looking at the KIFS  code I see some statements like
rotate1(x,y,z)
which really confuses me, as x, y and z are regular variables and there is no mention of angles. Can anyone help me here? I would like to try my hand at some of Knightly's code and I can't without knowing what the functions are.
Please, someone help.
Logged

fractalrebel
Fractal Lover
**
Posts: 211



WWW
« Reply #7 on: February 14, 2011, 12:19:19 AM »

I found the rotate code  embarrass
Logged

Pages: [1]   Go Down
  Print  
 
Jump to:  

Related Topics
Subject Started by Replies Views Last post
Delta Images Showcase (Rate My Fractal) Kali 6 1820 Last post June 27, 2012, 06:28:37 AM
by Kali
Improved delta-DE algorithm 3D Fractal Generation Buddhi 9 4099 Last post August 16, 2013, 08:25:40 PM
by Roquen
Hybri-Station on Delta 4 Images Showcase (Rate My Fractal) JoeFRAQ 1 2691 Last post October 31, 2014, 02:09:53 AM
by JohnVV
Orange Delta Images Showcase (Rate My Fractal) C.K. 4 1550 Last post January 11, 2015, 06:46:21 PM
by C.K.
Running derivative General Discussion lapinot 3 2683 Last post December 14, 2016, 07:54:53 PM
by DarkBeam

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.277 seconds with 26 queries. (Pretty URLs adds 0.011s, 2q)