  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.
