1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
# roundf.m4 serial 12
dnl Copyright (C) 2007-2011 Free Software Foundation, Inc.
dnl This file is free software; the Free Software Foundation
dnl gives unlimited permission to copy and/or distribute it,
dnl with or without modifications, as long as this notice is preserved.
AC_DEFUN([gl_FUNC_ROUNDF],
[
m4_divert_text([DEFAULTS], [gl_roundf_required=plain])
AC_REQUIRE([gl_MATH_H_DEFAULTS])
dnl Persuade glibc <math.h> to declare roundf().
AC_REQUIRE([gl_USE_SYSTEM_EXTENSIONS])
AC_CHECK_DECLS([roundf], , , [#include <math.h>])
if test "$ac_cv_have_decl_roundf" = yes; then
gl_CHECK_MATH_LIB([ROUNDF_LIBM], [x = roundf (x);])
if test "$ROUNDF_LIBM" != missing; then
dnl Test whether roundf() produces correct results. On mingw, for
dnl x = 1/2 - 2^-25, the system's roundf() returns a wrong result.
AC_REQUIRE([AC_PROG_CC])
AC_REQUIRE([AC_CANONICAL_HOST]) dnl for cross-compiles
AC_CACHE_CHECK([whether roundf works], [gl_cv_func_roundf_works],
[
save_LIBS="$LIBS"
LIBS="$LIBS $ROUNDF_LIBM"
AC_RUN_IFELSE([AC_LANG_SOURCE([[
#include <float.h>
#include <math.h>
int main()
{
/* 2^FLT_MANT_DIG. */
static const float TWO_MANT_DIG =
/* Assume FLT_MANT_DIG <= 3 * 31.
Use the identity n = floor(n/3) + floor((n+1)/3) + floor((n+2)/3). */
(float) (1U << (FLT_MANT_DIG / 3))
* (float) (1U << ((FLT_MANT_DIG + 1) / 3))
* (float) (1U << ((FLT_MANT_DIG + 2) / 3));
volatile float x = 0.5f - 0.5f / TWO_MANT_DIG;
exit (x < 0.5f && roundf (x) != 0.0f);
}]])], [gl_cv_func_roundf_works=yes], [gl_cv_func_roundf_works=no],
[case "$host_os" in
mingw*) gl_cv_func_roundf_works="guessing no";;
*) gl_cv_func_roundf_works="guessing yes";;
esac
])
LIBS="$save_LIBS"
])
case "$gl_cv_func_roundf_works" in
*no) ROUNDF_LIBM=missing ;;
esac
fi
if test "$ROUNDF_LIBM" = missing; then
REPLACE_ROUNDF=1
fi
m4_ifdef([gl_FUNC_ROUNDF_IEEE], [
if test $gl_roundf_required = ieee && test $REPLACE_ROUNDF = 0; then
AC_CACHE_CHECK([whether roundf works according to ISO C 99 with IEC 60559],
[gl_cv_func_roundf_ieee],
[
save_LIBS="$LIBS"
LIBS="$LIBS $ROUNDF_LIBM"
AC_RUN_IFELSE(
[AC_LANG_SOURCE([[
#ifndef __NO_MATH_INLINES
# define __NO_MATH_INLINES 1 /* for glibc */
#endif
#include <math.h>
]gl_FLOAT_MINUS_ZERO_CODE[
]gl_FLOAT_SIGNBIT_CODE[
static float dummy (float f) { return 0; }
int main (int argc, char *argv[])
{
float (*my_roundf) (float) = argc ? roundf : dummy;
/* Test whether roundf (-0.0f) is -0.0f. */
if (signbitf (minus_zerof) && !signbitf (my_roundf (minus_zerof)))
return 1;
return 0;
}
]])],
[gl_cv_func_roundf_ieee=yes],
[gl_cv_func_roundf_ieee=no],
[gl_cv_func_roundf_ieee="guessing no"])
LIBS="$save_LIBS"
])
case "$gl_cv_func_roundf_ieee" in
*yes) ;;
*) REPLACE_ROUNDF=1 ;;
esac
fi
])
else
HAVE_DECL_ROUNDF=0
fi
if test $HAVE_DECL_ROUNDF = 0 || test $REPLACE_ROUNDF = 1; then
dnl Find libraries needed to link lib/roundf.c.
AC_CHECK_DECLS([ceilf, floorf], , , [#include <math.h>])
if test "$ac_cv_have_decl_floorf" = yes \
&& test "$ac_cv_have_decl_ceilf" = yes; then
gl_FUNC_FLOORF_LIBS
gl_FUNC_CEILF_LIBS
if test "$FLOORF_LIBM" != '?' && test "$CEILF_LIBM" != '?'; then
AC_DEFINE([HAVE_FLOORF_AND_CEILF], [1],
[Define if the both the floorf() and ceilf() functions exist.])
ROUNDF_LIBM=
dnl Append $FLOORF_LIBM to ROUNDF_LIBM, avoiding gratuitous duplicates.
case " $ROUNDF_LIBM " in
*" $FLOORF_LIBM "*) ;;
*) ROUNDF_LIBM="$ROUNDF_LIBM $FLOORF_LIBM" ;;
esac
dnl Append $CEILF_LIBM to ROUNDF_LIBM, avoiding gratuitous duplicates.
case " $ROUNDF_LIBM " in
*" $CEILF_LIBM "*) ;;
*) ROUNDF_LIBM="$ROUNDF_LIBM $CEILF_LIBM" ;;
esac
else
ROUNDF_LIBM=
fi
else
ROUNDF_LIBM=
fi
fi
AC_SUBST([ROUNDF_LIBM])
])
|
