[musl] issue with acoshf

[Paul Zimmermann <Paul.Zimmermann@inria.fr>]

       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