Discussion :

[xav181]

bonjour,

voila j'y perd mon latin, j'effectue des calculs, généralement tout ce passe bien mais sur certaines valeurs, j'obtiens des valeurs NaN sans raison;

lorsque je fais:
double x = (wPos[i]+epsilon)/(wNeg[i]+epsilon);
cerr<<"x "<<x<<endl;
j'obtiens: x 0.0138675
mais lorsque je fais:
double val = log(x);
cerr<<val<<endl;
j'obtiens : nan.
Mais si je fais: cerr<<log(double((wPos[i]+epsilon)/(wNeg[i]+epsilon)))<<endl;
j'obtiens :
-4.27821

Bon voila je m'arrache les cheveux, si quelqu'un peut m'aider.
Merci

[NiamorH]

t'es sur que tu fais pas :

Code :

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

cer << log( val ) << endl;

à la place de :

Code :

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

cer << val << endl;

Moi je ne reproduis pas.

[xav181]

Non, c'est bien ce que je fais; j'ai copié le code. En fait, pour quasiment toutes les valeurs je n'ai pas de problème et pour cette valeur particulière j'ai cette erreur, j'ai l'impression que c'est peut-être un bug de la fonction log sous mac, je code xcode sous leopard.

[JolyLoic]

Tu pourrais nous donner un programme complet qui reproduise ça ? Avec entre autre les valeurs exactes des wPos, wNeg et epsilon ?

[xav181]

En fait, je pense pas que ce soit une erreur de log puisque l'instruction suivante passe:
cerr<<log(double((wPos[i]+epsilon)/(wNeg[i]+epsilon)))<<endl;
mais doit y avoir un truc qui m'échappe. Voici exactement tout ce que je fais:

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
const double epsilon = 0.0001;
valarray<double> wPos;
valarray<double> wNeg;
/.../je remplis les valarray
for(int i = 0; i < wPos.size();++i){
 
      cerr<<" wp "<<wPos[i]<<" wn "<<wNeg[i]<<endl;
      double wp = wPos[i];
      double wn = wNeg[i];
      double test = (wp+epsilon)/(wn+epsilon);
      double x = (wPos[i]+epsilon)/(wNeg[i]+epsilon);
      double val_0 = log(x);
      double val_1 = log(test);
      cerr<<"test "<<test<<endl;
      cerr<<"x "<<x<<endl;
      cerr<<"val 0"<<val_0<<endl;
      cerr<<"val 1"<<val_1<<endl;
     cerr<<log(double((wPos[i]+epsilon)/(wNeg[i]+epsilon)))<<endl;
}

Voila j'ai bug pour la valeur i = 0;
ou j'obtiens le resultat suivant:
wp 0 wn 0.00711111
test 0.0138675
x 0.0138675
val 0nan
val 1nan
-4.27821