----- Original Message ----- From: JackOfTradeZ <JackOfTradeZ@comcast.net> To: fractint@mailman.xmission.com Subject: [Fractint] Re: quaternion product correction (Russell Walsmith) Date: Mon, 21 Nov 2005 20:01:06 -0700
I think it will look like a gyroscope in precess, and stuff will be flying off from it from centrifugal force. And the weird thing will be seen also.
So when I figure it all out, the final viddie will be a pirouette about the vertical axis, simultaneous with morph, and from a distance where it can all be seen.
I will alert / post it when done, may take a week or so.
JoTz
Weird thing? This gets curiouser and curiouser... looking forward to the viddie. There are some strange creatures in this section of 3-space, for sure... ---------------------------PAR FILE------------------------------- A_Mantabrot {;Gerald Dobiasovsky's 3D envelope wrapped around ; a T-set variant ; Version 2002 Patchlevel 5 reset=2002 type=formula formulafile=fromfile formulaname=metat-set passes=3 center-mag=-0.158929/0.144849/0.8743172/1/-44.9999999999999858/-\ 1.23373533611470521e-014 params=180/90/-2/2/0.5/0.5/0/4/150/253 float=y maxiter=2000000000 outside=summ periodicity=0 colors=000555<10>222222111<3>111111000000zzz<64>jjjjjjiii<5>hhhh\ hh000ggg<156>555 } ---------------------------FRM FILE------------------------------- frm:MetaT-set {;periodicity=no, outside=summ ;maxit > p5real*(p5imag+1) ;-------------------------------------------- ;p1real: Rotation about x-axis (1st rotation) ;p1imag: Rotation about y-axis (2nd rotation) ;p2real: Far clipping plane ;p2imag: Near clipping plane ;p3real: Pixel^2 coefficient ;p3imag: ;p4real: z1(0) ;p4imag: Bailout ;p5real: Maxiter per slice ;p5imag: Number of slices - 1 ;-------------------------------------------- ; bailout = imag(p4), tiefnum = imag(p5) delta = (real(p2)-imag(p2))/tiefnum tmp = pi/180 rotXax = exp(flip(real(p1)*tmp)), rotYax = exp(flip(imag(p1)*tmp)) ; HPixXY = rotYax VPixZ = real(rotXax) VPixXY = flip(conj(rotYax)) NXY = VPixZ*VPixXY NZ = imag(conj(rotXax)) VPixXY = -NZ*VPixXY ; tmp = NXY*imag(p2) + HPixXY*real(pixel) + VPixXY*imag(pixel) cx = cx0 = real(tmp), cy = cy0 = imag(tmp) cz = cz0 = NZ*imag(p2) + VPixZ*imag(pixel) ;HPixZ -> 0 tmp = NXY*delta, dcx = real(tmp), dcy = imag(tmp) dcz = NZ*delta x1 = 0 y1 = 0 z1 = real(p4) j = m = i = 0: a = sqr(x1) + 2*y1*z1 b = sqr(z1) + 2*x1*y1 c = sqr(y1) + 2*x1*z1 x1 = imag(p3)*a + real(p3)*a^2 + cx y1 = imag(p3)*b + real(p3)*b^2 + cy z1 = imag(p3)*c + real(p3)*c^2 + cz IF (bailout >= (sqr(x1)+sqr(y1)+sqr(z1))) i = i + 1 ELSE i = 0 m = m + 1 cx = cx0 = cx0 + dcx cy = cy0 = cy0 + dcy cz = cz0 = cz0 + dcz x1 = 0, y1 = 0, z1 = real(p4) ENDIF z = m - j j = j + 1 tiefnum >= m && p5 >= i } ------------------------------END PAR----------------------------