The solution of large undamped gyroscopic eigensystems by a subspace iteration method