L'instruction suivante provoque une erreur lorsque le vecteur w_threshold est nul.
Cela ne devrait pas arriver, mais ça arrive !
hhi = int(np.round((np.linalg.norm(w_threshold, ord=1)/np.linalg.norm(w_threshold, ord=2))**2, decimals=0))
Suggestion: ajouter un epsilon à la norme L2 au dénominateur, quand celle-ci est nulle.