mailing list of musl libc
 help / color / mirror / code / Atom feed
* [musl] issue with exp10l
@ 2021-01-07 11:17 Paul Zimmermann
  2021-01-07 19:49 ` Rich Felker
  0 siblings, 1 reply; 4+ messages in thread
From: Paul Zimmermann @ 2021-01-07 11:17 UTC (permalink / raw)
  To: musl

       Hi,

I am extending my comparison of the accuracy of several mathematical libraries
to the "double extended precision" (long double on x86_64).

First I notice that Musl does not provide j0, j1, y0, and y1 for the long
double format. Do you confirm?

Then I got a segmentation fault using exp10l with NaN input with a non-zero
payload.

$ cat test_exp10.c
#define _GNU_SOURCE

typedef union { __uint128_t n; long double x; } union_t;

#include <stdio.h>
#include <math.h>
int main()
{
  union_t u;
  u.n = 16383UL;
  u.n = u.n << 64;
  u.n = u.n | 629329181547216221UL;
  /* u.n = 302213637488765131341149 */
  long double x = u.x;
  printf ("x=%La\n", x);
  fflush (stdout);
  long double y;
  y = exp10l (x);
  printf ("y=%La\n", y);
  fflush (stdout);
  return 0;
}

With glibc this works fine:

$ gcc -fno-builtin test_exp10.c -lm
$ ./a.out
x=nan
y=-nan

With Musl 1.2.1 I get:

$ ./a.out
x=nan
Segmentation fault

According to gdb, the issue is in pow10l:

Program received signal SIGSEGV, Segmentation fault.
0x000055555555d10e in pow10l ()
(gdb) where
#0  0x000055555555d10e in pow10l ()
#1  0x0000000080000000 in ?? ()
#2  0x0000000000003fff in ?? ()
#3  0x0000000000000000 in ?? ()

Best regards,
Paul


^ permalink raw reply	[flat|nested] 4+ messages in thread

end of thread, other threads:[~2021-01-08  9:07 UTC | newest]

Thread overview: 4+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2021-01-07 11:17 [musl] issue with exp10l Paul Zimmermann
2021-01-07 19:49 ` Rich Felker
2021-01-07 20:02   ` Rich Felker
2021-01-08  9:06     ` Paul Zimmermann

Code repositories for project(s) associated with this public inbox

	https://git.vuxu.org/mirror/musl/

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 NNTP newsgroup(s).