Large-scale Exact Diagonalization of spin models with non-Abelian symmetries for the study of Many-Body Localization