Discussion :

[ol9245]

Ma question est seulement pour l'esthétique.
Voici mono code à optimiser (lignes bleues). est-ce nécessaire de recopier la matrice d'origine ? la question n'est pas complètement triviale : dans mon application c'est une matrice creuse de 2600 x 2600 avec 260000 cases non vides.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
function [row, col] = find_special(M, row_indexes)
validateattributes(M, {'numeric'}, {'2d', 'nonnegative', 'nonempty'})
validateattributes(row_indexes, {'numeric'}, {'vector', 'integer', 'nonnegative', 'nonempty'})
% mettre à zéro les lignes ui ne sont pas ientifiées par row_indexes
M2 = sparse(size(M,1), size(M,2)) ;
M2(row_indexes, :) = M(row_indexes, :)

% faire le find sur ces lignes spécifiquement
[row, col] = find(M2) ;

end

exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
>> M = sparse(wilkinson(21)) ;
>> [row, col] = find_special(M, 7:10)
 
row =
 
     7
     7
     8
     7
     8
     9
     8
     9
    10
     9
    10
    10
 
 
col =
 
     6
     7
     7
     8
     8
     8
     9
     9
     9
    10
    10
    11
 
>>

[Jerome Briot]

Une alternative qui suppose que les valeurs dans row_indexes soient triées :

Code :

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

1
2
3
4
5
6
7
8
function [row, col] = find_special(M, row_indexes)
validateattributes(M, {'numeric'}, {'2d', 'nonnegative', 'nonempty'})
validateattributes(row_indexes, {'numeric'}, {'vector', 'integer', 'nonnegative', 'nonempty'})
 
[row, col] = find(M(row_indexes, :));
row(:) = row_indexes(row);
 
end

A tester.

[ol9245]

Très très fort !
merci
Ta présence sur le forum me donne vraiment envie de progresser et ça stimule l'activité de développement, qui est poutant assez barbante quand on ne sait pas la prendre du bon côté.