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

Dears I need some help with quaternonions. Mandelbulb3D uses quaternonion math that I do NOT know ...
What is the formula for;

(a + bi + cj + dk) * (e + fi + gj + hk) = fre() + i * fc1(...) + j * fc2(...) + k * fc3(...)

and

(a + bi + cj + dk) / (e + fi + gj + hk) = fre() + i * fc1(...) + j * fc2(...) + k * fc3(...)

I only am sure that (a + bi + cj + dk) + (e + fi + gj + hk) = (a+e) + (b+f) i + ... (Simple term by term sum)

I need float functions of float coefficients to be inserted in an assembly code.

 

I am not able to find formulas by myself. It is useful for some MB3D formulas; the second is the expression of 4D Mobius map (with a = (a1-a2) and stuff).

I also need;

(a + bi + cj ) * (d + ei + fj ) = fre() + i * fc1(...) + j * fc2(...)

and

(a + bi + cj ) / (d + ei + fj ) = fre() + i * fc1(...) + j * fc2(...)

(Quaternonions without fourth term for simplified code in 3D only)  

Darkbeam,

Your formula for the product (a₁ + b₁i + c₁j + d₁k)(a₂ + b₂i + c₂j + d₂k) looks right.

I believe the reciprocal of a quaternion is its conjugate divided by its squared norm. In other words, if z = a + bi + cj + dk, then

     z^-1 = z^* / ||z||² = (a - bi - cj - dk) / (a² + b² + c² + d²)

Here's a picky point about the quotient of two quaternions: It's not uniquely defined. We could interpret p/q as meaning either pq^-1 or q^-1p, and they're not generally equal due to the noncommutativity of quaternion multiplication.

The Wikipedia article on quaternions has loads of good information.

Regards,
Ken