Discussion :

[jozog47]

Bonjour à tous.

Je rencontre un petit souci dans le passage d'un tableau 2D en paramètres.

Je pars de ça dans mon main :

Code :

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

1
2
int ** R = malloc(m*sizeof(double*));
for ( l=0 ; l<m ; l++ ) { R[l]=malloc(m*sizeof(double)); }

Juste avant d'appeler une fonction qui travaille sur R, un printf("%d\n",R[0][0]); me renvoie 3.

Ensuite, je fais appel à une de mes fonctions (voici l'appel et le code de la dite fonction) :

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
double ** L = Choleski(R,m);
 
/* Pour alléger, tout le contenu du main n'est pas affiché.
Ici, je travaille un peu sur le tableau R, j'assigne des valeurs
à tous ses éléments avant de l'expédier à la fonction Choleski,
fonction source de mon problème vraisemblablement. */
 
double ** Choleski(int ** A, int m)
{
	int i, j, k;
	double sum;
 
	double ** L = malloc(m*sizeof(double*));
	for ( k=0 ; k<m ; k++ ) { L[k]=malloc(m*sizeof(double)); }
 
	// Initialization of L (0 above the diagonal of the matrix)
	for ( j=1 ; j<m ; j++ )
	{
		for ( i=0 ; i<j ; i++ ) { L[i][j]=0.0; }
	}
 
	// First colum of L
	L[0][0] = sqrt((double)A[0][0]);
        printf("%e\n",A[0][0]); /* CE PRINTF AFFICHE UN NOMBRE EN 10^-314 */
	for ( j=1 ; j<m ; j++ )
	{
		L[j][0] = ((double)A[0][j])/L[0][0];
	}
 
	// Other columns
	for ( i=1 ; i<m ; i++ )
	{
		printf("%d iterations remaining for Choleski-descomposition of RP.\n",m-1-i);
		sum=0.0;
		for ( k=0 ; k<i ; k++ ) { sum += pow(L[i][k],2.0); }
		L[i][i] = sqrt(((double)A[i][i])-sum);
		for ( j=i+1 ; j<m ; j++ )
		{
			sum=0.0;
			for ( k=0 ; k<i ; k++ ) { sum += L[i][k]*L[j][k]; }
			L[j][i] = (((double)A[i][j])-sum)/L[i][i];
		}
	}
 
	return L;
}

Le printf mis en évidence à la ligne 19 dans le code de la fonction Choleski me donne un nombre de l'ordre de 10^(-314), au lieu de 3... Tous les autres coefficients de ma matrice R, une fois passés en paramètres ont également cette valeur.
Il s'en suit pas mal de problèmes de NaN dans les divisions du reste de cette fonction.

Quelqu'un saurait me dire où se trouve mon erreur ?

Pas de lien vers la FAQ s'il vous plaît, j'avoue qu'elle ne m'a pas trop aidé et je souhaiterais savoir qu'est-ce qui, dans ce code, fait qu'il ne fonctionne pas...

Merci d'avance.

[Invité(e)]

Bonjour,
Est-ce voulu de définir R comme pointeur de pointeur sur int et de faire les allocations avec le type double ?

Pour éviter ce genre de confusion, on peut écrire

Code :

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

1
2
3
4
int ** R = malloc(m * sizeof *R);
for ( l=0 ; l<m ; l++ ) { 
    R[l] = malloc( m * sizeof *R[l]); 
}

[jozog47]

Ah, non, ce n'était pas voulu. En manipulant les deux, je me suis un peu embrouillé, c'est corrigé.

Dans ton petit algo mabu, sizeof *R est équivalent à sizeof(int) ?

La grosse erreur était dans mon printf aussi : j'ai utilisé un %e mais, comme le nombre étant déclaré en int, le chiffre affiché était aberrant.

[Invité(e)]

Dans ton petit algo mabu, sizeof *R est équivalent à sizeof(int) ?

sizeof *R va renvoyer la taille nécessaire au stockage d'un élément stocké dans R.
Si tu changes ce qui est pointé par R, tu n'as pas à modifier ton malloc

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
double *D;
int *I;
 
/* sizeof *D <=> sizeof(double) */
/* sizeof *I <=> sizeof(int) */;
 
/*allocation d'un tableau de 10 éléments double et 
  d'un tableau de 10 éléments int */
D = malloc(10*sizeof *D);
I = malloc(10*sizeof *I);

[jozog47]

Très bonne astuce, il faudra que j'y pense à l'avenir.