commons-issues mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From Sébastien Brisard (JIRA) <>
Subject [jira] [Commented] (MATH-738) Incomplete beta function I(x, a, b) is inaccurate for large values of a and/or b
Date Wed, 13 Mar 2013 06:22:13 GMT


Sébastien Brisard commented on MATH-738:

I don't think these improvements are sufficient (the issue is really with the incomplete beta
function, while I started with the beta function). The accuracy of the incomplete beta function
should first be assessed (for example with the tools I have provided). In any case, this issue
cannot be solved until 4.0, since there is a deprecated funtion which cannot be removed before

Improving the incomplete beta function is a lot of work. I have provided the references in
the javadoc (mainly the NSWC FORTRAN library). As I have already indicated, I no longer have
time to do it myself (I'm sorry I forgot to remove myself as an assignee for this ticket).
If someone is willing to take over, I'm happy to provide some help and references, but I just
can't find the time to do the job. Hope you will understand.
> Incomplete beta function I(x, a, b) is inaccurate for large values of a and/or b
> --------------------------------------------------------------------------------
>                 Key: MATH-738
>                 URL:
>             Project: Commons Math
>          Issue Type: Bug
>    Affects Versions: 3.0
>            Reporter: Sébastien Brisard
>              Labels: special-functions
>             Fix For: 3.2, 4.0
> This was first reported in MATH-718. The result of the current implementation of the
incomplete beta function I(x, a, b) is inaccurate when a and/or b are large-ish. 
> I've skimmed through [slatec|], GSL, [Boost|]
as well as NR. At first sight, neither uses the same method to compute this function. I think
[TOMS-708|] is probably the best option.

This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrators
For more information on JIRA, see:

View raw message