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
|
/* $OpenBSD: modf_test.c,v 1.2 2023/08/13 06:57:04 miod Exp $ */
#include <assert.h>
#include <math.h>
/* Test for bug introduced in 4.4BSD modf() on sparc */
/* Public domain, 2014, Tobias Ulmer <tobiasu@tmux.org> */
#define BIGFLOAT (5e15) /* Number large enough to trigger the "big" case */
void
modf_sparc(void)
{
double f, i;
f = modf(BIGFLOAT, &i);
assert(i == BIGFLOAT);
assert(f == 0.0);
/* Repeat, maybe we were lucky */
f = modf(BIGFLOAT, &i);
assert(i == BIGFLOAT);
assert(f == 0.0);
/* With negative number, for good measure */
f = modf(-BIGFLOAT, &i);
assert(i == -BIGFLOAT);
assert(f == 0.0);
}
/* Test for modf() behaviour on Inf and Nan */
/* Written by Willemijn Coene. Public domain */
void
modf_infnan(void)
{
double f, i;
f = modf(__builtin_inf(), &i);
assert(isinf(i));
assert(signbit(i) == 0);
assert(f == 0.0);
f = modf(-__builtin_inf(), &i);
assert(isinf(i));
assert(signbit(i) != 0);
assert(f == -0.0);
f = modf(NAN, &i);
assert(isnan(i));
assert(signbit(i) == 0);
assert(isnan(f));
assert(signbit(f) == 0);
f = modf(-NAN, &i);
assert(isnan(i));
assert(signbit(i) != 0);
assert(isnan(f));
assert(signbit(f) != 0);
}
int
main(void)
{
modf_sparc();
modf_infnan();
}
|
