Discussion :

[étoile de mer]

Bonsoir à tous,
J'ai un probleme avec ma fonction :

Code :

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

1
2
3
4
5
6
7
8
9
10
double info_mut(ChainonMot const *pc, ChainonMot const *pc2,int nblignes)
{
 double  IM=0;
 
 
   IM=((nbrLignesCommunes(pc,pc2)/nblignes)*log2((nbrLignesCommunes(pc,pc2)/nblignes)/(nombreLigne(pc)*nombreLigne(pc2))));
printf("%lf",IM);
//printf("%d",(nombreLigne(pc)*nombreLigne(pc2)));
    return  IM;
}

LE compilateur m'affiche NAN au lieu d'une certaine valeur numérique.
J'ai pensé à ce qu'il ya un souci de division par zéro, mais j'ai bien verifié que nombreLigne(pc)*nombreLigne(pc2) est différetnte de 0
ce que je trouve que bizarre
Une idée à cela svp?
Merci

[Jerome Briot]

A mon avis, le problème vient du terme (nbrLignesCommunes(pc,pc2)/nblignes) qui doit être nul => 0/0 => NaN

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
>> log2(0)
 
ans =
 
  -Inf
 
>> log2(0)/0
 
ans =
 
  -Inf
 
>> 0*log2(0)/0
 
ans =
 
   NaN

[étoile de mer]

Merci Dut,
aussi je vois qu'il ya un probleme à ce niveau

Code :

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

nbrLignesCommunes(pc,pc2)/nblignes

qui m'affiche égal à 0.000000
or quand je fais

Code :

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

printf("%i",nbrLignesCommunes(pc,pc2)), il m'affiche 3

et nblignes = 2
donc normallement ca vaut 0.5 et non pas 0.00000

[SofEvans]

Merci Dut,
aussi je vois qu'il ya un probleme à ce niveau

Code :

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

nbrLignesCommunes(pc,pc2)/nblignes

qui m'affiche égal à 0.000000
or quand je fais

Code :

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

printf("%i",nbrLignesCommunes(pc,pc2)), il m'affiche 3

et nblignes = 2
donc normallement ca vaut 0.5 et non pas 0.00000

C'est parce que tu fais une division d'entier et non une division de reel.
Cela veut dire que nbrLignesCommunes(pc,pc2) et nblignes sont des int.
Fais un simple cast pour avoir une division comme tu le souhaite.

Code :

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

((double) nbrLignesCommunes(pc,pc2))/nblignes

[étoile de mer]

Merci SofEvans,
mais même si je modifie la formule de cette façon :

Code :

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

1
2
3
4
5
6
7
8
9
double info_mut(ChainonMot const *pc, ChainonMot const *pc2,int nblignes)
{
 double  IM=0;
 
IM= (double)log2(nbrLignesCommunes(pc,pc2)/nblignes);
printf("%d  \n",IM);
 
    return  IM;
}

j'obtiens une valeur =0 au lieu de 0.6

[SofEvans]

Mauvais cast.

Si on regarde ceci :

Code :

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

1
2
 
IM= (double)log2(nbrLignesCommunes(pc,pc2)/nblignes);

on a ceci

Code :

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

1
2
 
(double)       log2(nbrLignesCommunes(pc,pc2)/nblignes)

Ce qui revient implicitement a ceci :

Code :

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

1
2
(double)    (   log2(nbrLignesCommunes(pc,pc2)/nblignes)   )

En resumé, tu demande au resultat de la division entiere de devenr un double.
Donc 0 deviendra 0.000, ce qui n'arrange pas l'affaire

Tu dois forcer le type en double sur une des deux composante de ta division.
On obtient ceci :

Code :

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

1
2
IM= log2( ((double)nbrLignesCommunes(pc,pc2))/nblignes);

ou ceci:

Code :

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

1
2
IM= log2(nbrLignesCommunes(pc,pc2)/((double)nblignes);

[étoile de mer]

et avec cette nouvelle formule

Code :

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

IM=((double)nbrLignesCommunes(pc,pc2) /nblignes)*log2(double)((nbrLignesCommunes(pc,pc2) /nblignes)/((nombreLigne(pc))*nombreLigne(pc2)));

il maffiche une errreeur
: expected expression before ‘double’|
Une idée?

[Jerome Briot]

Un petit problème de parenthèse :

Code :

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

log2((double)

[étoile de mer]

MErci Dut,
MAis le probleme de -INF persiste encore avec cette formule

Code :

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

IM=((double)nbrLignesCommunes(pc,pc2) /nblignes)*log2((double)((nbrLignesCommunes(pc,pc2) /nblignes)/((nombreLigne(pc))*nombreLigne(pc2))));

et le problème c'est que je suis sure que la valeur de nblignes est >0, donc d'où est ce que vient ce probleme de division par 0

[SofEvans]

Il y a deux division.

A ce stade, je te conseille de faire un printf de tout tes parametre pour eventuellement voir qui pourrait deranger. Sinon, le mieux serait de refaire la formule pas a pas, en parenthesant petit a petit pour eviter de coder la mauvaise formule.

[étoile de mer]

voilà je décortique ma fonction:

Code :

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

1
2
3
4
5
6
double  probabilit(ChainonMot const *pc, ChainonMot const *pc2,int nblignes)
{
    double  p=0;
    p=nbrLignesCommunes(pc,pc2)/nblignes;
    return p;
}

quand j'affiche la valeur de nbrLignesCommunes(pc,pc2) il m'affiche 2
aussi pour nblignes=3
donc il doit m'afficher 0.6666666
or là il m'affiche 0.0000000

[SofEvans]

Comme dit depuis toujours, tu fais une division d'entier !

Quant tu fais :

p=nbrLignesCommunes(pc,pc2)/nblignes;


En realite, p prend le nombre de fois que nbrLignesCommunes(pc,pc2) contient nblignes.

Plus explicitement, combien de fois 2 peut contenir 3 ? reponse : 0

Il est tout a fait normal que p, meme s'il est de type double prenne 0 et non 0.6666.

Pour resoudre ce probleme, meme conseil, meme remarque : caster un de deux membre de la division afin que la division soit une division de reel. Tu aura le resultat escompté.

[SofEvans]

Suite a MP avec cyrine, voici la formule :

IM(x,y)= probabilité (x,y)*LOG2(probabilité (x,y)/(probabilité(x )*probabilité(y )))

Si nous reprenons comme cela :





probabilité (x,y) -> obtenu grace a "nbrLignesCommunes(pc,pc2)/nblignes " (type INT)
probabilité(x) -> obtenue grace à nombreLigne(pc) (type INT)
probabilité(y) -> obtenue grace à nombreLigne(pc2) (type INT)

IM -> resultat (type DOUBLE)



Si nous reprenons calmement. ON va commencer par faire le log2

LOG2(probabilité (x,y)/(probabilité(x )*probabilité(y )))
reviendrait à
log2( (nbrLignesCommunes(pc,pc2)/nblignes) / (nombreLigne(pc)*nombreLigne(pc2)))

Seul probleme, les type :

log2((INT/INT)/(INT*INT)).
Ce qui peut s'ecrire :
log2(INT/INT)
Cela fait donc une division d'entier qui donne .... un INT. On a donc :

log2(INT)
Or, log2 n'est pas prevu pour avoir un INT, il fait donc une conversion implicite. De plus, tu perd de la precision. Il faut donc caster :

log2( (nbrLignesCommunes(pc,pc2)/((double)nblignes)) / (nombreLigne(pc)*nombreLigne(pc2)))


On a donc :
log2 ((INT/DOUBLE)/(INT*INT))
log2 (DOUBLE/INT)
log2 (DOUBLE)

On a donc ceci, et cela sera de type DOUBLE

Si on continue sur la formule :

IM(x,y)= probabilité (x,y)*LOG2(probabilité (x,y)/(probabilité(x )*probabilité(y )))

On a :

IM = (nbrLignesCommunes(pc,pc2)/nblignes) * log2((nbrLignesCommunes(pc,pc2)/((double)nblignes)) / (nombreLigne(pc)*nombreLigne(pc2))).

Selon, les type, on
DOUBLE = INT/INT * DOUBLE

Il faut encore une fois caster :

IM = (nbrLignesCommunes(pc,pc2)/((double)nblignes)) * log2((nbrLignesCommunes(pc,pc2)/((double)nblignes)) / (nombreLigne(pc)*nombreLigne(pc2))).


Cela devrait etre correcte.

La formule est normalement celle-ci :

Code :

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

1
2
 
IM = (nbrLignesCommunes(pc,pc2)/((double)nblignes))  *  log2((nbrLignesCommunes(pc,pc2)/((double)nblignes)) / (nombreLigne(pc)*nombreLigne(pc2)));

Cela marche-t'il ?

[étoile de mer]

Super, ca marche
Merci