Attachment 7383 Details for Bug 16519 – test program

test program

test6.c (text/x-c), 1.17 KB, created by Paul Zimmermann on 2014-02-01 09:13:06 UTC

(hide)

Description:

Filename:

MIME Type:

Creator: Paul Zimmermann

Created: 2014-02-01 09:13:06 UTC

Size: 1.17 KB

patch

obsolete

>#include <stdio.h>
>#include <math.h>
>#include <mpfr.h>
>#include <gnu/libc-version.h>
>#include <fenv.h>
>
>float
>ulp (float x, float y)
>{
>  float z = x - y, t, twoy;
>
>  twoy = y + y;
>  if (twoy != twoy + twoy) /* twoy is finite */
>    t = nextafterf (y, twoy);
>  else /* twoy is +Inf or -Inf */
>    t = nextafter (y, (float) 0);
>  t = y - t;
>  if (t < (float) 0)
>    t = -t;
>  /* now t = ulp(y) */
>  return z / t;
>}
>
>int
>main ()
>{
>  float x = 0x5.000000p-148;
>  float y, z;
>  mpfr_t t;
>  int inex, runderflow;
>
>  printf("GNU libc version: %s\n", gnu_get_libc_version ());
>  printf("GNU libc release: %s\n", gnu_get_libc_release ());
>  fesetround (FE_TONEAREST);
>  feclearexcept (FE_ALL_EXCEPT);
>  y = sinhf (x);
>  runderflow = fetestexcept (FE_UNDERFLOW);
>  mpfr_set_emin (-125 - 24 + 1);
>  mpfr_set_emax (128);
>  mpfr_init2 (t, 24);
>  mpfr_set_flt (t, x, MPFR_RNDN);
>  inex = mpfr_sinh (t, t, MPFR_RNDN);
>  inex = mpfr_subnormalize (t, inex, MPFR_RNDN);
>  z = mpfr_get_flt (t, MPFR_RNDD);
>  printf ("x=%a (%.9e)\n", x, x);
>  printf ("libc: sinhf(x)=%a (%.9e) underflow=%d\n", y, y, runderflow);
>  printf ("mpfr: sinhf(x)=%a (%.9e)\n", z, z);
>  printf ("ulp difference: %f\n", ulp (y, z));
>  return 0;
>}

Actions: View

Attachments on bug 16519: 7383