haskell / haskell/math-functions
invIncompleteBeta fails for some (exotic) inputs
Open
bug
- Dominant language
- Haskell
- Stars
- 45
- Forks
- 32
- PR merge metrics
- No merged PRs in 30d
Description
```
*Numeric.SpecFunctions.Internal> invIncompleteBeta 9 2 1e-300
*** Exception: Numeric.SpecFunctions.incompletBeta_: x out of [0,1] range. p=9.0 q=2.0 x=NaN
```
It fails in Halley step when derivative of incomplete beta underflows and becomes zero so no next approximation could be computed, It happens for large `a` and small `x`.
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.