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
156
157
158
159
160
161
162
163
|
/* Test of frexp() function.
$Id: frexp-01.c,v 1.2.2.1 2008/03/20 21:42:29 joerg_wunsch Exp $
*/
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "progmem.h"
#ifndef __AVR__
# define PRINTFLN(fmt, ...) \
printf ("\nLine %d: " fmt "\n", __LINE__, ##__VA_ARGS__)
# define EXIT(code) exit ((code) < 255 ? (code) : 100 + (code) % 100)
#else
# if defined(__AVR_ATmega128__)
/* ATmega128 has enough RAM for sprintf(), print to 0x2000 in XRAM. */
# define PRINTFLN(fmt, ...) \
sprintf ((char *)0x2000, "\nLine %d: " fmt "\n", __LINE__, ##__VA_ARGS__)
# else
/* small AVR */
# define PRINTFLN(args...)
# endif
# define EXIT exit
#endif
union lofl_u {
long lo;
float fl;
};
volatile union lofl_u vf = { 1 };
int ve = { 1 };
PROGMEM const struct { /* Table of test cases. */
union lofl_u x; /* argument */
union lofl_u f; /* fraction part: frexp(x,*pe) */
int e; /* exponent: *pe */
} t[] = {
/* zero */
{ { 0x00000000 }, { 0x00000000 }, 0 },
{ { 0x80000000 }, { 0x80000000 }, 0 },
/* subnormal */
{ { 0x00000001 }, { .fl = 0.5 }, -148 },
{ { 0x00000003 }, { .fl = 0.75 }, -147 },
{ { 0x00000100 }, { .fl = 0.5 }, -140 },
{ { 0x00010000 }, { .fl = 0.5 }, -132 },
{ { 0x007fffff }, { .fl = 0x0.fffffep+00 }, -126 },
{ { 0x80000001 }, { .fl = -0.5 }, -148 },
{ { 0x80000003 }, { .fl = -0.75 }, -147 },
{ { 0x80000100 }, { .fl = -0.5 }, -140 },
{ { 0x80010000 }, { .fl = -0.5 }, -132 },
{ { 0x807fffff }, { .fl = -0x0.fffffep+00 }, -126 },
/* integral normal */
{ { .fl = 0x000001 }, { .fl = 0.5 }, 1 },
{ { .fl = 0x400000 }, { .fl = 0.5 }, 23 },
{ { .fl = 0x7fffff }, { .fl = 0x0.fffffep+00 }, 23 },
{ { .fl = -0x000001 }, { .fl = -0.5 }, 1 },
{ { .fl = -0x400000 }, { .fl = -0.5 }, 23 },
{ { .fl = -0x7fffff }, { .fl = -0x0.fffffep+00 }, 23 },
/* very big */
{ { 0x7f000000 }, { .fl = 0.5 }, 128 },
{ { 0x7f7fffff }, { .fl = 0x0.ffffffp+00 }, 128 },
{ { 0xff000000 }, { .fl = -0.5 }, 128 },
{ { 0xff7fffff }, { .fl = -0x0.ffffffp+00 }, 128 },
};
int main ()
{
union lofl_u x, zf;
int ze;
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);
ze = pgm_read_word (& t[i].e);
vf.fl = frexp (x.fl, & ve);
if (vf.lo != zf.lo || ve != ze) {
PRINTFLN ("t[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 1);
}
}
/* Check all normal values with fractional 0.5 */
for (i = 1; i < 255; i++) {
/* Positive. */
x.lo = (long)i << 23;
vf.fl = frexp (x.fl, & ve);
if (vf.lo != 0x3f000000 || ve != i - 126) {
PRINTFLN ("t[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 100);
}
/* Negative. */
x.lo |= 0x80000000;
vf.fl = frexp (x.fl, & ve);
if ((unsigned long)vf.lo != 0xbf000000 || ve != i - 126) {
PRINTFLN ("t[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 100);
}
}
/* X is not a finite value. */
{
static PROGMEM const long tnf[] = {
0x7f800000, /* +Inf */
0xff800000, /* -Inf */
0x7f800001, /* NaN */
0x7f800100,
0x7f810000,
0x7fc00000,
0x7fffffff,
0xff800001,
0xff800100,
0xff810000,
0xffc00000,
0xffffffff,
};
for (i = 0; i < (int) (sizeof(tnf) / sizeof(tnf[0])); i++) {
x.lo = pgm_read_dword (& tnf[i]);
ve = 1;
vf.fl = frexp (x.fl, & ve);
#ifdef __AVR__
if (vf.lo != x.lo || ve) {
PRINTFLN ("tnf[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 200);
}
#else
if (isinf (x.fl)) {
if (vf.lo != x.lo || ve) {
PRINTFLN ("tnf[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 200);
}
} else {
if (!isnan (vf.fl) || ve) {
PRINTFLN ("tnf[%d]: %#lx %d", i, vf.lo, ve);
EXIT (i + 200);
}
}
#endif
}
}
#ifdef __AVR__
/* NULL is a legal address for Avr-libc realization: skip writing. */
{
unsigned char r1 = 1;
vf.fl = frexp (48, (int *)0);
/* Check __zero_reg__, as it is allocated by 0x0001 address. */
asm ("mov %0,r1 \n\t clr r1" : "=r"(r1));
if (r1) EXIT (300);
if (vf.fl != 0.75) EXIT (301);
}
#endif
return 0;
}
|
