Function quat_mul 

pub fn quat_mul(a: [f32; 4], b: [f32; 4]) -> [f32; 4]

Multiply two unit quaternions a and b (Hamilton product).

Format for both operands and the result: [x, y, z, w].