JuliaMath / JuliaMath/HypergeometricFunctions.jl
pFq falsely evaluates to NaN
- Dominant language
- Julia
- Stars
- 76
- Forks
- 12
- PR merge metrics
- No merged PRs in 30d
Description
Thanks for this great package. Has been really useful, however I noticed a potential bug when evaluating a specific value:
I got
```
julia> pFq((-4, -3, 151), (2, -153), -1.0)
NaN
```
while nearby values are computed correctly, such as
```
julia> pFq((-4, -3, 151), (2, -154), -1.0)
13.708301502419149
```
Mathematica would compute values
```
In[1]:= HypergeometricPFQ[{-4, -3, 151}, {2, -153}, -1.0]
Out[1]= 13.8431
In[2]:= HypergeometricPFQ[{-4, -3, 151}, {2, -154}, -1.0]
Out[2]= 13.7083
```
so matching the second value and providing a numerical value instead of NaN for the first case.
Increasing precision does not help
```
julia> pFq((-4, -3, 151), (2, -153), big(-1.0))
NaN
```
Julia version v"1.8.5", SpecialFunctions v2.3.0, HypergeometricFunctions v0.3.23.
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.