Discussion :

[LB_99]

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 :

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

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...

[bruno_pages]

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)

[LB_99]

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

Code :

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

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..

[bruno_pages]

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

Code :

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

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

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

[Joel F]

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

[3DArchi]

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

+++++++

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à :

Code :

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

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

Code :

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

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 :

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

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.