Discussion :

[ol9245]

Bonjour à tous,

J'ai la tête un peu dasn le pâté, donc je fais appel à vos lumières en cas de solution simple et triviale :
Je voudrais supprimer la boucle dans le code suivant. Note importante : le vecteur index possède des répétitions, comme dans le bout de code d'exemple.

Merci++ de votre aide.

Code :

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

1
2
3
4
5
6
index = [1 2 3 2 5 2 6] ; % c'est un exemple. dans la réalité index fait jusqu'à 10000 valeurs.
M = zeros(1, 10) ; % c'est un exemple. dans la réalité M est une matrice 1000 x 1000
a = 0.2 ;
for k=1:numel(index)
   M(index(k)) = M(index(k)) + a ;
end

[ol9245]

en attendant vos lumières, j'ai écrit cette hideuse béquille :

Code :

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

1
2
3
4
5
6
7
8
9
function M = ajouter (M, idx, a)
% ajoute à M la valeur a pour tous les index idx (avec répétition)
idx = sort(idx(:)) ; % trie les index à traiter
while numel(idx)>0
    df = [1; diff(idx)] ; % repère les doublons
    unefois = idx(df>0) ; % sélectionne les index sans doublons
    M(unefois) = M(unefois)+a ; % ajoute a à M sur les index non doublonnés
    idx = idx(df==0) ; % supprime les éléments traités
end

qui marche comme ça par exemple :

Code :

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

1
2
3
4
5
6
>> M = zeros(1,10) ;
M = ajouter (M, [5 4 3 2 1 5 1 2 3 2 3 4 4 5 5 6], 1)
 
M =
 
     2     3     3     3     4     1     0     0     0     0

peut-on faire plus élégant ? et plus rapide (éviter le tri, en particulier) ?
merci

[Jean Dumoncel]

Salut,

pour éviter la boucle :

Code :

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

1
2
3
4
5
6
7
index = [1 3 2 2 5 2 6] ;
M = zeros(1, 10) ;
a = 0.2 ;
 
[B,I,J] = unique(index);
h = hist(J,length(B));
MM(B) = M(B)+a.*h;

Par contre je vois pas comment tu pourrais éviter le tri (ici c'est la fonction unique qui s'en charge).

[ol9245]

Salut,

pour éviter la boucle :

Code :

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

1
2
3
4
5
6
7
index = [1 3 2 2 5 2 6] ;
M = zeros(1, 10) ;
a = 0.2 ;
 
[B,I,J] = unique(index);
h = hist(J,length(B));
MM(B) = M(B)+a.*h;

Par contre je vois pas comment tu pourrais éviter le tri (ici c'est la fonction unique qui s'en charge).

Bonne pioche. Ton code est plus de 10 fois plus rapide que le mien sur mes données

Je laisse le fil ouvert quelques heures au cas ou ...

[Jerome Briot]

N'arrivant pas bien à comprendre avec exactitude la forme des données de ton problème, voici quelques petites suggestions...

Code :

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

1
2
3
index = [1 2 3 2 5 2 6] ;
a = 0.2 ;
M = histc(index,1:10)*a;

ou

Code :

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

1
2
3
4
5
index = [1 2 3 2 5 2 6] ;
M = zeros(1,10);
a = 0.2 ;
v = 1:max(index);
M(v) = histc(index,v)*a;

ou encore

Code :

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

1
2
3
index = [1 2 3 2 5 2 6] ;
a = 0.2 ;
M = histc(index,1:max(index))*a;

Pour bien comprendre, la variable index contient les indices linéaires de certains éléments de la matrice M (1000x1000), c'est ça ?

[tubaas]

salut
une proposition différente mais moins efficace (en même temps c'est normal qu'une fonction comme HISTC soit plus optimisée...)

Code :

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

1
2
3
4
5
function M=toto(index,a)
M=a*arrayfun(@(x)countInd(x,index),1:10);
 
function y=countInd(x,index)
y=sum(index==x);

j'ai essayé différentes façons de compter mais les temps d'exécution restent similaires sur ton petit exemple