pub fn rotate_x<V: Float>(a: &[V; 4], angle: V) -> [V; 4]
Apply a rotation about the X-axis to this quaternion
Rotation about the X axis is c+i*s
(c+is) * a[] = ca - sa[i] + i.sa[r] - j.sa[k] + k.sa[j]