Probleme avec NaN et -Inf

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

Code:

---

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

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

Code:

---

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

---

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

Code:

---

nbrLignesCommunes(pc,pc2)/nblignes

---

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

Code:

---

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

---

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

Citation:

---

Envoyé par cyrine

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

Code:

---

nbrLignesCommunes(pc,pc2)/nblignes

---

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

Code:

---

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:

---

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

---

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

Code:

---

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

Mauvais cast.

Si on regarde ceci :

Code:

---

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

---

on a ceci

Code:

---

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

---

Ce qui revient implicitement a ceci :

Code:

---

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

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

Code:

---

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

---

ou ceci:

Code:

---

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

---

et avec cette nouvelle formule

Code:

---

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?

Un petit problème de parenthèse :

Code:

---

log2((double)

---

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

Code:

---

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

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.

voilà je décortique ma fonction:

Code:

---

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

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

Suite a MP avec cyrine, voici la formule :

Citation:

---

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 :

Citation:

---

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:

---

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

---

Cela marche-t'il ?

Super, ca marche :ccool:
Merci :)