Comment puis-je contourner une erreur d'arrondi qui provoque une boucle infinie dans Perl's Statistics :: Descriptive?

Question

J'utilise la bibliothèque Statistics::Descriptive en Perl pour calculer les distributions de fréquence et faire face à un problème d'erreur d'arrondi en virgule flottante. Je passe deux valeurs, 0,205 et 0,205, (tirées d'autres nombres et sprintées vers celles-ci) au module stats et je lui demande de calculer la distribution de fréquence mais il se bloque dans une boucle infinie.

Parcourant avec un débogueur je peux voir que ça fait:

my $interval = $self->{sample_range}/$partitions; 

my $iter = $self->{min}; 

while (($iter += $interval) < $self->{max}) { 

    $bins{$iter} = 0; 

    push @k, $iter; ##Keep the "keys" unstringified 

} 

$ self-> sample_range (La plage est max-min) est de retour 2.77555756156289e-17 plutôt que 0 comme j'attends . Cela signifie que la boucle ((min + = plage) < max)) entre dans une boucle infinie (à toutes fins utiles).

DB < 8> print $ self -> {max};
0.205
DB < 9> print $ self -> {min};
0.205
DB < 10> print $ self -> {max} - $ self -> {min};
2.77555756156289e-17

Cela ressemble à un problème d'arrondi. Je ne sais pas comment résoudre ça de mon côté, et je ne suis pas sûr que l'édition de la bibliothèque soit une bonne idée. Je cherche des suggestions d'une solution de contournement ou d'une alternative.

Cheers, Neil

Answer 1

Je suis le Statistics :: Descriptive maintainer. En raison de sa nature numérique, de nombreux problèmes d'arrondis ont été signalés. Je crois que celui-ci a été corrigé dans une version ultérieure à celle que vous utilisiez que j'ai publiée récemment, en utilisant la multiplication pour les divisions au lieu de + =.

Veuillez utiliser the most up-to-date version du CPAN, et cela devrait être mieux.

Answer 2

Pas exactement un problème d'arrondi; vous pouvez voir les valeurs plus précises avec quelque chose comme

printf("%.18g %.18g", $self->{max}, $self->{min}); 

Me semble qu'il ya une faille dans le module où il assume la plage d'échantillon peut être divisé en partitions $ morceaux; parce que le point flottant n'a pas une précision infinie, ce n'est pas toujours possible. Dans votre cas, les valeurs min et max sont exactement des valeurs représentables adjacentes, il ne peut donc pas y avoir plus d'une partition. Je ne sais pas exactement ce que le module utilise pour les partitions, donc je ne suis pas sûr de l'impact que cela pourrait avoir. Un autre problème possible dans le module est qu'il utilise des nombres comme des clés de hachage, les stringifie implicitement ce qui arrondit légèrement la valeur.

Vous avez sans doute un certain succès dans le blanchiment de vos données par stringization avant de le nourrir au module:

$data = 0+"$data"; 

Cela permettra au moins faire en sorte que deux numéros (avec la précision d'impression par défaut) apparaissent égale sont en fait égal.

Answer 3

Cela ne devrait pas provoquer une boucle infinie. Ce qui ferait que cette boucle soit infinie serait $self->{sample_range}/$partitions est 0.