[
https://issues.apache.org/jira/browse/MATH-849?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=13439271#comment-13439271
]
Sébastien Brisard edited comment on MATH-849 at 8/22/12 3:44 PM:
-----------------------------------------------------------------
First shot. The relative accuracy is now down to 3 ulps. Unit tests relating to
MATH-735 now fail; will investigate.
Please review, especially the references to NSWC library. Thanks!
[EDIT] The files are attached as both patch {{.patch}} and drop-in replacement
of classes {{.zip}}.
was (Author: celestin):
First shot. The relative accuracy is now down to 3 ulps. Unit tests
relating to MATH-735 now fail; will investigate.
Please review, especially the references to NSWC library. Thanks!
> Precision improvement in class Gamma
> ------------------------------------
>
> Key: MATH-849
> URL: https://issues.apache.org/jira/browse/MATH-849
> Project: Commons Math
> Issue Type: Improvement
> Reporter: Sébastien Brisard
> Assignee: Sébastien Brisard
> Priority: Minor
> Labels: accuracy, special-functions
> Fix For: 3.1
>
> Attachments: MATH-849.patch, MATH-849.zip
>
>
> New unit tests (see below) show that (on the set of values tested) the
> current accuracy of {{Gamma.logGamma(double)}} is about 130 ulps.
> The [NSWC Library of Mathematical
> Functions|http://www.dtic.mil/docs/citations/ADA476840] proposes an
> implementation of this function with much higher accuracy: on the same set of
> values, the error goes down to about 3 ulps.
> This library is available
> [here|http://www.ualberta.ca/CNS/RESEARCH/Software/NumericalNSWC/site.html],
> and is license-free [see this
> thread|http://markmail.org/thread/fg4y4dfdnfty5hov]. It is proposed to port
> this implementation in Commons-Math. This in turn could lead to a resolution
> of MATH-738.
--
This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrators:
https://issues.apache.org/jira/secure/ContactAdministrators!default.jspa
For more information on JIRA, see: http://www.atlassian.com/software/jira
