From mboxrd@z Thu Jan  1 00:00:00 1970
X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on inbox.vuxu.org
X-Spam-Level: 
X-Spam-Status: No, score=-3.3 required=5.0 tests=MAILING_LIST_MULTI,
	RCVD_IN_DNSWL_MED,RCVD_IN_MSPIKE_H3,RCVD_IN_MSPIKE_WL autolearn=ham
	autolearn_force=no version=3.4.4
Received: (qmail 13762 invoked from network); 5 Feb 2021 07:18:18 -0000
Received: from mother.openwall.net (195.42.179.200)
  by inbox.vuxu.org with ESMTPUTF8; 5 Feb 2021 07:18:18 -0000
Received: (qmail 25724 invoked by uid 550); 5 Feb 2021 07:18:14 -0000
Mailing-List: contact musl-help@lists.openwall.com; run by ezmlm
Precedence: bulk
List-Post: <mailto:musl@lists.openwall.com>
List-Help: <mailto:musl-help@lists.openwall.com>
List-Unsubscribe: <mailto:musl-unsubscribe@lists.openwall.com>
List-Subscribe: <mailto:musl-subscribe@lists.openwall.com>
List-ID: <musl.lists.openwall.com>
Reply-To: musl@lists.openwall.com
Received: (qmail 25689 invoked from network); 5 Feb 2021 07:18:14 -0000
X-IronPort-AV: E=Sophos;i="5.81,154,1610406000"; 
   d="scan'208";a="372157746"
Date: Fri, 05 Feb 2021 08:18:02 +0100
Message-Id: <mwim77oudx.fsf@tomate.loria.fr>
From: Paul Zimmermann <Paul.Zimmermann@inria.fr>
To: musl@lists.openwall.com
Subject: [musl] issue with acoshf

       Hi,

while updating to my comparison of the accuracy of mathematical functions [1],
I have noticed an issue with acoshf in musl-1.2.2:

$ cat test_acosh_musl.c
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int
main (int argc, char *argv[])
{
  float x = -0x1.1e6ae8p+5;
  float y;
  y = acoshf (x);
  printf ("x=%a y=%a\n", x, y);
}

With gcc I get NaN as expected:

$ gcc -fno-builtin test_acosh_musl.c -lm
$ ./a.out
x=-0x1.1e6ae8p+5 y=-nan

With musl-1.2.2 I get -0x1.2f63acp+3:

$ gcc -fno-builtin test_acosh_musl.c $FILES
$ ./a.out
x=-0x1.1e6ae8p+5 y=-0x1.2f63acp+3

Please can someone confirm?

Best regards,
Paul

[1] https://members.loria.fr/PZimmermann/papers/accuracy.pdf