GNU Scientific Library - Bugs: bug #45925, Incomplete Gamma Functions flipped?
You are not allowed to post comments on this tracker with your current authentication level.
bug #45925: Incomplete Gamma Functions flipped?
Submitter: | Patrick Alken <psa> | ||
Submitted: | Fri 11 Sep 2015 10:07:55 AM UTC | ||
Category: | Runtime error | Severity: | 3 - Normal |
Operating System: | Status: | None | |
Assigned to: | None | Open/Closed: | Open |
Release: |
Thu 19 Nov 2015 10:22:06 PM UTC, comment #1: |
Patrick Alken <psa>![]() |
Fri 11 Sep 2015 10:07:55 AM UTC, original submission:
from vsakharuk =at= hbk =dot= com
|
Patrick Alken <psa>![]() |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
I've investigated bug #45925 "Incomplete Gamma Functions flipped?" and I believe there is in fact no problem with gsl. The bug author compares the (upper) incomplete gamma function with the function P(a,x) and correctly notes that the values disagree. However the comparison should actually be done with the function Q(a,x) rather than P(a,x), and here the values agree. The definitions for the upper incomplete gamma function and Q(a,x) are contained in the gsl reference document under the section on "Incomplete Gamma Functions". I recommend that no code changes be made and that this bug be closed.
-Bill Maier