Discussion :

[me_myself]

Bonjour,

Question toute bête sans doute, mais j'ai par exemple:

double x,y,z;
z = x - y;

comment gérer l'underflow sur z ?

D'avance merci,
bv

(j'utilise MS VisualStudio 2005)

[souviron34]

comment gérer l'underflow sur z ?

Que veux-tu dire ?

[me_myself]

J'ai une DLL écrite en C qui est interfacée à Excel.
Je transmets donc des 'doubles' vers la DLL où des calculs sont effectués puis renvoyés vers la feuille de calcul.
Il se trouve que ces calculs génèrent parfois un underflow.


Merci d'avance
bv

[souviron34]

je repete : que veux-tu dire par "underflow" ????????????????????

[me_myself]

oops sorry...

une valeur trop proche de zéro, l'exact opposé de l'overflow, nombre trop grand pour la capacité d'un 'double'.

[Mat.M]

J'ai une DLL écrite en C qui est interfacée à Excel.

bv

"l'underflow" depuis que je fais de l'informatique c'est la première fois que je vosi ça de ma vie.
Ton problème c'est un problème de format et taille de données différents entre Excel/Visual Basic et le C.
Aller sur www.codeproject.com ou dans le MSDN il ya une table de différence entre Visual Basic et C/c++

[me_myself]

rien de bien extraordinaire pourtant.

je travaille sur un algorithme mathématique d'optimisation. On recherche donc une convergence entre une valeur théorique et une valeur réelle et dans certains cas, elle est très faible.
La plupart des fonctions comme sqrt par exemple de la librairie math, traitent les exceptions classiques y compris l'underflow, et je voudrais faire de même avec une simple soustraction.

[souviron34]

bon on commence qa voir plus clair..

Mais l'expression "underflow" est incomprise....

Maintenant pour ton probleme :

je signale que, dans les fichiers d'entetes de la biblotheque standard, tu as le fichier "limits.h", qui definit les valeurs max, min et de detection de chaque type de variable.

Tu as donc un DBL_EPSILON (ou EPSILON_DBL, je neme souviens plus) qui donne la plus petite variation mesurable entre 2 doubles.

Vu que c'est de l'ordre de 10^-13 (ou -17), m'etonnerait que pour des finances ce soit ton ordre de grandeur.....

Donc je ferais simplement une multiplication par 10000, ou quelque chose comme ca, si la valeur est < 1.

Sinon, si tu veux tester par rapport a zero, tu fais

Code :

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

if ( (a - b) == DBL_EPSILON )

ce qui equivaut a

Code :

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

if ( (a - b) == 0 )

Ou bien, si tu veux tester par rapport a une petite valeur :

Code :

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

1
2
3
4
 
#define ACCURACY 0.000001
 
if ( fabs(a-b) <= ACCURACY )

[me_myself]

Merci.

Essentiellement, ce qui m'interesse est de capturer l'erreur et de ramener alors la différence des 2 nombres à zéro. L'erreur est surtout visible dans Excel. Elle peut bien sûr être traitée en VBA, mais je cherche à ce que le traitement soit dans la DLL en C.

[souviron34]

donne un exemple, parce que la on nage un peu

[corentin59]

Code :

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

1
2
3
4
 
#define ACCURACY 0.000001
 
if ( fabs(a-b) <= ACCURACY )

Je pense que c'est la bonne approche. De toute manière, si tu travailles sur la convergence d'une variable, il faut obligatoirement que tu définisses un seuil en-dessous duquel tu considères qu'il y a convergence. A toi de choisir ce seuil plus grand que la valeur limite.

[diogene]

Je pense que c'est la bonne approche. De toute manière, si tu travailles sur la convergence d'une variable, il faut obligatoirement que tu définisses un seuil en-dessous duquel tu considères qu'il y a convergence. A toi de choisir ce seuil plus grand que la valeur limite.

Je crois qu'il est plus logique de comparer l'erreur relative plutôt que l'erreur absolue qui n'a pas grand sens compte tenu de la représentation des flottants

Code :

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

1
2
3
#define ACCURACY 0.000001
....
if ( fabs(a-b) <= ACCURACY * (fabs(a)+fabs(b)))

[souviron34]

Je crois qu'il est plus logique de comparer l'erreur relative plutôt que l'erreur absolue qui n'a pas grand sens compte tenu de la représentation des flottants

Code :

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

1
2
3
#define ACCURACY 0.000001
....
if ( fabs(a-b) <= ACCURACY * (fabs(a)+fabs(b)))

tout depend.. Dans la plupart des cas que j'ai eu a traiter, l'erreur absolue etait la regle :

• si Delta temperature < 0.5 degres
  
  
• si Delta euro < 0.01
  
  
• etc etc...

il n'y a que dans certains cas que tu veux un traitement relatif...

Mais par exemple pour un algo de convergence de resolution d'une formule de maths, si tu n'es pas en maths pures mais en physique ou chimie, tu SAIS que ta precision sera au maxi de 10^-3, -4, ou -6...

et ca fait des operations en moins (dont des divisions) pour du temps reel....

[corrector]

double x,y,z;
z = x - y;

comment gérer l'underflow sur z ?

Gérer comment? Quelle gestion souhaites-tu? Quel est le contexte?

Si tu veux dire détecter l'underflow, il faut préciser ce que tu entends par underflow :

1. quand le résultat mathématique est trop petit, la soustraction donne 0 (implémentation non-IEEE)
2. quand le résultat mathématique est trop petit, résultat de la soustraction perd en précision relative, sa mantisse est raccourcie (implémentation IEEE)

Le test a effectuer, dans chaque cas :

1. x != y && z == 0
2. z != 0 && fabs(z) < DBL_MIN

[corrector]

Tu as donc un DBL_EPSILON (ou EPSILON_DBL, je neme souviens plus) qui donne la plus petite variation mesurable entre 2 doubles.

Tu es sûr?

Donc je ferais simplement une multiplication par 10000, ou quelque chose comme ca, si la valeur est < 1.

Pourquoi faire cela?

Sinon, si tu veux tester par rapport a zero, tu fais

Code :

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

if ( (a - b) == DBL_EPSILON )

ce qui equivaut a

Code :

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

if ( (a - b) == 0 )

Ah bon? DBL_EPSILON = 0?

[souviron34]

Tu es sûr?

oui

La norme C99 précise que ce doit être dans le fichier <float.h>, et avant, c'était assez souvent dans le fichier <limits.h> ou <maths.h>.

Tu as les liens ici :

http://www.developpez.net/forums/sho...5&postcount=21

Pourquoi faire cela?

Parce que un nombre flottant multiplié par 10000 par exemple, et représentant des euros, te donne en eniter la valeur du 10 millième d'euro !!! et que je veux bien qu'on tente d'être précis, mais à moins d'extra, pour ce que tu souhaites faire, ça me semble nettement suffisant..

Ah bon? DBL_EPSILON = 0?

Non..

Mais

smallest number x such that 1.0 + x != 1.0",