operateur éviter copie

Bonjour,

J'ai l'opérateur suivant. Mon problème est que mon code tourne très lentement.
Je penses que cela est due au fait que je créé un vecteur à chaque utilisation de *.
Mais je ne vois pas comment on peut l'éviter.

Code:

---

1 2 3 4 5 6 7 8 9 10 11

vecteur tridiag::operator*(const vecteur V){ int N=V.size(); vecteur Y(N);   Y(0) = ann[0]*V(0) + a1n[0]*V(1); for (int i=1; i<N-1; ++i) Y(i) = an1[i-1]*V(i-1)+ann[i]*V(i)+a1n[i]*V(i+1);   Y(N-1) = an1[N-2]*V(N-2)+ann[N-1]*V(N-1); return Y; }

---

Merci d'avance.

Citation:

---

Envoyé par xavierdestev

Je penses que cela est due au fait que je créé un vecteur à chaque utilisation de *.

---

La seule façon d'en être sur est le profiling.

Code:

---

vecteur tridiag::operator*(const vecteur V) {

---

Il y à ici une copie inutile, tu peux passer une ref en paramètre.

Tu peux aussi éviter la création d'un nouveau vecteur, si tu en possède déja un (si possible à la bonne taille)

Code:

---

1 2 3 4 5 6 7 8

void dot(const vecteur& V, vecteur& result) { int N=V.size(); if(result.size() != N) { result.resize(N); } result(0) = ...; ... }

---

Merci pour ta réponse.
Il y à ici une copie inutile, tu peux passer une ref en paramètre.

--->Si je fais cela j'obtiens, cpp/vecteur.cpp:57:17: attention : reference to local variable ‘U’ returned

Tu peux aussi éviter la création d'un nouveau vecteur, si tu en possède déja un (si possible à la bonne taille)

-->Je préfères éviter d'avoir recours à une méthode plutot qu'à la surdéfinition d'un opérateur.

Citation:

---

Envoyé par xavierdestev

cpp/vecteur.cpp:57:17: attention : reference to local variable ‘U’ returned

---

Code:

---

1 2 3 4 5 6 7 8 9 10 11

vecteur tridiag::operator*(const vecteur& V){ int N=V.size(); vecteur Y(N);   Y(0) = ann[0]*V(0) + a1n[0]*V(1); for (int i=1; i<N-1; ++i) Y(i) = an1[i-1]*V(i-1)+ann[i]*V(i)+a1n[i]*V(i+1);   Y(N-1) = an1[N-2]*V(N-2)+ann[N-1]*V(N-1); return Y; }

---

Le paramètre qui est passé par référence, pas la valeur de retour. ;)

Code:

---

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

vecteur tridiag::operator*(vecteur const& V) const{ assert(! V.empty()); vecteur Y(V.size()); // je tenterai en premier. Dans mes souvenirs, il y a ainsi plus de chances de permettre une copy-elision const size_t N=V.size() -1; // const n'est pas une optim; mais le -1 peut l'être   Y(0) = ann[0]*V(0) + a1n[0]*V(1); for (size_t i=1; i<N; ++i) Y(i) = an1[i-1]*V(i-1)+ann[i]*V(i)+a1n[i]*V(i+1);   Y(N) = an1[N-1]*V(N-2)+ann[N]*V(N); return Y; }

---

Prendre le vecteur de retour par référence ne vaut le coup que si la fonction est appelée de multiples fois avec la même variable en paramètre.
Il faut faire confiance au compilo sinon pour profiter de l'élision de copie et des rvalue-reférences sinon -- penser à écrire un constructeur de déplacement dans la classe vecteur.

Après, il est toujours possible de passer d'une notation indices à l'arithmétique de pointeurs : 1- si le profiling dit bien que le temps est perdu dans la fonction et 2- si le compilo est plus efficace avec cette technique plutôt que l'autre -- des fois ça change.