Bonjour,
J'ai la boucle suivante:
Je l'ai vectorisée comme ceci:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4 for n=2:N-1 xPoint(n)=a(n-1)/b(n)*max(x(n-1)-x(n),0)^(3/2); end
Est-ce la meilleur façon de faire ? Par que c'est certes plus rapide que la boucle for mais ça gache de la mémoire avec les matrices intermédiaires.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6 a2=a(1:end-1); b2=b(2:end-1); x2=x(1:end-2); x3=x(2:end-1); xPoint=a./b.*max(x2-x3,0).^(3/2);
Merci.
Partager