Re: [Rd] R-devel on FreeBSD: new C99 functions don't build

From: Rainer Hurling <rhurlin_at_gwdg.de>
Date: Tue, 15 May 2012 21:45:13 +0200

On 15.05.2012 20:49 (UTC+1), Murray Stokely wrote:
> On Tue, May 15, 2012 at 10:05 AM, Rainer Hurling<rhurlin@gwdg.de> wrote:
>> About April 25th, there had been some changes within R-devel's
>> src/nmath/pnbeta.c (and probably some other relevant places) and now

>> building R-devel on FreeBSD 10.0-CURRENT (amd64) with gcc-4.6.4 and
>> math/R-devel (selfmade forked port from math/R) fails like this:
>
>> It seems, that at least one new C99 function (log1pl) is introduced in
>> R-devel, see
>>
>> src/nmath/pnbeta.c:l95
>> return (double) (log_p ? log1pl(-ans) : (1 - ans));
>
> AFAIK, Bruce Evans is not happy with the numerical accuracy of other
> open-source implementations of log1pl and so has blocked their
> inclusion in FreeBSD pending work on a better implementation.
>
> Can you put a conditional FreeBSD check here and use log1p instead of
> log1pl instead as a workaround?
>
> I can admire the insistence on correctness from the FreeBSD libm
> maintainers for their technical purity, but it can be a bit of a pain
> for things like this.
>
> - Murray

I read about this discussion and in principle I concur with your opinion. As a scientist I tend to expect greatest possible correctness from a statistical routine, especially when it uses long double format.

As a quick and dirty workaround I applied the following patch:

It builds and installs fine now and I hope there are no side effects ...

Thank you for the quick answer and your advice, Rainer



R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue 15 May 2012 - 19:53:17 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 18 May 2012 - 15:16:32 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive