JuliaDiff / JuliaDiff/ForwardDiff.jl
SpecialFunctions.gamma_inc goes into infinite recursion
Open
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 1k
- Forks
- 160
- PR merge metrics
- No merged PRs in 30d
Description
The function SpecialFunctions.gamma_inc goes into infinite recursion when used with ForwardDiff.Duals. MWE:
using ForwardDiff, SpecialFunctions
julia> f = x -> gamma_inc(x, 1.0, 1)
julia> f(ForwardDiff.Dual(1.0, 1.0))
ERROR: StackOverflowError:
Stacktrace:
[1] gamma_inc(::ForwardDiff.Dual{Nothing,Float64,1}, ::Float64, ::Int64) (repeats 80000 times)
julia> VERSION
v"1.2.0"
julia> Pkg.status("ForwardDiff")
- ForwardDiff v0.10.6
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reproducing the supplied Julia MWE using ForwardDiff.Dual and SpecialFunctions.gamma_inc. Trace the gamma_inc dispatch involved in the stack overflow, then add or update a regression test showing that the call terminates correctly and preserves the expected result.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- tooling
- Issue type
- Bug
- Difficulty
- 3/5
- Estimated time
- 1-2 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 35/100