Function quat_mul_prim 

pub(crate) fn quat_mul_prim<B: Backend, const D: usize>(
    a: F<B, D>,
    b: F<B, D>,
) -> F<B, D>

Quaternion Hamilton product on the packed trailing (w, x, y, z) axis (primitive port of rotation::quat_mul).