[ https://issues.apache.org/jira/browse/NUMBERS52?page=com.atlassian.jira.plugin.system.issuetabpanels:alltabpanel
]
Gilles updated NUMBERS52:

Description:
This was first reported in MATH718. The result of the current implementation of the incomplete
beta function I(x, a, b) is inaccurate when a and/or b are largeish.
I've skimmed through [slatechttp://www.netlib.org/slatec/fnlib/betai.f], GSL, [Boosthttp://www.boost.org/doc/libs/1_38_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/sf_beta/ibeta_function.html]
as well as NR. At first sight, neither uses the same method to compute this function. I think
[TOMS708http://www.netlib.org/toms/708] is probably the best option.
_Issue moved from MATH project on January 27, 2018 (concerned implementation was moved to
module {{commonsnumbersgamma}} of "Commons Numbers")._
was:
This was first reported in MATH718. The result of the current implementation of the incomplete
beta function I(x, a, b) is inaccurate when a and/or b are largeish.
I've skimmed through [slatechttp://www.netlib.org/slatec/fnlib/betai.f], GSL, [Boosthttp://www.boost.org/doc/libs/1_38_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/sf_beta/ibeta_function.html]
as well as NR. At first sight, neither uses the same method to compute this function. I think
[TOMS708http://www.netlib.org/toms/708] is probably the best option.
> Incomplete beta function I(x, a, b) is inaccurate for large values of a and/or b
> 
>
> Key: NUMBERS52
> URL: https://issues.apache.org/jira/browse/NUMBERS52
> Project: Commons Numbers
> Issue Type: Bug
> Reporter: Sébastien Brisard
> Priority: Major
> Labels: specialfunctions
>
> This was first reported in MATH718. The result of the current implementation of the
incomplete beta function I(x, a, b) is inaccurate when a and/or b are largeish.
> I've skimmed through [slatechttp://www.netlib.org/slatec/fnlib/betai.f], GSL, [Boosthttp://www.boost.org/doc/libs/1_38_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/sf_beta/ibeta_function.html]
as well as NR. At first sight, neither uses the same method to compute this function. I think
[TOMS708http://www.netlib.org/toms/708] is probably the best option.
> _Issue moved from MATH project on January 27, 2018 (concerned implementation was moved
to module {{commonsnumbersgamma}} of "Commons Numbers")._

This message was sent by Atlassian JIRA
(v7.6.3#76005)
