Function quat_rotate_vec 

pub fn quat_rotate_vec(q: [f64; 4], v: [f64; 3]) -> [f64; 3]

Rotate a vector v by unit quaternion q.

Uses sandwich product: v’ = q * [v,0] * q*.