[parfor] Où puis-je l'utiliser correctement dans mon algorithme ?

**Globoxx** · 23/11/2013, 19h40

Bonsoir,

J'ai écris un algorithme qui permet, en gros, de trouver le point le plus proche de chaque point dans l'espace.
Le fait est que cet algorithme est lui-même imbriqué dans 2 boucles for. Au final le temps de calcul se compte en heures.
J'aimerais savoir où je pourrais utiliser [parfor] dans mon code.

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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
 
% Algorithme
 
        VectorO = zeros(N,4);
        NormesO = zeros(N,1);
        All_Best_ways = zeros(N,1);
        Point = zeros(N,3);
 
        for i = 1 : N
            VectorO(i,1:3) = P(i,1:3) - Origine(1,1:3);
            NormesO(i) = sqrt((VectorO(i,1))^2+(VectorO(i,2))^2+(VectorO(i,3))^2);
        end
        Emplacement_Best_way = find(NormesO == min(NormesO(:)));
        if length(Emplacement_Best_way) ~= 1 % Si il y a plusieurs plus courts chemins
            Emplacement_Best_way(2:end) = []; % On ne prend que le premier
        end
        Best_way = min(NormesO);
        All_Best_ways(1) = Best_way;
        P(Emplacement_Best_way,4) = 1; % Le système est visité
        Point(1,:) = P(Emplacement_Best_way,1:3);
 
        if Point(1,:) == Terre(1:3) % Si la Terre est atteinte
            Jusqua_Terre = sum(All_Best_ways);
        end
 
        Vector = zeros(N-1,4);
        Normes = zeros(N-1,1);
        for j = 2 : N
            for i = 1 : N
                if i ~= Emplacement_Best_way
                    if P(i,4) == 0
                        Vector(i,1:3) = P(i,1:3) - P(Emplacement_Best_way,1:3);
                        Normes(i) = sqrt((Vector(i,1))^2+(Vector(i,2))^2+(Vector(i,3))^2);
                    else
                        Normes(i) = inf;
                    end
                else
                    Normes(i) = inf;
                end
            end
            Emplacement_Best_way = find(Normes == min(Normes(:)));
            if length(Emplacement_Best_way) ~= 1 % Si il y a plusieurs plus courts chemins
                Emplacement_Best_way(2:end) = []; % On ne prend que le premier
            end
            Best_way = min(Normes);
            All_Best_ways(j) = Best_way;
            P(Emplacement_Best_way,4) = 1; % Le système est visité
            Point(j,:) = P(Emplacement_Best_way,1:3);
 
            if Point(j,:) == Terre(1:3) % Si la Terre est atteinte
                Jusqua_Terre = sum(All_Best_ways);
                break % On arrête l'algorithme
            end
        end

Merci d'avance pour votre aide !

**Jerome Briot** · 25/11/2013, 12h02

Avant d'utiliser parfor, il y a quelques améliorations à apporter à ton code.

Le bloc suivant :

Code :

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

1
2
3
4
5
6
Emplacement_Best_way = find(Normes == min(Normes(:)));
if length(Emplacement_Best_way) ~= 1 % Si il y a plusieurs plus courts chemins
    Emplacement_Best_way(2:end) = []; % On ne prend que le premier
end
Best_way = min(Normes);
All_Best_ways(j) = Best_way;

devrait simplement s'écrire :

Code :

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

[All_Best_ways(j),Emplacement_Best_way] = min(Normes);

As-tu besoin des valeurs stockées dans Vector à la fin de ton code ?

Si ce n'est pas le cas, tu n'as pas besoin de lui attribuer une taille Nx3. Il suffirait de faire comme ceci :

Code :

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

1
2
Vector = P(i,1:3) - P(Emplacement_Best_way,1:3);
Normes(i) = sqrt(Vector(1)^2+Vector(2)^2+Vector(3)^2);

Pourquoi mets tu N-1 éléments dans la ligne suivante :

Code :

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

Normes = zeros(N-1,1);

Ne serait-ce pas plutôt N ?

Ensuite, tu remplis Normes avec des valeurs Inf. Initialise donc ce vecteur avec la fonction inf :

Code :

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

Normes = inf(N-1,1);

Ceci te permettra de simplifier la boucle interne sur i :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
for i = 1 : N
        if i ~= Emplacement_Best_way
            if P(i,4) == 0
                Vector(i,1:3) = P(i,1:3) - P(Emplacement_Best_way,1:3);
                Normes(i) = sqrt((Vector(i,1))^2+(Vector(i,2))^2+(Vector(i,3))^2);
            else
                Normes(i) = inf;
            end
        else
            Normes(i) = inf;
        end
    end

Ensuite, je ne suis pas certain que la condition suivante fasse ce que tu espères :

Code :

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

if Point(j,:) == Terre(1:3)

Je ferais plutôt

Code :

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

if all(Point(j,:) == Terre(1:3))

ou

Code :

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

if Point(j,1) == Terre(1) && Point(j,2) == Terre(2) && Point(j,3) == Terre(3)

**Globoxx** · 25/11/2013, 13h32

Merci beaucoup pour ta réponse !

J'ai donc simplifié mon code avec tes 3 premières propositions et cela fonctionne parfaitement.
Par contre, pour ce qui est d'écrire :

Code :

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

Normes = inf(N,1);

Mon code ne fonctionne plus correctement, après quelques itérations, la valeur de ligne i de All_Best_ways reste la même jusqu'à la fin.

Enfin, la condition :

Code :

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

if Point(j,:) == Terre(1:3)

Me permet de tester si le point j actuel se trouve à certaines coordonnées.
Si je mets un all, il va tester toutes les valeurs présentes dans Point et pas le point actuel non ?

**Jerome Briot** · 25/11/2013, 14h02

Envoyé par Globoxx

Par contre, pour ce qui est d'écrire :

Code :

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

Normes = inf(N,1);

Mon code ne fonctionne plus correctement, après quelques itérations, la valeur de ligne i de All_Best_ways reste la même jusqu'à la fin.

Bizarre, il n'y a pas de raison à cela

Envoyé par Globoxx

Me permet de tester si le point j actuel se trouve à certaines coordonnées.

Essai ceci :

Code :

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

1
2
3
4
>> if [1 2 3] == [1 2 3], disp('OK'); end
OK
>> if [0 2 3] == [1 2 3], disp('OK'); end
>> if [1 4 5] == [1 2 3], disp('OK'); end

Le test ne vérifie que l'égalité des premières valeurs de chaque terme.
Il faut plutôt utiliser all ou any selon le besoin

Envoyé par Globoxx

Si je mets un all, il va tester toutes les valeurs présentes dans Point et pas le point actuel non ?

Il va seulement tester les valeurs à la ligne j

**Jerome Briot** · 25/11/2013, 14h07

Peux-tu poster le code avec les modifications ?

**Globoxx** · 25/11/2013, 17h13

Le voici :

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
36
37
38
39
40
41
42
43
44
45
46
47
 
% Algorithme
 
        NormesO = zeros(N,1);
        All_Best_ways = zeros(N,1);
        Point = zeros(N,3);
 
        for i = 1 : N
            VectorO = P(i,1:3) - Origine(1,1:3);
            NormesO(i) = sqrt((VectorO(1))^2+(VectorO(2))^2+(VectorO(3))^2);
        end
        Emplacement_Best_way = find(NormesO == min(NormesO(:)));
        if length(Emplacement_Best_way) ~= 1 % Si il y a plusieurs plus courts chemins
            Emplacement_Best_way(2:end) = []; % On ne prend que le premier
        end
        Best_way = min(NormesO);
        All_Best_ways(1) = Best_way;
        P(Emplacement_Best_way,4) = 1; % Le système est visité
        Point(1,:) = P(Emplacement_Best_way,1:3);
 
        if all(Point(1,:) == Terre(1:3)) % Si la Terre est atteinte
            Jusqua_Terre = sum(All_Best_ways);
        end
 
        Normes = zeros(N,1);
        for j = 2 : N
            for i = 1 : N
                if i ~= Emplacement_Best_way
                    if P(i,4) == 0
                        Vector = P(i,1:3) - P(Emplacement_Best_way,1:3);
                        Normes(i) = sqrt((Vector(1))^2+(Vector(2))^2+(Vector(3))^2);
                    else
                        Normes(i) = inf;
                    end
                else
                    Normes(i) = inf;
                end
            end
            [All_Best_ways(j),Emplacement_Best_way] = min(Normes);
            P(Emplacement_Best_way,4) = 1; % Le système est visité
            Point(j,:) = P(Emplacement_Best_way,1:3);
 
            if all(Point(j,:) == Terre(1:3)) % Si la Terre est atteinte
                Jusqua_Terre = sum(All_Best_ways);
                break % On arrête l'algorithme
            end
        end

Finalement all fonctionne très bien merci !
Mais pour le Normes = inf(N,1)... Je sèche.

[parfor] Où puis-je l'utiliser correctement dans mon algorithme ?

MATLAB

Vue hybride

Discussions similaires

Partager

Partager