Added the possibility to start from a random unitary for the Edmiston-Ruedenberg localization procedure