Optimisation d'un bout de code

Bonjour,
Je traite un nombre important de données et je me trouve confrontée à un temps de calcul assez élevé, alors je me demandais si on ne pouvait pas mieux faire.

Voici la partie qui me bouffe le plus de temps:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

P=rand(65000,3); V1=rand(90,3); for i=1:3:length(P)-2 gi=zeros(3,3); gi(1,:)=P(i,:); gi(2,:)=P(i+1,:); gi(3,:)=P(i+2,:); [ci, ia, ib]=intersect(gi, V1,'rows'); if size(ia)~=0 P(i,:)= [0 0 0]; P(i+1,:)=[0 0 0]; P(i+2,:)=[0 0 0]; end   end

---

Je répète se code au minimum 100 fois alors j'aimerais bien savoir s'il y a moyen de faire autrement en gaspillant moins de temps?

Merci d'avance :fleur:

Bonjour,

Utilise plutôt la fonction NUMEL ou ISEMPTY à la place de SIZE dans if size(ia)~=0Ensuite tu peux simplifier ta boucle comme ceci:

Code:

---

1 2 3 4 5 6

for i=1:3:length(P)-2 [ci, ia, ib] = intersect(P(i:i+2,:), V1,'rows'); if numel(ia)~=0 P(i:i+2,:) = 0; end end

---

Sinon pour optimiser le temps d'exécution, je te conseille d'appliquer la fonction INTERSECT sur tout le vecteur P et de travailler ensuite sur les indices, notamment avec FLOOR.
Bonne continuation.

Merci Winjerome :ccool:

En suivant ton conseil j'ai remplacé par ceci:

Code:

---

1 2 3 4 5

[ci, ia, ib]=intersect(P, V1,'rows'); if numel(ia)~=0 P(ia,: )=0; end

---

Mais j'ai pas utilsé Floor et je ne sais pas si ce que j'ai fait est correct ou pas, qu'en pensez vous?

Citation:

---

Envoyé par Décembre

je ne sais pas si ce que j'ai fait est correct ou pas, qu'en pensez vous?

---

Cela dépend de tes besoins...
Dans le code de ton premier message, tu fais la vérification 3 par 3 (si au moins l'un des 3 est contenu dans V1, tu mets les 3 à zéro), alors qu'ici tu n'en tiens pas compte...

Citation:

---

Dans le code de ton premier message, tu fais la vérification 3 par 3 (si au moins l'un des 3 est contenu dans V1, tu mets les 3 à zéro), alors qu'ici tu n'en tiens pas compte

---

Effectivement, en faidsant vite j'ai omis de prendre en considération la vérification 3 par 3 qui est très importante :oops:

J'ai essayé ceci, mais le résultat que j'obtiens est différent de ce que j'avais avant (en utilisant le code lent), il y a certainement une erreur mais j'arrive pas à la voir !

Code:

---

1 2 3 4 5 6 7 8

[ci, ia, ib]=intersect(P, V1,'rows');   for i=1:3:length(P)-2 if (ia(:)==i)|(ia(:)==i+1)|(ia(:)==i+2) P(i:i+2,:)=0; end end

---

Voyez-vous le problème?

Il te manque l'utilisation de la fonction ANY:

Code:

---

if any(ia(:)==i)|(ia(:)==i+1)|(ia(:)==i+2)

---

Sinon tu vérifies si toutes les valeurs du vecteur sont vraies.

Ou comme indiqué précédemment avec FLOOR:

Code:

---

1 2 3 4 5

[ci, ia, ib]=intersect(P, V1,'rows'); ia = 3*floor((ia-1)/3)+1; P(ia,:) = 0; P(ia+1,:) = 0; P(ia+2,:) = 0;

---

Merci Winjerome :lahola: