`gamma(0)` returns the wrong value.

Open
#411 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Assessment

Difficulty
5/5
Estimated time
Over a week
Newbie friendliness
25/100
Issue type
Bug
Clarity
Mostly clear
Activity status
Stale
Tech stack
julia
Domain
backend

Research direction

Start by reviewing the gamma(0.0) and complex gamma cases shown in the issue, along with the linked opposite issue #163. Determine the intended behavior for real zero and signed complex zero, then identify the relevant implementation and tests in the repository. Done means the zero cases follow the agreed mathematical behavior and regression coverage exists.

Written by the indexing model from the issue text.

Description

This issue is the opposite issue of https://github.com/JuliaMath/SpecialFunctions.jl/issues/163. Currently gamma(0.0) returns Inf as does gamma(0.0+0.0im), and gamma(0.0-0.0im) This isn't the correct limit because

gamma(0. + 1e-16im)
0.6123233995736787 - 1.0000000000000034e16im

and

julia> gamma(0.0-1e-16im)
0.6123233995736787 + 1.0000000000000034e16im

which really makes it seem like gamma of complex 0s should return an appropriate complex infinity and gamma(0.0) should error.

Dominant language
Julia
Stars
384
Forks
115
Avg merge
2d 9h
Merged PRs (30d)
2

Getting set up

This project ships no dev container, Dockerfile or contributing guide, so setting up is up to you: start from its README, and see our first-contribution guide for the general steps.

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from JuliaMath/SpecialFunctions.jl

All issues in JuliaMath/SpecialFunctions.jl

Similar issues

More Julia issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.