Re: tolerance in binopdf.m

[Jordi Gutiérrez Hermoso]

>> On 9/21/2011 3:07 PM, Ben Abbott wrote:
>>> Thus, it looks to me as if Apple has a different implementation of
>>> lgamma (?).
>>>
>>> As this is not a bug in Octave, I'm inclined to add a tolerance
>>> for MaOS.
>>>
>>> However, I'm curious about what Apple did.
>>>
>>> I'm using Xcode 4.1. The sources are below.
>>>
>>>      http://www.opensource.apple.com/
>>>
>>> Does anyone have an idea of where to look to find the sources for
>>> lgamma?

> On Sep 21, 2011, at 11:01 AM, Marco atzeri wrote:
>>
>> Could be here ?
>>
>> http://opensource.apple.com/source/Libm/
>> http://opensource.apple.com/source/Libm/Libm-2026/Source/Intel/xmm_erfgamma.c

On 21 September 2011 10:25, Ben Abbott <address@hidden> wrote:
> Yes, that is it!

Who knows why, but their implementation does differ slightly from the
one in GNU's libc:

    
http://www.eglibc.org/cgi-bin/viewcvs.cgi/trunk/libc/sysdeps/ieee754/dbl-64/e_lgamma_r.c?rev=3747&view=auto

HTH,
- Jordi G. H.