tableau const, mais modifié ?!

Bonjour,
Je me suis remis récemment au C++, et j'ai quelques soucis avec le passage en argument des tableaux.

Pour résumer je code un classe Matrice, qui effectue un stockage creux, et j'ai un problème dans la fonction membre qui effectue le produit de cette matrice par un vecteur (tableau de double).

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29

void MATRICE::mult(const double * xx, double * vv) { int ii,j; double valx;   /*for(ii=0;ii<10;ii++) { cout<<xx[ii]<<" "; } cout<<endl;*/   for(j=0 ; j<n ; j++,vv++) *vv=0;   for(j=0;j<n;j++) { valx = xx[j]; for(ii=tableau_p[j]; ii<tableau_p[j+1];ii++) { vv[tableau_i[ii]]+= valx * tableau_valeur[ii]; } } /* for(ii=0;ii<10;ii++) { cout<<xx[ii]<<" "; }*/     }

---

Le problème est que mon tableau xx est modifié entre les deux affichages (qui sont commentés ici). Je ne vois pas trop pourquoi...

il suffit que vv pointe dans xx, soit avant l'appel, soit plus probablement à cause de l'incrémentation de vv lors de la mise à zéro qui semble louche et qui provoque un débordement dans la suite de la fonction (il aurait fallut écrire vv[j] = 0 et laisser vv inchangé, ou remplacer la boucle par un memset)

Merci, en fait j'avais mal alloué les tableaux :

Code:

---

1 2 3 4 5

vv = new double(n);   au lieu de   vv= new double[n];

---

Je ne comprend d'ailleurs pas pourquoi le code compilait..

le code compilait parce que vv = new double(n) est tout à faut légal, c'est l'équivalent de

Code:

---

1 2	vv = new double; *vv = n;

---

vous êtes sure que la modification de vv via le premier for est correcte ?

et bon, std::vector et ces problemes disparaissent ...

Citation:

---

Envoyé par Joel F

et bon, std::vector et ces problemes disparaissent ...

---

+++++++

Citation:

---

Envoyé par LB_99

Le problème est que mon tableau xx est modifié entre les deux affichages (qui sont commentés ici). Je ne vois pas trop pourquoi...

---

Ecrasement du à un problème d'indice a été ma première idée. Du coup, j'ai lu ton code et là :

Citation:

---

Envoyé par LB_99

Code:

---

1 2 3 4 5 6 7 8 9 10 11

for(j=0 ; j<n ; j++,vv++) *vv=0;   for(j=0;j<n;j++) { valx = xx[j]; for(ii=tableau_p[j]; ii<tableau_p[j+1];ii++) { vv[tableau_i[ii]]+= valx * tableau_valeur[ii]; } }

---

---

Comme le dit bruno, tu as modifié ton pointeur (pas la valeur pointée, mais l'adresse) de vv. Si n est la taille allouée pour vv, alors le pointeur désigne la fin de la zone mémoire valide et va commencer à jardiner :aie:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

for(j=0 ; j<n ; j++) vv[j]=0.;   for(j=0;j<n;j++) { valx = xx[j]; for(ii=tableau_p[j]; ii<tableau_p[j+1];ii++) { vv[tableau_i[ii]]+= valx * tableau_valeur[ii]; } }

---

Note que const peut aussi s'utiliser sur l'adresse d'un pointeur et non sur la valeur pointée. Tu peux alors surveiller que tu ne changes pas cette adresse dans ta fonction :

Code:

---

void fonction(double const * const xx, double * const vv)

---

Enfin, pour mettre au point, je mettrai des assert sur les indices pour détecter ces erreurs.