Discussion :

[ol9245]

Bonjour,

Tout est dans le titre !

j'ai une matrice creuse M (sparse matrix) de taille (300000 x 300000).
Les éléments non nuls de M sont tous positifs.
quelle est l'instruction Matlab qui me donne pour chaque ligne de M la valeur du plus petit élément non nul dans la ligne ?

Merci++

[ol9245]

Bon, c'était urgent.
J'ai trouvé un workaround.
j'ai fait ma matrice creuse avec les inverses des valeurs, j'ai pris le max, et j'ai reinversé en sortie.
je suis toujours preneur d'une solution plus rapide/élégante pour optimiser ma routine.

[Jerome Briot]

Il doit y avoir un moyen en utilisant SPFUN

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
M = [5 0 4 0 0
     6 0 0 0 1
     4 3 0 0 6
     0 0 0 0 0
     0 7 0 2 4]
 
sM = sparse(M);
 
X = spfun(@(x) 1./x,sM);
[a,b]=max(X,[],2);
idx = sub2ind(size(sM),1:size(sM,1),b.');
minsM = sM(idx)

ou :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
M = [5 0 4 0 0
     6 0 0 0 1
     4 3 0 0 6
     0 0 0 0 0
     0 7 0 2 4]
 
sM = sparse(M);
 
X = spfun(@(x) 1./x,sM);
X = spfun(@(x) 1./x,max(X,[],2));

[ol9245]

Il doit y avoir un moyen en utilisant SPFUN

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
M = [5 0 4 0 0
     6 0 0 0 1
     4 3 0 0 6
     0 0 0 0 0
     0 7 0 2 4]
 
sM = sparse(M);
 
X = spfun(@(x) 1./x,sM);
[a,b]=max(X,[],2);
idx = sub2ind(size(sM),1:size(sM,1),b.');
minsM = sM(idx)

ou :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
M = [5 0 4 0 0
     6 0 0 0 1
     4 3 0 0 6
     0 0 0 0 0
     0 7 0 2 4]
 
sM = sparse(M);
 
X = spfun(@(x) 1./x,sM);
X = spfun(@(x) 1./x,max(X,[],2));

Merci DUT. Toujours aussi pertinent !
je ne connaissais pas spfun.
mais il a une structure bizare qui lui donne des résultats inatendus sur certaines fonctions. Aparemment, il applique la fonction passée en argument sur le vecteur colonne des valeurs non nulles, puis affecte ce résultat à la matrice creuse de départ elle-même. D'ou des résultats totalement inatendus pour :

sparse_matrix = spfun(@(x) min(x,[],2), sparse_matrix) ; % rien de changé

et :

sparse_matrix = spfun(@min, sparse_matrix) ; % tous les éléments non nuls son remplacés par la plus petite valeur non nulle de la matrice complète