comment se passer des boucles for

Bonjour,

J'ai un grand handicap avec matlab !!! je ne sais pas comment éviter les boucles for qui sont très gourmandes en temps de calcul :aie:

J'essaie désespérément de remplacer cette boucle :

Code:

---

1 2 3 4 5 6 7 8 9 10

P=[x y z]; for i=1:length(w) S=find(L==w(i)); xp(j)= P(S(1),1,:); yp(j)= P(S(1),2,:); zp(j)= P(S(1),3,:); j=j+1;   end

---

mais je n'y arrive pas :arf:

Pouvez-vous m'orienter sans me donner la solution pour que je puisse enfin me passer une bonne fois pour toute des boucles for :help:

Il ne faut pas tout de suite vouloir supprimer la boucle.

En premier lieu, il faut chercher à l'optimiser.

Voici deux pistes :

1. Utiliser la pré-allocation de mémoire
2. Bien lire la documentation de la fonction FIND qui peut prendre plus d'un argument d'entrée

Merci pour ta réponse,

Citation:

---

Utiliser la pré-allocation de mémoire

---

Je ne connais pas exactement la taille de mes matrices à l'avance !!!

Citation:

---

Envoyé par Décembre

Je ne connais pas exactement la taille de mes matrices à l'avance !!!

---

Dans le code que tu nous montres, leurs tailles sont connues... d'ailleurs la variable j ne sert à rien :?

Sinon, il arrive parfois que l'on ne puisse pas connaitre la taille des variables à l'avance. Mais on peut estimer une taille maximale. Et donc faire la pré-allocation sur cette taille maximale. Les variables finales contiendront donc quelques données inutiles mais faciles à gérer.

Citation:

---

Envoyé par Décembre

Pouvez-vous m'orienter sans me donner la solution pour que je puisse enfin me passer une bonne fois pour toute des boucles for

---

Et si tu nous disais que représente les variables w, L x y z ? et que veux tu faire avec cette boucle ?

une copie de "whos" par exemple !

Voici la copie de whos

Code:

---

1 2 3 4 5 6 7

Name Size Bytes Class x 71385x1 571080 double array y 71385x1 571080 double array z 71385x1 571080 double array w 65x1 520 double array L 1x71385 571080 double array

---

Citation:

---

Dans le code que tu nous montres, leurs tailles sont connues... d'ailleurs la variable j ne sert à rien

---

Si j'enlève le j je risque d'écraser les valeurs de xp, yp et zp.

Merci beaucoup pour votre aide :fleur:

Citation:

---

Envoyé par Décembre

Si j'enlève le j je risque d'écraser les valeurs de xp, yp et zp.

---

Oui mais tu peux bien remplacer j par i, non ?

Citation:

---

Oui mais tu peux bien remplacer j par i, non ?

---

Ouiiiiiiii tu as raison :merci:

J'aurais aimée me passer de la précédente boucle for mais apparemment c'est impossible :traine:

Dans le but d'optimiser mon code, je vais tenter de remplacer la boucle suivante:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

for i=1:3:length(L) a=[L(i) L(i+1) L(i+2)]; a=sort(a); M(j,:)=[a(1) a(2)]; j=j+1; M(j,:)=[a(1) a(3)]; j=j+1; M(j,:)=[a(2) a(3)]; j=j+1; end

---

Pouvez-vous m'orienter une autre fois :calim2:

Tu traites les valeurs de L par paquets de 3... dans ce cas, redimensionne L en une matrice 3xN (ou Nx3) appliques SORT sur la matrice en dehors de la boucle

Deux choses :

* Note bien ce que t'as dit Dut : la fonction FIND peut etre utilisee avec plus qu'un seul argument ;-)

* Pourquoi ecrire P(S(1), 1, : ) alors qu'a priori P a seulement 2 dimensions ? :roll:

Et enfin, tu connais la fonction ACCUMARRAY ? Elle peut t'etre utile...

Salut,

Je n'arrive pas à trouver ce qui cloche tout parait parfait pour moi mais le résultat est différent !!

voici ce que j'ai fait en suivant tes précieux conseils monsieur Dut:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

L1=reshape(L,length(L)/3,3); a=sort(L1,2); j=1; for i=1:length(L)/3 M(j,:)=[a(i,1) a(i,2)]; j=j+1; M(j,:)=[a(i,1) a(i,3)]; j=j+1; M(j,:)=[a(i,2) a(i,3)]; j=j+1; end

---

ce qui me trouble c'est que la matrice M que j'obtiens n'est pas la même 8O

si vous arrivez à trouver ou ça cloche je vous prie de m'aider :help:

Remplace

Citation:

---

Envoyé par Décembre

Code:

---

1 2	L1=reshape(L,length(L)/3,3);

---

---

par

Code:

---

1 2	L1=reshape(L,3,length(L)/3).';

---

Tu peux aussi te passer de l'opérateur de transposition en modifiant les indices dans la suite du code.

Merci Dut :applo: ça marche grâce à toi:mrgreen:



J'ai gagné 4 secondes, mais ce n'est encore pas suffisant alors peux-tu m'expliquer s'il te plait que veux-tu dire par là:

Citation:

---

Tu peux aussi te passer de l'opérateur de transposition en modifiant les indices dans la suite du code.

---

merci encore :fleur:

C'est assez difficile de parler d'optimisation sans connaitre ni la nature exacte, ni les dimensions du problème ...

Essaye toujours ceci :

Code:

---

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

L1 = reshape(L,3,[]); a = sort(L1,1); j = 1; M = zeros(numel(L),2); for i = 1:size(a,2) M(j,1) = a(1,i) ; M(j,2) = a(2,i); j=j+1; M(j,1) = a(1,i); M(j,2) = a(3,i); j=j+1; M(j,1) = a(2,i); M(j,2) = a(3,i); j=j+1; end

---

Je ne sais pas comment te remercier Dut :merci:

Cette boucle s'exécutait en 60 secondes et grâce à toi elle ne dure que 5 secondes :lahola:

Merci,merci, merci,merci, merci !!!

:fleur: