Discussion :

[Décembre]

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 :

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
 
 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

[Invité]

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 :

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

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.

[Décembre]

Merci Winjerome

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

Code :

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

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?

[Invité]

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...

[Décembre]

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

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 :

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

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?

[Invité]

Il te manque l'utilisation de la fonction ANY:

Code :

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

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 :

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

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;

[Décembre]

Merci Winjerome