Numerical diagonalization of fermion mass matrices
  The diagonalization of general mass matrices is a more delicate problem when eigenvalue degeneracies exist. In this case, often overlooked in the literature, some difficulties arise related to the freedom in the choice of basis in degenerate subspaces. Here two simple algorithms are developed to deal with quark and neutrino mass matrices with arbitrary degeneracies. 
J. A. Aguilar-Saavedra