Logo by kr0mat1k - 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: Support us via Flattr FLATTR Link
 
*
Welcome, Guest. Please login or register. March 28, 2024, 12:44:42 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: A different 3D Mandelbrot  (Read 4287 times)
0 Members and 1 Guest are viewing this topic.
pseudogenius
Guest
« on: May 09, 2010, 07:49:21 PM »

Hey,

I came up with this formula using my own set of rules for hypercomplex numbers.
It is in chaospro format.

parameter real bailout;
real a,b,c,d,j,k,l,m,costheta,sintheta,pa,pb,pc,pd;
   void init(void)
   {
      a=part_r(pixel);
      b=part_i(pixel);
      c=part_j(pixel);
      d=part_k(pixel);
      pa=a;
      pb=b;
      pc=c;
      pd=d;
      costheta=(a^2-b^2)/(a^2+b^2);
      sintheta=(2*a*b)/(a^2+b^2);
   }
   void loop(void)
   {
      j=a^2-b^2-costheta*(c^2)+2*sintheta*c*d+costheta*(d^2);
      k=2*a*b-sintheta*(c^2)-2*costheta*c*d+sintheta*(d^2);
      l=2*a*c-2*b*d;
      m=2*a*d+2*b*c;
      a=j;
      b=k;
      c=l;
      d=m;
      j=a^2-b^2-costheta*(c^2)+2*sintheta*c*d+costheta*(d^2);
      k=2*a*b-sintheta*(c^2)-2*costheta*c*d+sintheta*(d^2);
      l=2*a*c-2*b*d;
      m=2*a*d+2*b*c;
      a=j;
      b=k;
      c=l;
      d=m;
      j=a^2-b^2-costheta*(c^2)+2*sintheta*c*d+costheta*(d^2)+pa;
      k=2*a*b-sintheta*(c^2)-2*costheta*c*d+sintheta*(d^2)+pb;
      l=2*a*c-2*b*d+pc;
      m=2*a*d+2*b*c+pd;
      a=j;
      b=k;
      c=l;
      d=m;
      costheta=(a^2-b^2)/(a^2+b^2);
      sintheta=(2*a*b)/(a^2+b^2);
   }
   bool bailout(void)
   {
      return(j^2+k^2+l^2+m^2<bailout);
   }

It is power 8, so that is why the formula repeats in the loop.

So far I see infinite detail, but choaspro isn't letting me do a lot(it takes forever to render).

I have some renders of it.
The second image is what you get when you remove

      costheta=(a^2-b^2)/(a^2+b^2);
      sintheta=(2*a*b)/(a^2+b^2);

from the end of the loop.

Tell me what you think.


* JComplexP8Large.jpg (102.74 KB, 1280x800 - viewed 616 times.)

* JComplexP8AngleRemoved.jpg (41.21 KB, 915x675 - viewed 610 times.)
Logged
reesej2
Guest
« Reply #1 on: May 10, 2010, 01:10:30 AM »

Interesting gridlike pattern on it.

I think I understand your equations, but what was the rationale behind them? It doesn't really look like something you'd just make up on the spot.
Logged
pseudogenius
Guest
« Reply #2 on: May 13, 2010, 02:21:09 AM »

Well, I used the following rules:

<br />{j_{\theta}}^2=-cos (2\theta)-sin(2\theta)i<br />ij_{\theta}=j_{\theta+\frac{\pi}{2}}<br />j_{\theta+\pi}=-j_{\theta}<br />

They're based on the idea that you can multiply higher order complex numbers by adding the angles of each complex number and multiplying the lengths.

I got the expression in my formula from (a+bi+cj_{\theta}+dj_{\theta+\frac{\pi}{2}})^2

Expanding, and then taking each part I get the above formula with some trigonometric functions, mainly cos(2*theta) and sin(2*theta)

Recognizing that theta=atan(b/a) and using trig identities I get

cos(2*theta)=(a^2-b^2)/(a^2+b^2)
sin(2*theta)=(2*a*b)/(a^2+b^2)

which are the costheta and sintheta terms you see in the formula
Logged
Pages: [1]   Go Down
  Print  
 
Jump to:  

Related Topics
Subject Started by Replies Views Last post
Stripe Average Mandelbrot 3 Images Showcase (Rate My Fractal) Zephitmaal 0 1511 Last post January 06, 2012, 09:52:58 AM
by Zephitmaal
Mandelbrot In Stripes2 Images Showcase (Rate My Fractal) Ilionstar 0 998 Last post February 06, 2012, 01:33:44 AM
by Ilionstar
Having fun with scilab and the mandelbrot set Programming ker2x 1 5615 Last post February 22, 2012, 10:09:32 PM
by ker2x
Rhythmic Gymnastic Mandelbrot FractInt Gallery simon.snake 0 1167 Last post May 06, 2012, 06:46:51 PM
by simon.snake
Mandelbrot in SQL Programming jhaig 2 2452 Last post June 14, 2012, 02:23:54 PM
by jhaig

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