Discussion :

[Cpowa]

Bonjour,

J'ai la boucle suivante:

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

Je l'ai vectorisée comme ceci:

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

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.


Merci.

[Jean Dumoncel]

Bonjour,

dans ta dernière ligne il faut remplacer a et b par a2 et b2, non?

De toute manière, tu n'as pas besoin de ces variables, tu peux écrire directement :

Code :

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

xPoint=a(1:end-1)./b(2:end-1).*max(x(1:end-2)-x(2:end-1),0).^(3/2);

[Cpowa]

Bonjour,

dans ta dernière ligne il faut remplacer a et b par a2 et b2, non?

Oui, pardon.

De toute manière, tu n'as pas besoin de ces variables, tu peux écrire directement :

Code :

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

xPoint=a(1:end-1)./b(2:end-1).*max(x(1:end-2)-x(2:end-1),0).^(3/2);

Ok, merci.