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
|
#include <math.h>
#include <float.h>
#include <assert.h>
int main()
{
double d;
#ifndef _MSC_VER
// first check constants
assert(isnormal(FLT_MAX));
assert(isinf(HUGE_VAL));
assert(isinf(HUGE_VALF));
// assert(isinf(HUGE_VALL));
assert(isinf(INFINITY));
assert(isnan(NAN));
// check +
assert(isinf(INFINITY+INFINITY));
assert(isnan(-INFINITY+INFINITY));
assert(INFINITY+INFINITY>0);
assert(isnan(NAN+d));
assert(isnan(NAN+INFINITY));
// check -
assert(isnan(INFINITY-INFINITY));
assert(isinf(-INFINITY-INFINITY));
assert(-INFINITY-INFINITY<0);
assert(isnan(NAN-d));
assert(isnan(NAN-INFINITY));
// check *
assert(isinf(INFINITY*INFINITY));
assert(isinf(-INFINITY*INFINITY));
assert(INFINITY*INFINITY>0);
assert(-INFINITY*INFINITY<0);
assert(isnan(NAN*d));
assert(isnan(NAN*INFINITY));
assert(isnan(INFINITY*0));
assert(signbit(1.0*-0.0));
assert(!signbit(1.0*0.0));
// check /
assert(isnan(INFINITY/INFINITY));
assert(isnan(-INFINITY/INFINITY));
assert(isinf(INFINITY/0));
assert(0.0/INFINITY==0);
assert(1.0/INFINITY==0);
assert(signbit(-1.0/INFINITY));
assert(signbit(1.0/-INFINITY));
assert(INFINITY/-2<0);
assert(isinf(1.0/0.0));
assert(isinf(INFINITY/2));
assert(isnan(0.0/0.0));
assert(isnan(NAN/d));
assert(isnan(NAN/INFINITY));
assert(signbit(-0.0/1));
#endif
}
|
