JuliaMath / JuliaMath/HypergeometricFunctions.jl

pFq falsely evaluates to NaN

Open
#73 2 comments 0 reactions 0 assignees View on GitHub
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.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.