Function quat_mul 

pub fn quat_mul(q1: [f64; 4], q2: [f64; 4]) -> [f64; 4]

Multiply two quaternions q1 * q2 (Hamilton product).

Both inputs are [x, y, z, w].