tableau de double - problème de remplissage valeur nulle

Version imprimable

Bonjour,

J'ai un problème avec une fonction un peu compliquée qui retourne un double, que je stocke dans un tableau. J'appelle la fonction dans une boucle.
Sauf qu'à certains moments où elle est censée me retourner 0, au lieu de mettre 0 dans mon tableau, elle recopie la valeur de la case précédente...

A noter que ma fonction est reprise d'une fonction originale qui elle fonctionne très bien. Ma seule modif de cette fonction originale : j'ai supprimé les fprintf et modifier les paramètres.

Voilà l'originale qui fonctionne :

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
void
RandomExcursionsVariant(int n)
{
	int		i, p, J, x, constraint, count, *S_k;
	int		stateX[18] = { -9, -8, -7, -6, -5, -4, -3, -2, -1, 1, 2, 3, 4, 5, 6, 7, 8, 9 };
	double	p_value;
 
	if ( (S_k = (int *)calloc(n, sizeof(int))) == NULL ) {
		fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\tRANDOM EXCURSIONS VARIANT: Insufficent memory allocated.\n");
		return;
	}
	J = 0;
	S_k[0] = 2*(int)epsilon[0] - 1;
	for ( i=1; i<n; i++ ) {
		S_k[i] = S_k[i-1] + 2*epsilon[i] - 1;
		if ( S_k[i] == 0 )
			J++;
	}
	if ( S_k[n-1] != 0 )
		J++;
 
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t\tRANDOM EXCURSIONS VARIANT TEST\n");
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t--------------------------------------------\n");
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\tCOMPUTATIONAL INFORMATION:\n");
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t--------------------------------------------\n");
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t(a) Number Of Cycles (J) = %d\n", J);
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t(b) Sequence Length (n)  = %d\n", n);
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t--------------------------------------------\n");
 
	constraint = (int)MAX(0.005*pow(n, 0.5), 500);
	if (J < constraint) {
		fprintf(stats[TEST_RND_EXCURSION_VAR], "\n\t\tWARNING:  TEST NOT APPLICABLE.  THERE ARE AN\n");
		fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t\t  INSUFFICIENT NUMBER OF CYCLES.\n");
		fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t---------------------------------------------\n");
		for ( i=0; i<18; i++ )
			fprintf(results[TEST_RND_EXCURSION_VAR], "%f\n", 0.0);
	}
	else {
		for ( p=0; p<=17; p++ ) {
			x = stateX[p];
			count = 0;
			for ( i=0; i<n; i++ )
				if ( S_k[i] == x )
					count++;
			p_value = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2))));
 
			if ( isNegative(p_value) || isGreaterThanOne(p_value) )
				fprintf(stats[TEST_RND_EXCURSION_VAR], "\t\t(b) WARNING: P_VALUE IS OUT OF RANGE.\n");
			fprintf(stats[TEST_RND_EXCURSION_VAR], "%s\t\t", p_value < ALPHA ? "FAILURE" : "SUCCESS");
			fprintf(stats[TEST_RND_EXCURSION_VAR], "(x = %2d) Total visits = %4d; p-value = %f\n", x, count, p_value);
			fprintf(results[TEST_RND_EXCURSION_VAR], "%f\n", p_value); fflush(results[TEST_RND_EXCURSION_VAR]);
		}
	}
	fprintf(stats[TEST_RND_EXCURSION_VAR], "\n"); fflush(stats[TEST_RND_EXCURSION_VAR]);
	free(S_k);
}

Voilà ma fonction dont le problème est décrit ci-dessus :

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
30
31
32
33
34
void excursionsVariantTest(char fichierATester[], int n, unsigned char *epsilon, double pvalues[])
{
	int		i, p, J, x, constraint, count, *S_k;
	int		stateX[18] = { -9, -8, -7, -6, -5, -4, -3, -2, -1, 1, 2, 3, 4, 5, 6, 7, 8, 9 };
 
	if ( (S_k = (int *)calloc(n, sizeof(int))) == NULL ) {
		exit(0);
	}
	J = 0;
	S_k[0] = 2*(int)epsilon[0] - 1;
	for ( i=1; i<n; i++ ) {
		S_k[i] = S_k[i-1] + 2*epsilon[i] - 1;
		if ( S_k[i] == 0 )
			J++;
	}
	if ( S_k[n-1] != 0 )
		J++;
 
	constraint = (int)MAX(0.005*pow(n, 0.5), 500);
	if (J < constraint) {
 
	}
	else {
		for ( p=0; p<=17; p++ ) {
			x = stateX[p];
			count = 0;
			for ( i=0; i<n; i++ )
				if ( S_k[i] == x )
					count++;
			pvalues[p] = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2))));
	}
	free(S_k);
}
}

Merci d'avance.

La ligne 45 d'un coté et 30 de l'autre sont différentes :
Code:

1 2 45 : pvalues[p] = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2)))); 30 : p_value = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2))));
Normal ?

12/07/2012, 15h25
douze_12

Oui c'est normal, au lieu d'écrire les valeurs dans un fichier je les stocke dans un tableau pvalues, du coup dans la version originale on a un double "pvalue" alors que moi j'ai mis un tableau de double "pvalues".

Merci.
12/07/2012, 17h09
Neckara
Bonjour,
Code:

1 2 3 if ( (S_k = (int *)calloc(n, sizeof(int))) == NULL ) { exit(0); }
Si calloc retourne NULL, c'est qu'il y a une erreur, pourquoi dans ce cas là retourner 0 qui signifie que le programme c'est déroulé correctement ?
Code:

1 2 3 4 if (J < constraint) { } else {
A quoi sert ce if ? Pourquoi ne pas mettre directement if(J >= constraint) ?
Code:

1 2 3 free(S_k); } }
Si J < constraint alors tu ne libère pas S_k, tu aurais alors une fuite de mémoire ?

Essaye de bien indenter ton code aussi.

Si j'ai bien compris c'est à la ligne :

Code:

pvalues[p] = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2))));

que tu as un problème? As-tu essayé d'afficher la valeur retournée par erfc ?
12/07/2012, 20h05
Trademark

Je pencherais pour un problème d'initialisation, après avoir vaguement lu le code original, il écrit parfois 0.0. Or si tu te retrouves dans ce cas, il se peut que ton tableau n'ait pas été initialisé à zéro et donc te renvoie des résultats différents. Essaye en l'initialisant.
13/07/2012, 17h06
douze_12

Merci pour vos réponses.
Même après initialisation, même problème.

Neckara, c'est bien cette ligne :

Code:

pvalues[p] = erfc(fabs(count-J)/(sqrt(2.0*J*(4.0*fabs(x)-2))));

J'ai ajouté un printf juste derrière, et en fait au lieu d'afficher 0, c'est comme s'il n'exécutait pas la fonction. Dans mon cas, la fonction est appelée dans une boucle

Code:

for ( i=0; i<4; i++ )

.

Dans la version originale, cela renvoie :

Citation:

0.708754
0.111570
0.000000
0.227275

Dans ma version, cela ne renvoie plus que 3 valeurs :

Citation:

0.708754
0.111570
0.227275

Il n'entre pas dans la boucle, peut-être à cause de la contrainte...?