* lib/math.in.h (setpayloadl): New declaration.
* lib/setpayloadl.c: New file.
* m4/math_h.m4 (gl_MATH_H): Test whether setpayloadl is declared.
(gl_MATH_H_REQUIRE_DEFAULTS): Initialize GNULIB_SETPAYLOADL.
(gl_MATH_H_DEFAULTS): Initialize HAVE_SETPAYLOADL.
* modules/math (Makefile.am): Substitute GNULIB_SETPAYLOADL,
HAVE_SETPAYLOADL.
* modules/setpayloadl: New file.
* doc/posix-functions/setpayloadl.texi: Mention the new module.
+2024-04-16 Bruno Haible <bruno@clisp.org>
+
+ setpayloadl: New module.
+ * lib/math.in.h (setpayloadl): New declaration.
+ * lib/setpayloadl.c: New file.
+ * m4/math_h.m4 (gl_MATH_H): Test whether setpayloadl is declared.
+ (gl_MATH_H_REQUIRE_DEFAULTS): Initialize GNULIB_SETPAYLOADL.
+ (gl_MATH_H_DEFAULTS): Initialize HAVE_SETPAYLOADL.
+ * modules/math (Makefile.am): Substitute GNULIB_SETPAYLOADL,
+ HAVE_SETPAYLOADL.
+ * modules/setpayloadl: New file.
+ * doc/posix-functions/setpayloadl.texi: Mention the new module.
+
2024-04-16 Bruno Haible <bruno@clisp.org>
setpayloadf: Add tests.
@url{https://www.gnu.org/software/libc/manual/html_node/FP-Bit-Twiddling.html}.
@end ifnotinfo
-Gnulib module: ---
+Gnulib module: setpayloadl
Portability problems fixed by Gnulib:
@itemize
+@item
+This function is missing on all non-glibc platforms:
+glibc 2.24, macOS 11.1, FreeBSD 14.0, NetBSD 10.0, OpenBSD 6.7, Minix 3.1.8, AIX 7.1, HP-UX 11.31, IRIX 6.5, Solaris 11.4, Cygwin 2.9, mingw, MSVC 14, Android 9.0.
@end itemize
Portability problems not fixed by Gnulib:
@itemize
-@item
-This function is missing on all non-glibc platforms:
-glibc 2.24, macOS 11.1, FreeBSD 14.0, NetBSD 10.0, OpenBSD 6.7, Minix 3.1.8, AIX 7.1, HP-UX 11.31, IRIX 6.5, Solaris 11.4, Cygwin 2.9, mingw, MSVC 14, Android 9.0.
@end itemize
# endif
#endif
+#if @GNULIB_SETPAYLOADL@
+# if !@HAVE_SETPAYLOADL@
+_GL_FUNCDECL_SYS (setpayloadl, int, (long double *, long double));
+# endif
+_GL_CXXALIAS_SYS (setpayloadl, int, (long double *, long double));
+_GL_CXXALIASWARN (setpayloadl);
+#elif defined GNULIB_POSIXCHECK
+# undef setpayloadl
+# if HAVE_RAW_DECL_SETPAYLOADL
+_GL_WARN_ON_USE (setpayloadl, "setpayloadl is unportable - "
+ "use gnulib module setpayloadl for portability");
+# endif
+#endif
+
#if @GNULIB_TOTALORDERF@
# if @REPLACE_TOTALORDERF@
--- /dev/null
+/* Construct a quiet NaN 'long double' with a given payload.
+ Copyright 2024 Free Software Foundation, Inc.
+
+ This file is free software: you can redistribute it and/or modify
+ it under the terms of the GNU Lesser General Public License as
+ published by the Free Software Foundation, either version 3 of the
+ License, or (at your option) any later version.
+
+ This file is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU Lesser General Public License for more details.
+
+ You should have received a copy of the GNU Lesser General Public License
+ along with this program. If not, see <https://www.gnu.org/licenses/>. */
+
+/* Written by Bruno Haible. */
+
+#include <config.h>
+
+/* Specification. */
+#include <math.h>
+
+#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE
+
+int
+setpayloadl (long double *result, long double payload)
+{
+ return setpayload ((double *) result, payload);
+}
+
+#else
+
+# include <float.h>
+# include <stdint.h>
+
+# include "signed-nan.h"
+# include "snan.h"
+
+/* 2^(LDBL_MANT_DIG-1). */
+# define TWO_LDBL_MANT_DIG \
+ ((long double) (1U << ((LDBL_MANT_DIG - 1) / 4)) \
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 1) / 4)) \
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 2) / 4)) \
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 3) / 4)))
+
+int
+setpayloadl (long double *result, long double payload)
+{
+# if (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 106 || LDBL_MANT_DIG == 113) \
+ && defined LDBL_EXPBIT0_WORD && defined LDBL_EXPBIT0_BIT
+ if (
+# if defined __hppa || defined __mips__ || defined __sh__
+ /* A zero payload is not allowed, because that would denote Infinity.
+ Cf. snan.h. */
+ payload > 0.0L
+# else
+ payload >= 0.0L
+# endif
+# if LDBL_MANT_DIG == 64 /* on i386, x86_64, ia64, m68k */
+ && payload < 4611686018427387904.0L /* exp2l (LDBL_MANT_DIG - 2) */
+ && payload == (long double) (int64_t) payload
+# endif
+# if LDBL_MANT_DIG == 106 /* on powerpc, powerpc64, powerpc64le */
+ && payload < 2251799813685248.0L /* exp2l (DBL_MANT_DIG - 2) */
+ && payload == (long double) (int64_t) payload
+# endif
+# if LDBL_MANT_DIG == 113 /* on alpha, arm64, loongarch64, mips64, riscv64, s390x, sparc64 */
+ && payload < 2596148429267413814265248164610048.0L /* exp2l (LDBL_MANT_DIG - 2) */
+ && payload == +(+(payload + TWO_LDBL_MANT_DIG) - TWO_LDBL_MANT_DIG)
+# endif
+ )
+ {
+ memory_long_double x;
+ x.value = positive_NaNl ();
+
+# if LDBL_MANT_DIG == 64 /* on i386, x86_64, ia64, m68k */
+ uint64_t pl = (int64_t) payload;
+# if LDBL_EXPBIT0_WORD == 2 && LDBL_EXPBIT0_BIT == 0 /* on i386, x86_64, ia64 */
+ x.word[1] = (x.word[1] & 0xC0000000U) | (uint32_t) (pl >> 32);
+ x.word[0] = (uint32_t) pl;
+# elif LDBL_EXPBIT0_WORD == 0 && LDBL_EXPBIT0_BIT == 16 /* on m68k */
+ x.word[0] = (x.word[0] & 0xFFFFC000U) | (uint32_t) (pl >> 48)
+ x.word[1] = (uint32_t) (pl >> 16);
+ x.word[2] = (uint32_t) (pl << 16);
+# else
+# error "Please port gnulib setpayloadl.c to your platform!"
+# endif
+# endif
+# if LDBL_MANT_DIG == 106 /* on powerpc, powerpc64, powerpc64le */
+ uint64_t pl = (int64_t) payload;
+# if LDBL_EXPBIT0_BIT == 20
+ uint32_t pl_hi = (uint32_t) (pl >> 32);
+ uint32_t pl_lo = (uint32_t) pl;
+ x.word[LDBL_EXPBIT0_WORD] = (x.word[LDBL_EXPBIT0_WORD] & 0xFFF80000U)
+ | pl_hi;
+ x.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 1 : -1)] = pl_lo;
+# else
+# error "Please port gnulib setpayloadl.c to your platform!"
+# endif
+# endif
+# if LDBL_MANT_DIG == 113 /* on alpha, arm64, loongarch64, mips64, riscv64, s390x, sparc64 */
+# if LDBL_EXPBIT0_BIT == 16
+ memory_long_double pl;
+ pl.value = payload + TWO_LDBL_MANT_DIG;
+ x.word[LDBL_EXPBIT0_WORD] = (x.word[LDBL_EXPBIT0_WORD] & 0xFFFFC000U)
+ | (pl.word[LDBL_EXPBIT0_WORD] & 0x00003FFFU);
+ x.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 1 : -1)] =
+ pl.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 1 : -1)];
+ x.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 2 : -2)] =
+ pl.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 2 : -2)];
+ x.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 3 : -3)] =
+ pl.word[LDBL_EXPBIT0_WORD + (LDBL_EXPBIT0_WORD == 0 ? 3 : -3)];
+# else
+# error "Please port gnulib setpayloadl.c to your platform!"
+# endif
+# endif
+
+ *result = x.value;
+ return 0;
+ }
+ else
+ {
+ *result = 0.0L;
+ return -1;
+ }
+# else
+# error "Please port gnulib setpayloadl.c to your platform!"
+# endif
+}
+
+#endif
# math_h.m4
-# serial 128
+# serial 129
dnl Copyright (C) 2007-2024 Free Software Foundation, Inc.
dnl This file is free software; the Free Software Foundation
dnl gives unlimited permission to copy and/or distribute it,
modf modff modfl powf
remainder remainderf remainderl
rint rintf rintl round roundf roundl
- setpayload setpayloadf
+ setpayload setpayloadf setpayloadl
sinf sinl sinhf sqrtf sqrtl
tanf tanl tanhf totalorder totalorderf totalorderl trunc truncf truncl])
])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_ROUNDL])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SETPAYLOAD])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SETPAYLOADF])
+ gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SETPAYLOADL])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SIGNBIT])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SINF])
gl_MODULE_INDICATOR_INIT_VARIABLE([GNULIB_SINL])
HAVE_RINTL=1; AC_SUBST([HAVE_RINTL])
HAVE_SETPAYLOAD=1; AC_SUBST([HAVE_SETPAYLOAD])
HAVE_SETPAYLOADF=1; AC_SUBST([HAVE_SETPAYLOADF])
+ HAVE_SETPAYLOADL=1; AC_SUBST([HAVE_SETPAYLOADL])
HAVE_SINF=1; AC_SUBST([HAVE_SINF])
HAVE_SINL=1; AC_SUBST([HAVE_SINL])
HAVE_SINHF=1; AC_SUBST([HAVE_SINHF])
-e 's/@''GNULIB_ROUNDL''@/$(GNULIB_ROUNDL)/g' \
-e 's/@''GNULIB_SETPAYLOAD''@/$(GNULIB_SETPAYLOAD)/g' \
-e 's/@''GNULIB_SETPAYLOADF''@/$(GNULIB_SETPAYLOADF)/g' \
+ -e 's/@''GNULIB_SETPAYLOADL''@/$(GNULIB_SETPAYLOADL)/g' \
-e 's/@''GNULIB_SIGNBIT''@/$(GNULIB_SIGNBIT)/g' \
-e 's/@''GNULIB_SINF''@/$(GNULIB_SINF)/g' \
-e 's/@''GNULIB_SINL''@/$(GNULIB_SINL)/g' \
-e 's|@''HAVE_RINTL''@|$(HAVE_RINTL)|g' \
-e 's|@''HAVE_SETPAYLOAD''@|$(HAVE_SETPAYLOAD)|g' \
-e 's|@''HAVE_SETPAYLOADF''@|$(HAVE_SETPAYLOADF)|g' \
+ -e 's|@''HAVE_SETPAYLOADL''@|$(HAVE_SETPAYLOADL)|g' \
-e 's|@''HAVE_SINF''@|$(HAVE_SINF)|g' \
-e 's|@''HAVE_SINL''@|$(HAVE_SINL)|g' \
-e 's|@''HAVE_SINHF''@|$(HAVE_SINHF)|g' \
--- /dev/null
+Description:
+setpayloadl function: construct a quiet NaN with a given payload
+
+Files:
+lib/setpayloadl.c
+m4/mathfunc.m4
+m4/setpayload.m4
+m4/exponentl.m4
+
+Depends-on:
+math
+extensions
+setpayload [test $HAVE_SETPAYLOADL = 0 && test $HAVE_SAME_LONG_DOUBLE_AS_DOUBLE = 1]
+float [test $HAVE_SETPAYLOADL = 0 && test $HAVE_SAME_LONG_DOUBLE_AS_DOUBLE = 0]
+stdint [test $HAVE_SETPAYLOADL = 0 && test $HAVE_SAME_LONG_DOUBLE_AS_DOUBLE = 0]
+signed-nan [test $HAVE_SETPAYLOADL = 0 && test $HAVE_SAME_LONG_DOUBLE_AS_DOUBLE = 0]
+snan [test $HAVE_SETPAYLOADL = 0 && test $HAVE_SAME_LONG_DOUBLE_AS_DOUBLE = 0]
+
+configure.ac:
+gl_FUNC_SETPAYLOADL
+gl_CONDITIONAL([GL_COND_OBJ_SETPAYLOADL], [test $HAVE_SETPAYLOADL = 0])
+gl_MATH_MODULE_INDICATOR([setpayloadl])
+
+Makefile.am:
+if GL_COND_OBJ_SETPAYLOADL
+lib_SOURCES += setpayloadl.c
+endif
+
+Include:
+<math.h>
+
+Link:
+$(SETPAYLOADL_LIBM)
+
+License:
+LGPL
+
+Maintainer:
+all