Discussion :

[L. Sierra]

Salut à tous,

Lorsque j'utilise des algorithmes linéaires de résolution de systèmes d'équation linéaires (genre Gradient Conjugué), il arrive que les valeurs tendent vers l'infini .

Je voudrais donc pouvoir "détecter" si une valeur réelle de type double ou float est "valide" (c'est à dire pas #IND ou #NAN) pour lancer des exceptions dans le cas contraire. Existe-t-il un moyen ?

Merci d'avance !

L.S.

[yan]

il y as isnan(double x) et finite(double x)

[Mat007]

Salut,

Pour NaN (ou IND qui est je crois un 'quiet' NaN chez Microsoft) une solution simple (et je pense portable) est :

Code :

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

1
2
3
 
    if( value != value )
        ;

Fascinant hein ?

Sinon pour les nombres infinis il suffit de faire :

Code :

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

1
2
3
4
5
 
    if( value == std::numeric_limits< float >::infinity() )
        ;
    if( value == - std::numeric_limits< float >::infinity() )
        ;

Si c'est vraiment des cas d'infinis obtenus par division par zéro je pense qu'il s'agit plutôt des seconds.

edit : typo malheureuse, il fallait lire 'quiet NaN' et non 'quite NaN'

MAT.

[L. Sierra]

Ok sympa, je vais essayer tout ça .

Merci A+

[nicroman]

Je voudrai pas dire de bétise, mais il me semble bien qu'une exception est levée lors d'une division par 0.
Ce qui se traduit par un 'crash' du programme

Me gourre-je ?

[Mat007]

une exception est levée lors d'une division par 0

C'est en C que ça plante fourbement, non ?
En tous cas le code suivant est tout à fait valide et ne lève pas d'exception :

Code :

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

1
2
3
4
5
6
 
BOOST_AUTO_TEST_CASE( DivideByZeroEqualsInfinity )
{
    const float zero = 0.f;
    BOOST_CHECK_EQUAL( std::numeric_limits< float >::infinity(), 1.f / zero );
}

Bon là ça produit un warning sous vc80 et si on met 1.f / 0.f c'est carrément une erreur

MAT.