A harmonic restarted Arnoldi algorithm for calculating eigenvalues and determining multiplicity