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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
|
/* Test of modf() function.
$Id: modf-01.c,v 1.1 2007/02/05 21:35:58 dmix Exp $
*/
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "progmem.h"
union lofl_u {
long lo;
float fl;
};
volatile union lofl_u vf = { 1 };
double vn = { 1 };
PROGMEM const struct { /* Table of test cases. */
union lofl_u x; /* argument */
union lofl_u f; /* fraction part: modf(x,*pn) */
union lofl_u n; /* integral part: *pn */
} t[] = {
/* zero */
{ { 0x00000000 }, { 0x00000000 }, { .fl = 0 } },
{ { 0x80000000 }, { 0x80000000 }, { .fl = -0.0 } },
/* subnormal */
{ { 0x00000001 }, { 0x00000001 }, { .fl = 0 } },
{ { 0x00000100 }, { 0x00000100 }, { .fl = 0 } },
{ { 0x00010000 }, { 0x00010000 }, { .fl = 0 } },
{ { 0x007fffff }, { 0x007fffff }, { .fl = 0 } },
{ { 0x80000001 }, { 0x80000001 }, { .fl = -0.0 } },
{ { 0x80000100 }, { 0x80000100 }, { .fl = -0.0 } },
{ { 0x80010000 }, { 0x80010000 }, { .fl = -0.0 } },
{ { 0x807fffff }, { 0x807fffff }, { .fl = -0.0 } },
/* integral normal */
{ { .fl = 0x000001 }, { .fl = 0 }, { .fl = 0x000001 } },
{ { .fl = 0x000002 }, { .fl = 0 }, { .fl = 0x000002 } },
{ { .fl = 0x000004 }, { .fl = 0 }, { .fl = 0x000004 } },
{ { .fl = 0x000008 }, { .fl = 0 }, { .fl = 0x000008 } },
{ { .fl = 0x000010 }, { .fl = 0 }, { .fl = 0x000010 } },
{ { .fl = 0x000020 }, { .fl = 0 }, { .fl = 0x000020 } },
{ { .fl = 0x000040 }, { .fl = 0 }, { .fl = 0x000040 } },
{ { .fl = 0x000080 }, { .fl = 0 }, { .fl = 0x000080 } },
{ { .fl = 0x000100 }, { .fl = 0 }, { .fl = 0x000100 } },
{ { .fl = 0x000200 }, { .fl = 0 }, { .fl = 0x000200 } },
{ { .fl = 0x000400 }, { .fl = 0 }, { .fl = 0x000400 } },
{ { .fl = 0x000800 }, { .fl = 0 }, { .fl = 0x000800 } },
{ { .fl = 0x001000 }, { .fl = 0 }, { .fl = 0x001000 } },
{ { .fl = 0x002000 }, { .fl = 0 }, { .fl = 0x002000 } },
{ { .fl = 0x004000 }, { .fl = 0 }, { .fl = 0x004000 } },
{ { .fl = 0x008000 }, { .fl = 0 }, { .fl = 0x008000 } },
{ { .fl = 0x010000 }, { .fl = 0 }, { .fl = 0x010000 } },
{ { .fl = 0x020000 }, { .fl = 0 }, { .fl = 0x020000 } },
{ { .fl = 0x040000 }, { .fl = 0 }, { .fl = 0x040000 } },
{ { .fl = 0x080000 }, { .fl = 0 }, { .fl = 0x080000 } },
{ { .fl = 0x100000 }, { .fl = 0 }, { .fl = 0x100000 } },
{ { .fl = 0x200000 }, { .fl = 0 }, { .fl = 0x200000 } },
{ { .fl = 0x400000 }, { .fl = 0 }, { .fl = 0x400000 } },
{ { .fl = 0x7fffff }, { .fl = 0 }, { .fl = 0x7fffff } },
{ { .fl = -0x000001L }, { .fl = -0.0 }, { .fl = -0x000001L } },
{ { .fl = -0x000002L }, { .fl = -0.0 }, { .fl = -0x000002L } },
{ { .fl = -0x000080L }, { .fl = -0.0 }, { .fl = -0x000080L } },
{ { .fl = -0x000100L }, { .fl = -0.0 }, { .fl = -0x000100L } },
{ { .fl = -0x008000L }, { .fl = -0.0 }, { .fl = -0x008000L } },
{ { .fl = -0x010000L }, { .fl = -0.0 }, { .fl = -0x010000L } },
{ { .fl = -0x400000L }, { .fl = -0.0 }, { .fl = -0x400000L } },
{ { .fl = -0x7fffffL }, { .fl = -0.0 }, { .fl = -0x7fffffL } },
/* normal too big */
{ { .fl = 0x800000 }, { .fl = 0 }, { .fl = 0x800000 } },
{ { .fl = 0x800001 }, { .fl = 0 }, { .fl = 0x800001 } },
{ { 0x7e000000 }, { .fl = 0 }, { 0x7e000000 } },
{ { 0x7f000000 }, { .fl = 0 }, { 0x7f000000 } },
{ { 0x7f7fffff }, { .fl = 0 }, { 0x7f7fffff } },
{ { .fl = -0x800000L }, { .fl = -0.0 }, { .fl = -0x800000L } },
{ { .fl = -0x800001L }, { .fl = -0.0 }, { .fl = -0x800001L } },
{ { 0xfe000000 }, { .fl = -0.0 }, { 0xfe000000 } },
{ { 0xff000000 }, { .fl = -0.0 }, { 0xff000000 } },
{ { 0xff7fffff }, { .fl = -0.0 }, { 0xff7fffff } },
/* Inf */
{ { 0x7f800000 }, { .fl = 0 }, { 0x7f800000 } },
{ { 0xff800000 }, { .fl = -0.0 }, { 0xff800000 } },
/* integral and fraction */
{ { .fl= 0x1.000002p+0 }, { .fl= 0x0.000002p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000004p+0 }, { .fl= 0x0.000004p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000008p+0 }, { .fl= 0x0.000008p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000010p+0 }, { .fl= 0x0.000010p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000020p+0 }, { .fl= 0x0.000020p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000040p+0 }, { .fl= 0x0.000040p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000080p+0 }, { .fl= 0x0.000080p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000100p+0 }, { .fl= 0x0.000100p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000200p+0 }, { .fl= 0x0.000200p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000400p+0 }, { .fl= 0x0.000400p+0 }, { .fl= 1 } },
{ { .fl= 0x1.000800p+0 }, { .fl= 0x0.000800p+0 }, { .fl= 1 } },
{ { .fl= 0x1.001000p+0 }, { .fl= 0x0.001000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.002000p+0 }, { .fl= 0x0.002000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.004000p+0 }, { .fl= 0x0.004000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.008000p+0 }, { .fl= 0x0.008000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.010000p+0 }, { .fl= 0x0.010000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.020000p+0 }, { .fl= 0x0.020000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.040000p+0 }, { .fl= 0x0.040000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.080000p+0 }, { .fl= 0x0.080000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.100000p+0 }, { .fl= 0x0.100000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.200000p+0 }, { .fl= 0x0.200000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.400000p+0 }, { .fl= 0x0.400000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.800000p+0 }, { .fl= 0x0.800000p+0 }, { .fl= 1 } },
{ { .fl= 0x1.fffffep+0 }, { .fl= 0x0.fffffep+0 }, { .fl= 1 } },
{ { .fl= 0x2.000004p+0 }, { .fl= 0x0.000004p+0 }, { .fl= 2 } },
{ { .fl= 0x2.000400p+0 }, { .fl= 0x0.000400p+0 }, { .fl= 2 } },
{ { .fl= 0x2.800000p+0 }, { .fl= 0x0.800000p+0 }, { .fl= 2 } },
{ { .fl= 0x2.fffffcp+0 }, { .fl= 0x0.fffffcp+0 }, { .fl= 2 } },
{ { .fl= 0x7fffff.8p+0 }, { .fl= 0x0.8p+0 }, { .fl= 0x7fffffp+0 } },
{ { .fl= -0x1.000002p+0 }, { .fl= -0x0.000002p+0 }, { .fl= -1 } },
{ { .fl= -0x1.800000p+0 }, { .fl= -0x0.800000p+0 }, { .fl= -1 } },
{ { .fl= -0x1.fffffep+0 }, { .fl= -0x0.fffffep+0 }, { .fl= -1 } },
{ { .fl= -0x2.000004p+0 }, { .fl= -0x0.000004p+0 }, { .fl= -2 } },
{ { .fl= -0x2.000400p+0 }, { .fl= -0x0.000400p+0 }, { .fl= -2 } },
{ { .fl= -0x2.800000p+0 }, { .fl= -0x0.800000p+0 }, { .fl= -2 } },
{ { .fl= -0x2.fffffcp+0 }, { .fl= -0x0.fffffcp+0 }, { .fl= -2 } },
{ { .fl= -0x7fffff.8p+0 }, { .fl= -0x0.8p+0 }, { .fl= -0x7fffffp+0 } },
};
void x_exit (int index)
{
#ifndef __AVR__
fprintf (stderr, "t[%d]: %#lx\n", index - 1, vf.lo);
#endif
exit (index ? index : -1);
}
int main ()
{
union lofl_u x, zf, zn, vn1;
int i;
for (i = 0; i < (int) (sizeof(t) / sizeof(t[0])); i++) {
x.lo = pgm_read_dword (& t[i].x);
zf.lo = pgm_read_dword (& t[i].f);
zn.lo = pgm_read_dword (& t[i].n);
vf.fl = modf (x.fl, & vn);
vn1.fl = vn;
if ((vf.lo != zf.lo) || (vn1.lo != zn.lo))
x_exit (i+1);
}
return 0;
}
|
