Discussion :

[douze_12]

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 :

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

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

[Bktero]

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

Code :

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

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 ?

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

[Neckara]

Bonjour,

Code :

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

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 :

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

1
2
3
4
if (J < constraint) {
 
	}
        else {

A quoi sert ce if ? Pourquoi ne pas mettre directement if(J >= constraint) ?

Code :

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

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 :

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

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 ?

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

[douze_12]

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

Neckara, c'est bien cette ligne :

Code :

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

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 :

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

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

.

Dans la version originale, cela renvoie :

0.708754
0.111570
0.000000
0.227275

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

0.708754
0.111570
0.227275

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