A simultaneous coordinate relaxation algorithm for large, sparse matrix eigenvalue problems