Lenteur d'exécution et bug programme

Bonjour,

J'utilise un code qui me permet en connaissant la matrice de rotation (3x3) et celle de translation (vecteur 3x1) de recaler des points sur des autres. Le problème est que matrice de points Pinterb est de taille 3x70000 (points en 3D) et que lorsque j'utilise la transformation pour recaler chacun de mes points dans une boucle for, le programme rame et plante finalement... Voici le code que j'utilise ci-dessous :

Code:

---

1 2 3

for i=1:size(Pinterpb,2) %Calcul des points recalés Q(:,i)=inv(R)*(Pinterpb(:,i)-q(5:7)); end

---

Pinterpb est ma matrice de points 3x70000 à recaler, R la matrice de rotation et q le vecteur translation : c'est une transformation inverse en fait que j'applique.

Merci de votre aide s'il existe un moyen pour optimiser et ne pas faire planter l'ordi !

1ère chose à faire : as-tu penser à pré-allouer ta matrice :

Code:

---

1 2	Q=zeros(size(Pinterb));

---

Merci, impeccable :king: , j'avais oublié de le faire pour celle-ci...

2ème chose :

Code:

---

1 2 3 4 5 6

invR=inv(R); Q=zeros(size(Pinterb)); for i=1:size(Pinterpb,2) %Calcul des points recalés Q(:,i)=invR*(Pinterpb(:,i)-q(5:7)); end

---

Ou mieux, l'utilisation de l'opérateur \ (backslash)