bug package hypergeo

[Stéphane Laurent]

Bonjour,

Le package hypergeo fournit une implémentation remarquable de la fonction hypergéométrique de Gauss.

Hélas j'ai découvert des cas où le résultat n'est pas correct. Par exemple

Code : Tout sélectionner

> x <- 30; hypergeo(x+6,11,x+12.5,0.7)
> [1] -34603008+0i

alors que le résultat théorique est un nombre positif.

Je l'ai signalé à l'auteur du package et je ne comprends pas totalement sa réponse: I believe that the software is correct but your example exposes a limitation of IEEE arithmetic.

Qu'est-ce que IEEE arithmetic ? Et y aurait-il une solution si le problème est bien dû à cela ?