Riemannian Barzilai–Borwein method for the hermitian quaternion eigenvalue problem