Bonjour à tous,
Toujours dans l'exploration de Numpy...
Je trouve pas de fonction pour trouver le rang d'une matrice... Pour l'instant je fais :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
 
A = eye(3)
tolerance = 1e-16
s,v,d = linalg.svd(A)
rang = sum( where(v>tolerance, 1, 0))
Ce qui revient à compter le nombre d'éléments non nuls sur la diagonale de V issu de la SVD.
N'y-a-t il pas plus simple, et surtout moins consommateur de puissance de calcul?
Merci...