From: C Howland <cc1964t@gmail.com>
To: newlib@sourceware.org
Subject: Re: Fw: issue with tgamma
Date: Fri, 5 Jan 2024 12:34:43 -0500 [thread overview]
Message-ID: <CANk6obTZSyDx7BU9tdTjfWuCLdT2+jSHANME9XD8ZScM9T5MGA@mail.gmail.com> (raw)
In-Reply-To: <BN2P110MB154427B901AD06FB47A831F59A66A@BN2P110MB1544.NAMP110.PROD.OUTLOOK.COM>
[-- Attachment #1: Type: text/plain, Size: 1626 bytes --]
> *From:* Paul Zimmermann <Paul.Zimmermann@inria.fr>
> *Sent:* Friday, January 5, 2024 5:49 AM
> *To:* newlib@sourceware.org <newlib@sourceware.org>
> *Subject:* issue with tgamma
>
> Hi,
>
> I discovered the following issue with tgamma: for
> x=-0x1.13a5ccd87c9bbp+1008
> it yields +Inf instead of NaN:
>
> Checking tgamma with newlib-4.4.0.20231231
> tgamma 0 -1 -0x1.13a5ccd87c9bbp+1008 [-nan] [inf] inf inf
> libm gives inf
> mpfr gives -nan
>
> This issue was present in earlier versions of Newlib, but we only found it
> recently since we improved our testing framework.
>
> Paul
>
> This strictly is OK, as the standards allow either a range or domain
error for negative integers, and domain errors may return an
implementation-defined value. It most likely has been this way since the
1993 date in the source file, which includes
* 5. Special Cases
* lgamma(2+s) ~ s*(1-Euler) for tiny s
* lgamma(1)=lgamma(2)=0
* lgamma(x) ~ -log(x) for tiny x
* lgamma(0) = lgamma(inf) = inf
* lgamma(-integer) = +-inf
(The tgamma function ends up using the lgamma function.)
POSIX does have an implied preference of sorts for the NaN return, so
changing it could possibly be considered a small improvement even if not
necessary--although it is strange that they define the result for argument
0 differently, calling for INF there, seems inconsistent.
(It's somewhat interesting that the MPFR return value you show is
-NaN. I suppose they chose that to help reflect it being due to a negative
argument.)
Craig
next prev parent reply other threads:[~2024-01-05 17:34 UTC|newest]
Thread overview: 3+ messages / expand[flat|nested] mbox.gz Atom feed top
2024-01-05 10:49 Paul Zimmermann
[not found] ` <BN2P110MB154427B901AD06FB47A831F59A66A@BN2P110MB1544.NAMP110.PROD.OUTLOOK.COM>
2024-01-05 17:34 ` C Howland [this message]
2024-01-08 14:45 ` Fw: " Paul Zimmermann
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=CANk6obTZSyDx7BU9tdTjfWuCLdT2+jSHANME9XD8ZScM9T5MGA@mail.gmail.com \
--to=cc1964t@gmail.com \
--cc=newlib@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).