From abd383584b16dd0fb1bbf40e4ece65ebe7b839ec Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Fri, 10 Sep 2021 20:56:22 +0000 Subject: Add narrowing square root functions This patch adds the narrowing square root functions from TS 18661-1 / TS 18661-3 / C2X to glibc's libm: fsqrt, fsqrtl, dsqrtl, f32sqrtf64, f32sqrtf32x, f32xsqrtf64 for all configurations; f32sqrtf64x, f32sqrtf128, f64sqrtf64x, f64sqrtf128, f32xsqrtf64x, f32xsqrtf128, f64xsqrtf128 for configurations with _Float64x and _Float128; __f32sqrtieee128 and __f64sqrtieee128 aliases in the powerpc64le case (for calls to fsqrtl and dsqrtl when long double is IEEE binary128). Corresponding tgmath.h macro support is also added. The changes are mostly similar to those for the other narrowing functions previously added, so the description of those generally applies to this patch as well. However, the not-actually-narrowing cases (where the two types involved in the function have the same floating-point format) are aliased to sqrt, sqrtl or sqrtf128 rather than needing a separately built not-actually-narrowing function such as was needed for add / sub / mul / div. Thus, there is no __nldbl_dsqrtl name for ldbl-opt because no such name was needed (whereas the other functions needed such a name since the only other name for that entry point was e.g. f32xaddf64, not reserved by TS 18661-1); the headers are made to arrange for sqrt to be called in that case instead. The DIAG_* calls in sysdeps/ieee754/soft-fp/s_dsqrtl.c are because they were observed to be needed in GCC 7 testing of riscv32-linux-gnu-rv32imac-ilp32. The other sysdeps/ieee754/soft-fp/ files added didn't need such DIAG_* in any configuration I tested with build-many-glibcs.py, but if they do turn out to be needed in more files with some other configuration / GCC version, they can always be added there. I reused the same test inputs in auto-libm-test-in as for non-narrowing sqrt rather than adding extra or separate inputs for narrowing sqrt. The tests in libm-test-narrow-sqrt.inc also follow those for non-narrowing sqrt. Tested as followed: natively with the full glibc testsuite for x86_64 (GCC 11, 7, 6) and x86 (GCC 11); with build-many-glibcs.py with GCC 11, 7 and 6; cross testing of math/ tests for powerpc64le, powerpc32 hard float, mips64 (all three ABIs, both hard and soft float). The different GCC versions are to cover the different cases in tgmath.h and tgmath.h tests properly (GCC 6 has _Float* only as typedefs in glibc headers, GCC 7 has proper _Float* support, GCC 8 adds __builtin_tgmath). --- NEWS | 7 + manual/arith.texi | 19 + math/Makefile | 6 +- math/Versions | 9 + math/auto-libm-test-out-narrow-sqrt | 4485 ++++++++++++++++++++ math/bits/mathcalls-narrow.h | 3 + math/gen-tgmath-tests.py | 5 +- math/libm-test-driver.c | 18 +- math/libm-test-narrow-sqrt.inc | 56 + math/math-narrow.h | 145 +- math/math.h | 6 + math/tgmath.h | 56 + math/w_sqrt_compat.c | 6 + math/w_sqrt_template.c | 8 + math/w_sqrtl_compat.c | 4 + sysdeps/generic/math-narrow-alias.h | 135 + sysdeps/generic/math-type-macros-double.h | 5 + sysdeps/generic/math-type-macros-float.h | 4 + sysdeps/generic/math-type-macros-float128.h | 5 + sysdeps/generic/math-type-macros-ldouble.h | 5 + sysdeps/generic/math-type-macros.h | 7 + sysdeps/i386/fpu/w_sqrt.c | 2 + sysdeps/i386/fpu/w_sqrt_compat.c | 2 + sysdeps/ia64/fpu/e_sqrt.S | 1 + sysdeps/ia64/fpu/libm-symbols.h | 1 + sysdeps/ieee754/dbl-64/s_f32xsqrtf64.c | 1 + sysdeps/ieee754/dbl-64/s_fsqrt.c | 34 + sysdeps/ieee754/float128/float128_private.h | 5 + sysdeps/ieee754/float128/s_f32sqrtf128.c | 6 + sysdeps/ieee754/float128/s_f64sqrtf128.c | 10 + sysdeps/ieee754/float128/s_f64xsqrtf128.c | 2 + sysdeps/ieee754/ldbl-128/s_dsqrtl.c | 37 + sysdeps/ieee754/ldbl-128/s_f64xsqrtf128.c | 38 + sysdeps/ieee754/ldbl-128/s_fsqrtl.c | 33 + sysdeps/ieee754/ldbl-128ibm-compat/Versions | 4 + sysdeps/ieee754/ldbl-128ibm/s_dsqrtl.c | 28 + sysdeps/ieee754/ldbl-128ibm/s_fsqrtl.c | 28 + sysdeps/ieee754/ldbl-96/s_dsqrtl.c | 33 + sysdeps/ieee754/ldbl-96/s_fsqrtl.c | 31 + sysdeps/ieee754/ldbl-opt/Makefile | 4 +- sysdeps/ieee754/ldbl-opt/nldbl-dsqrt.c | 28 + sysdeps/ieee754/ldbl-opt/nldbl-fsqrt.c | 28 + sysdeps/ieee754/soft-fp/s_dsqrtl.c | 67 + sysdeps/ieee754/soft-fp/s_fsqrt.c | 54 + sysdeps/ieee754/soft-fp/s_fsqrtl.c | 53 + sysdeps/mach/hurd/i386/libm.abilist | 13 + sysdeps/powerpc/fpu/libm-test-ulps | 3 + .../le/fpu/multiarch/float128-ifunc-macros.h | 2 + .../powerpc64/le/fpu/multiarch/float128-ifunc.h | 6 + sysdeps/unix/sysv/linux/aarch64/libm.abilist | 13 + sysdeps/unix/sysv/linux/alpha/libm.abilist | 13 + sysdeps/unix/sysv/linux/arc/libm.abilist | 6 + sysdeps/unix/sysv/linux/arm/be/libm.abilist | 6 + sysdeps/unix/sysv/linux/arm/le/libm.abilist | 6 + sysdeps/unix/sysv/linux/csky/libm.abilist | 6 + sysdeps/unix/sysv/linux/hppa/libm.abilist | 6 + sysdeps/unix/sysv/linux/i386/libm.abilist | 13 + sysdeps/unix/sysv/linux/ia64/libm.abilist | 13 + sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist | 6 + sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist | 6 + sysdeps/unix/sysv/linux/microblaze/be/libm.abilist | 6 + sysdeps/unix/sysv/linux/microblaze/le/libm.abilist | 6 + sysdeps/unix/sysv/linux/mips/mips32/libm.abilist | 6 + sysdeps/unix/sysv/linux/mips/mips64/libm.abilist | 13 + sysdeps/unix/sysv/linux/nios2/libm.abilist | 6 + .../sysv/linux/powerpc/powerpc32/fpu/libm.abilist | 6 + .../linux/powerpc/powerpc32/nofpu/libm.abilist | 6 + .../sysv/linux/powerpc/powerpc64/be/libm.abilist | 6 + .../sysv/linux/powerpc/powerpc64/le/libm.abilist | 15 + sysdeps/unix/sysv/linux/riscv/rv32/libm.abilist | 13 + sysdeps/unix/sysv/linux/riscv/rv64/libm.abilist | 13 + sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist | 13 + sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist | 13 + sysdeps/unix/sysv/linux/sh/be/libm.abilist | 6 + sysdeps/unix/sysv/linux/sh/le/libm.abilist | 6 + sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist | 13 + sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist | 13 + sysdeps/unix/sysv/linux/x86_64/64/libm.abilist | 13 + sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist | 13 + 79 files changed, 5736 insertions(+), 92 deletions(-) create mode 100644 math/auto-libm-test-out-narrow-sqrt create mode 100644 math/libm-test-narrow-sqrt.inc create mode 100644 sysdeps/generic/math-narrow-alias.h create mode 100644 sysdeps/ieee754/dbl-64/s_f32xsqrtf64.c create mode 100644 sysdeps/ieee754/dbl-64/s_fsqrt.c create mode 100644 sysdeps/ieee754/float128/s_f32sqrtf128.c create mode 100644 sysdeps/ieee754/float128/s_f64sqrtf128.c create mode 100644 sysdeps/ieee754/float128/s_f64xsqrtf128.c create mode 100644 sysdeps/ieee754/ldbl-128/s_dsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-128/s_f64xsqrtf128.c create mode 100644 sysdeps/ieee754/ldbl-128/s_fsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-128ibm/s_dsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-128ibm/s_fsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-96/s_dsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-96/s_fsqrtl.c create mode 100644 sysdeps/ieee754/ldbl-opt/nldbl-dsqrt.c create mode 100644 sysdeps/ieee754/ldbl-opt/nldbl-fsqrt.c create mode 100644 sysdeps/ieee754/soft-fp/s_dsqrtl.c create mode 100644 sysdeps/ieee754/soft-fp/s_fsqrt.c create mode 100644 sysdeps/ieee754/soft-fp/s_fsqrtl.c diff --git a/NEWS b/NEWS index 5b014fabbf..838381f562 100644 --- a/NEWS +++ b/NEWS @@ -19,6 +19,13 @@ Major new features: in various downstream distributions. The locale is not built into glibc, and must be installed. +* functions that round their results to a narrower type, and + corresponding macros, are added from TS 18661-1:2014, TS + 18661-3:2015 and draft ISO C2X: + + - fsqrt, fsqrtl, dsqrtl and corresponding fMsqrtfN, fMsqrtfNx, + fMxsqrtfN and fMxsqrtfNx functions. + Deprecated and removed features, and other changes affecting compatibility: [Add deprecations, removals and changes affecting compatibility here] diff --git a/manual/arith.texi b/manual/arith.texi index 75eaf67fe7..edf1fd71f8 100644 --- a/manual/arith.texi +++ b/manual/arith.texi @@ -2270,6 +2270,25 @@ function without any intermediate rounding to the type of the arguments. @end deftypefun +@deftypefun float fsqrt (double @var{x}) +@deftypefunx float fsqrtl (long double @var{x}) +@deftypefunx double dsqrtl (long double @var{x}) +@deftypefunx _FloatM fMsqrtfN (_Float@var{N} @var{x}) +@deftypefunx _FloatM fMsqrtfNx (_Float@var{N}x @var{x}) +@deftypefunx _FloatMx fMxsqrtfN (_Float@var{N} @var{x}) +@deftypefunx _FloatMx fMxsqrtfNx (_Float@var{N}x @var{x}) +@standards{TS 18661-1:2014, math.h} +@standardsx{fMsqrtfN, TS 18661-3:2015, math.h} +@standardsx{fMsqrtfNx, TS 18661-3:2015, math.h} +@standardsx{fMxsqrtfN, TS 18661-3:2015, math.h} +@standardsx{fMxsqrtfNx, TS 18661-3:2015, math.h} +@safety{@prelim{}@mtsafe{}@assafe{}@acsafe{}} +These functions, from TS 18661-1:2014 and TS 18661-3:2015, return the +square root of @var{x}, rounded once to the return type of the +function without any intermediate rounding to the type of the +arguments. +@end deftypefun + @node Complex Numbers @section Complex Numbers @pindex complex.h diff --git a/math/Makefile b/math/Makefile index ceb1eb2085..987e9841f9 100644 --- a/math/Makefile +++ b/math/Makefile @@ -91,7 +91,7 @@ libm-compat-calls = \ w_lgammaF_r_compat w_lgammaF_compat2 w_expF_compat \ w_lgamma_compatF k_standardF -libm-narrow-fns = add div mul sub +libm-narrow-fns = add div mul sqrt sub libm-narrow-types-basic = s_fF s_f32xFf64 libm-narrow-types-ldouble-yes = s_fFl s_dFl libm-narrow-types-float128-yes = s_f32Ff128 s_f64Ff128 s_f64xFf128 @@ -299,7 +299,7 @@ libm-test-funcs-noauto = canonicalize ceil cimag conj copysign cproj creal \ significand totalorder totalordermag trunc ufromfp \ ufromfpx compat_totalorder compat_totalordermag libm-test-funcs-compat = compat_totalorder compat_totalordermag -libm-test-funcs-narrow = add div mul sub +libm-test-funcs-narrow = add div mul sqrt sub libm-test-funcs-all = $(libm-test-funcs-auto) $(libm-test-funcs-noauto) libm-test-c-auto = $(foreach f,$(libm-test-funcs-auto),libm-test-$(f).c) libm-test-c-noauto = $(foreach f,$(libm-test-funcs-noauto),libm-test-$(f).c) @@ -664,6 +664,7 @@ CFLAGS-s_ctanh.c += -fno-builtin-ctanhl CFLAGS-s_dadd.c += -fno-builtin-daddl CFLAGS-s_ddiv.c += -fno-builtin-ddivl CFLAGS-s_dmul.c += -fno-builtin-dmull +CFLAGS-s_dsqrt.c += -fno-builtin-dsqrtl CFLAGS-s_dsub.c += -fno-builtin-dsubl CFLAGS-s_erf.c += -fno-builtin-erfl CFLAGS-s_erfc.c += -fno-builtin-erfcl @@ -687,6 +688,7 @@ CFLAGS-s_fmul.c += -fno-builtin-fmull CFLAGS-s_frexp.c += -fno-builtin-frexpl CFLAGS-s_fromfp.c += -fno-builtin-fromfpl CFLAGS-s_fromfpx.c += -fno-builtin-fromfpxl +CFLAGS-s_fsqrt.c += -fno-builtin-fsqrtl CFLAGS-s_fsub.c += -fno-builtin-fsubl CFLAGS-s_gamma.c += -fno-builtin-gammal CFLAGS-s_getpayload.c += -fno-builtin-getpayloadl diff --git a/math/Versions b/math/Versions index 5b3c321490..b121b454cd 100644 --- a/math/Versions +++ b/math/Versions @@ -593,4 +593,13 @@ libm { GLIBC_2.32 { exp10f; } + GLIBC_2.35 { + # Functions not involving _Float64x or _Float128, for all configurations. + fsqrt; fsqrtl; dsqrtl; + f32sqrtf32x; f32sqrtf64; f32xsqrtf64; + # Functions involving _Float64x or _Float128, for some configurations. + f32sqrtf64x; f32sqrtf128; + f32xsqrtf64x; f32xsqrtf128; f64sqrtf64x; f64sqrtf128; + f64xsqrtf128; + } } diff --git a/math/auto-libm-test-out-narrow-sqrt b/math/auto-libm-test-out-narrow-sqrt new file mode 100644 index 0000000000..5e995b0780 --- /dev/null +++ b/math/auto-libm-test-out-narrow-sqrt @@ -0,0 +1,4485 @@ +sqrt 0 += sqrt downward binary32:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest binary32:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero binary32:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward binary32:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt downward binary64:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest binary64:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero binary64:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward binary64:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt downward intel96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest intel96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero intel96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward intel96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt downward m68k96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest m68k96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero m68k96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward m68k96:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt downward binary128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest binary128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero binary128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward binary128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt downward ibm128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt tonearest ibm128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt towardzero ibm128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : += sqrt upward ibm128:arg_fmt(0,0,0,0) 0x0p+0 : 0x0p+0 : +sqrt -0 += sqrt downward binary32:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest binary32:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero binary32:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward binary32:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt downward binary64:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest binary64:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero binary64:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward binary64:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt downward intel96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest intel96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero intel96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward intel96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt downward m68k96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest m68k96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero m68k96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward m68k96:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt downward binary128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest binary128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero binary128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward binary128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt downward ibm128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt tonearest ibm128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt towardzero ibm128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : += sqrt upward ibm128:arg_fmt(0,0,0,0) -0x0p+0 : -0x0p+0 : +sqrt 2209 += sqrt downward binary32:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest binary32:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero binary32:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward binary32:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt downward binary64:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest binary64:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero binary64:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward binary64:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt downward intel96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest intel96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero intel96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward intel96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt downward m68k96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest m68k96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero m68k96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward m68k96:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt downward binary128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest binary128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero binary128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward binary128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt downward ibm128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt tonearest ibm128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt towardzero ibm128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : += sqrt upward ibm128:arg_fmt(11,1,0,12) 0x8.a1p+8 : 0x2.fp+4 : +sqrt 4 += sqrt downward binary32:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest binary32:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero binary32:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward binary32:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt downward binary64:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest binary64:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero binary64:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward binary64:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt downward intel96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest intel96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero intel96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward intel96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt downward m68k96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest m68k96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero m68k96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward m68k96:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt downward binary128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest binary128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero binary128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward binary128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt downward ibm128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt tonearest ibm128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt towardzero ibm128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : += sqrt upward ibm128:arg_fmt(2,1,2,1) 0x4p+0 : 0x2p+0 : +sqrt 2 += sqrt downward binary32:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e6p+0 : inexact += sqrt tonearest binary32:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e6p+0 : inexact += sqrt towardzero binary32:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e6p+0 : inexact += sqrt upward binary32:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e8p+0 : inexact += sqrt downward binary64:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bccp+0 : inexact += sqrt tonearest binary64:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcdp+0 : inexact += sqrt towardzero binary64:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bccp+0 : inexact += sqrt upward binary64:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcdp+0 : inexact += sqrt downward intel96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt tonearest intel96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt towardzero intel96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt upward intel96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc90ap+0 : inexact += sqrt downward m68k96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt tonearest m68k96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt towardzero m68k96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908p+0 : inexact += sqrt upward m68k96:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc90ap+0 : inexact += sqrt downward binary128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea95p+0 : inexact += sqrt tonearest binary128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea95p+0 : inexact += sqrt towardzero binary128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea95p+0 : inexact += sqrt upward binary128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea96p+0 : inexact += sqrt downward ibm128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea8p+0 : inexact += sqrt tonearest ibm128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea8p+0 : inexact += sqrt towardzero ibm128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ea8p+0 : inexact += sqrt upward ibm128:arg_fmt(1,1,1,1) 0x2p+0 : 0x1.6a09e667f3bcc908b2fb1366ebp+0 : inexact +sqrt 0.25 += sqrt downward binary32:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest binary32:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero binary32:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward binary32:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt downward binary64:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest binary64:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero binary64:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward binary64:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt downward intel96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest intel96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero intel96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward intel96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt downward m68k96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest m68k96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero m68k96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward m68k96:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt downward binary128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest binary128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero binary128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward binary128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt downward ibm128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt tonearest ibm128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt towardzero ibm128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : += sqrt upward ibm128:arg_fmt(-2,1,-2,1) 0x4p-4 : 0x8p-4 : +sqrt 6642.25 += sqrt downward binary32:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest binary32:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero binary32:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward binary32:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt downward binary64:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest binary64:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero binary64:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward binary64:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt downward intel96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest intel96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero intel96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward intel96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt downward m68k96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest m68k96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero m68k96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward m68k96:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt downward binary128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest binary128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero binary128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward binary128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt downward ibm128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt tonearest ibm128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt towardzero ibm128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : += sqrt upward ibm128:arg_fmt(12,2,-2,15) 0x1.9f24p+12 : 0x5.18p+4 : +sqrt 15190.5625 += sqrt downward binary32:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest binary32:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero binary32:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward binary32:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt downward binary64:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest binary64:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero binary64:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward binary64:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt downward intel96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest intel96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero intel96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward intel96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt downward m68k96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest m68k96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero m68k96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward m68k96:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt downward binary128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest binary128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero binary128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward binary128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt downward ibm128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt tonearest ibm128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt towardzero ibm128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : += sqrt upward ibm128:arg_fmt(13,3,-4,18) 0x3.b569p+12 : 0x7.b4p+4 : +sqrt 0.75 += sqrt downward binary32:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d7p-4 : inexact += sqrt tonearest binary32:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d7p-4 : inexact += sqrt towardzero binary32:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d7p-4 : inexact += sqrt upward binary32:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d8p-4 : inexact += sqrt downward binary64:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c2655p-4 : inexact += sqrt tonearest binary64:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c2655p-4 : inexact += sqrt towardzero binary64:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c2655p-4 : inexact += sqrt upward binary64:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c26558p-4 : inexact += sqrt downward intel96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539dp-4 : inexact += sqrt tonearest intel96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539ep-4 : inexact += sqrt towardzero intel96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539dp-4 : inexact += sqrt upward intel96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539ep-4 : inexact += sqrt downward m68k96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539dp-4 : inexact += sqrt tonearest m68k96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539ep-4 : inexact += sqrt towardzero m68k96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539dp-4 : inexact += sqrt upward m68k96:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539ep-4 : inexact += sqrt downward binary128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c5cp-4 : inexact += sqrt tonearest binary128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c5cp-4 : inexact += sqrt towardzero binary128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c5cp-4 : inexact += sqrt upward binary128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c5c8p-4 : inexact += sqrt downward ibm128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c4p-4 : inexact += sqrt tonearest ibm128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c4p-4 : inexact += sqrt towardzero ibm128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c4p-4 : inexact += sqrt upward ibm128:arg_fmt(-1,2,-2,2) 0xcp-4 : 0xd.db3d742c265539d92ba16b83c8p-4 : inexact +sqrt 0x1.fffffffffffffp+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0x1p+512 : inexact += sqrt downward intel96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact += sqrt downward binary128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbfffffffffffff78p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbfffffffffffff8p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbfffffffffffff78p+508 : inexact += sqrt upward binary128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbfffffffffffff8p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffffffffffffcp+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffbffffffffffffcp+508 : inexact += sqrt upward ibm128:arg_fmt(1023,53,971,53) 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffcp+508 : inexact +sqrt 0x1.ffffffffffffbp+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffe8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffe8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffe8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.ffffffffffffp+508 : inexact += sqrt downward intel96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffecp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffecp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffecp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffecp+508 : inexact += sqrt downward binary128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff378p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff38p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff378p+508 : inexact += sqrt upward binary128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff38p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffffp+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff4p+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffffp+508 : inexact += sqrt upward ibm128:arg_fmt(1023,50,971,53) 0xf.fffffffffffd8p+1020 : 0xf.fffffffffffebffffffffffff4p+508 : inexact +sqrt 0x1.ffffffffffff7p+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffd8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffd8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffd8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffep+508 : inexact += sqrt downward intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdcp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdcp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdcp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdcp+508 : inexact += sqrt downward binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd778p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd78p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd778p+508 : inexact += sqrt upward binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd78p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd4p+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd8p+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd4p+508 : inexact += sqrt upward ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffffb8p+1020 : 0xf.fffffffffffdbfffffffffffd8p+508 : inexact +sqrt 0x1.ffffffffffff3p+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffc8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffc8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffc8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffdp+508 : inexact += sqrt downward intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffccp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffccp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffccp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffccp+508 : inexact += sqrt downward binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffab78p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffab8p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffab78p+508 : inexact += sqrt upward binary128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffab8p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffa8p+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffacp+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffa8p+508 : inexact += sqrt upward ibm128:arg_fmt(1023,49,971,53) 0xf.fffffffffff98p+1020 : 0xf.fffffffffffcbfffffffffffacp+508 : inexact +sqrt 0x1.fffffffffffefp+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffb8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffb8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffb8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffcp+508 : inexact += sqrt downward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbcp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbcp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbcp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbcp+508 : inexact += sqrt downward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6f78p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6f8p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6f78p+508 : inexact += sqrt upward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6f8p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6cp+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff7p+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff6cp+508 : inexact += sqrt upward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff78p+1020 : 0xf.fffffffffffbbfffffffffff7p+508 : inexact +sqrt 0x1.fffffffffffebp+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffa8p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffa8p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffa8p+508 : inexact += sqrt upward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffbp+508 : inexact += sqrt downward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffacp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffacp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffacp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffacp+508 : inexact += sqrt downward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff2378p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff238p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff2378p+508 : inexact += sqrt upward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff238p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff2p+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff24p+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff2p+508 : inexact += sqrt upward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff58p+1020 : 0xf.fffffffffffabfffffffffff24p+508 : inexact +sqrt 0x1.fffffffffffe7p+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt tonearest binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffd8p+60 : inexact += sqrt upward binary64:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward intel96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt tonearest m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt towardzero m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffp+60 : inexact += sqrt upward m68k96:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffep+60 : inexact += sqrt downward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt tonearest binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5ff8p+60 : inexact += sqrt upward binary128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt tonearest ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt towardzero ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff5cp+60 : inexact += sqrt upward ibm128:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffff7fffffdfffffeffffff6p+60 : inexact += sqrt downward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt tonearest binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : plus_infty : inexact overflow errno-erange += sqrt towardzero binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffp+124 : inexact overflow errno-erange-ok += sqrt upward binary32:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : plus_infty : inexact overflow errno-erange += sqrt downward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff98p+508 : inexact += sqrt tonearest binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff98p+508 : inexact += sqrt towardzero binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff98p+508 : inexact += sqrt upward binary64:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffffap+508 : inexact += sqrt downward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffp+508 : inexact += sqrt tonearest intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9cp+508 : inexact += sqrt towardzero intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffp+508 : inexact += sqrt upward intel96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9cp+508 : inexact += sqrt downward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffp+508 : inexact += sqrt tonearest m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9cp+508 : inexact += sqrt towardzero m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffp+508 : inexact += sqrt upward m68k96:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9cp+508 : inexact += sqrt downward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec778p+508 : inexact += sqrt tonearest binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec78p+508 : inexact += sqrt towardzero binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec778p+508 : inexact += sqrt upward binary128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec78p+508 : inexact += sqrt downward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec4p+508 : inexact += sqrt tonearest ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec8p+508 : inexact += sqrt towardzero ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec4p+508 : inexact += sqrt upward ibm128:arg_fmt(1023,48,971,53) 0xf.fffffffffff38p+1020 : 0xf.fffffffffff9bffffffffffec8p+508 : inexact +sqrt 0x1.fffffffffffe3p+1023 += sqrt downward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt tonearest binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt towardzero binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0xf.fffffp+60 : inexact += sqrt upward binary32:arg_fmt(127,24,104,24) 0xf.fffffp+124 : 0x1p+64 : inexact += sqrt downw