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
|
/* Test of hypot() function.
$Id: hypot-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 v = { .lo = 1 };
PROGMEM const struct { /* Table of test cases */
union lofl_u x, y; /* args */
union lofl_u z; /* result */
} t[] = {
/* NaN and NaN --> NaN */
{ { 0x7f800001 }, { 0x7f800100 }, { 0x7fc00000 } },
{ { 0x7f810001 }, { 0x7fc00000 }, { 0x7fc00000 } },
{ { 0x7fffffff }, { 0xff800001 }, { 0x7fc00000 } },
{ { 0xff800100 }, { 0x7f800100 }, { 0x7fc00000 } },
/* NaN and finite --> NaN */
{ { 0x7f800001 }, { 0x3f800000 }, { 0x7fc00000 } },
{ { 0x3f800000 }, { 0x7f800001 }, { 0x7fc00000 } },
/* Inf & NaN --> Inf */
{ { 0x7f800000 }, { 0x7f800001 }, { 0x7f800000 } },
{ { 0x7f800000 }, { 0x7f800100 }, { 0x7f800000 } },
{ { 0x7f800000 }, { 0x7f810000 }, { 0x7f800000 } },
{ { 0x7f800000 }, { 0x7fc00000 }, { 0x7f800000 } },
{ { 0x7f800000 }, { 0x7fffffff }, { 0x7f800000 } },
{ { 0x7f800000 }, { 0xffffffff }, { 0x7f800000 } },
{ { 0xff800000 }, { 0x7fc00000 }, { 0x7f800000 } },
{ { 0x7fc00000 }, { 0x7f800000 }, { 0x7f800000 } },
{ { 0x7fc00000 }, { 0xff800000 }, { 0x7f800000 } },
/* 0 and 0 */
{ { 0x00000000 }, { 0x00000000 }, { 0x00000000 } },
{ { 0x80000000 }, { 0x00000000 }, { 0x00000000 } },
{ { 0x00000000 }, { 0x80000000 }, { 0x00000000 } },
{ { 0x80000000 }, { 0x80000000 }, { 0x00000000 } },
/* hypot(0,*) */
{ { 0 }, { 0x00000001 }, { 0x00000001 } },
{ { 0 }, { 0x007fffff }, { 0x007fffff } },
{ { 0 }, { 0x00800000 }, { 0x00800000 } },
{ { 0 }, { 0x3f800000 }, { 0x3f800000 } },
{ { 0 }, { 0x7f7fffff }, { 0x7f7fffff } },
{ { 0 }, { 0x7f800000 }, { 0x7f800000 } }, /* Inf */
/* hypot(*,0) */
{ { 0x00000001 }, { 0 }, { 0x00000001 } },
{ { 0x007fffff }, { 0 }, { 0x007fffff } },
{ { 0x00800000 }, { 0 }, { 0x00800000 } },
{ { 0x3f800000 }, { 0 }, { 0x3f800000 } },
{ { 0x7f7fffff }, { 0 }, { 0x7f7fffff } },
{ { 0x7f800000 }, { 0 }, { 0x7f800000 } }, /* Inf */
/* overflow */
{ { 0x7f7fffff }, { 0x7f7fffff }, { 0x7f800000 } },
{ { 0xff7fffff }, { 0xff7fffff }, { 0x7f800000 } },
/* { { .fl= 2.4061596e+38 }, { .fl= 2.4061596e+38 }, { 0x7f800000 } }, */
{ { .fl= 2.4061598e+38 }, { .fl= 2.4061598e+38 }, { 0x7f800000 } },
{ { .fl= 3 }, { .fl= 4 }, { .fl= 5 } },
{ { .fl= -3 }, { .fl= 4 }, { .fl= 5 } },
{ { .fl= 3 }, { .fl= -4 }, { .fl= 5 } },
{ { .fl= -3 }, { .fl= -4 }, { .fl= 5 } },
/* very small args */
{ { 0x00000001 }, { 0x00000001 }, { 0x00000001 } }, /* 1.4 */
{ { 0x00000001 }, { 0x00000002 }, { 0x00000002 } }, /* 2.2 */
{ { 0x00000002 }, { 0x00000001 }, { 0x00000002 } }, /* 2.2 */
{ { 0x00000002 }, { 0x00000002 }, { 0x00000003 } }, /* 2.8 */
/* very big args */
{ { .fl= 2.4061596e+38 }, { .fl= 2.4061595e+38 }, { 0x7f7fffff } },
{ { .fl= 2.4061595e+38 }, { .fl= 2.4061596e+38 }, { 0x7f7fffff } },
/* x >> y */
{ { .fl= 1 }, { .fl= 0x0.8p-9 }, { 0x3f800004 } },
{ { .fl= 1 }, { .fl= 0x0.8p-10 }, { 0x3f800001 } },
{ { .fl= 1 }, { .fl= 0x0.8p-11 }, { 0x3f800000 } },
{ { .fl= 1 }, { .fl= 0x0.8p-12 }, { 0x3f800000 } },
{ { .fl= 1 }, { .fl= 0x0.8p-13 }, { 0x3f800000 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-10 }, { 0x3f000011 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-11 }, { 0x3f000005 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-12 }, { 0x3f000002 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-13 }, { 0x3f000001 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-14 }, { 0x3f000001 } },
{ { .fl= 0x0.800001p0 }, { .fl= 0x0.ffffffp-15 }, { 0x3f000001 } },
/* x << y */
{ { .fl= 0x0.8p-9 }, { .fl= 1 }, { 0x3f800004 } },
{ { .fl= 0x0.8p-10 }, { .fl= 1 }, { 0x3f800001 } },
{ { .fl= 0x0.8p-11 }, { .fl= 1 }, { 0x3f800000 } },
{ { .fl= 0x0.8p-12 }, { .fl= 1 }, { 0x3f800000 } },
{ { .fl= 0x0.8p-13 }, { .fl= 1 }, { 0x3f800000 } },
{ { .fl= 0x0.ffffffp-10 }, { .fl= 0x0.800001p0 }, { 0x3f000011 } },
{ { .fl= 0x0.ffffffp-11 }, { .fl= 0x0.800001p0 }, { 0x3f000005 } },
{ { .fl= 0x0.ffffffp-12 }, { .fl= 0x0.800001p0 }, { 0x3f000002 } },
{ { .fl= 0x0.ffffffp-13 }, { .fl= 0x0.800001p0 }, { 0x3f000001 } },
{ { .fl= 0x0.ffffffp-14 }, { .fl= 0x0.800001p0 }, { 0x3f000001 } },
{ { .fl= 0x0.ffffffp-15 }, { .fl= 0x0.800001p0 }, { 0x3f000001 } },
/* a set of scales */
{ { 0x00000003 }, { 0x00000004 }, { 0x00000005 } },
{ { .fl= 0x0.3p-140 }, { .fl= 0x0.4p-140 }, { .fl= 0x0.5p-140 } },
{ { .fl= 0x0.3p-120 }, { .fl= 0x0.4p-120 }, { .fl= 0x0.5p-120 } },
{ { .fl= 0x0.3p-100 }, { .fl= 0x0.4p-100 }, { .fl= 0x0.5p-100 } },
{ { .fl= 0x0.3p-80 }, { .fl= 0x0.4p-80 }, { .fl= 0x0.5p-80 } },
{ { .fl= 0x0.3p-60 }, { .fl= 0x0.4p-60 }, { .fl= 0x0.5p-60 } },
{ { .fl= 0x0.3p-40 }, { .fl= 0x0.4p-40 }, { .fl= 0x0.5p-40 } },
{ { .fl= 0x0.3p-20 }, { .fl= 0x0.4p-20 }, { .fl= 0x0.5p-20 } },
{ { .fl= 0x0.3p+0 }, { .fl= 0x0.4p+0 }, { .fl= 0x0.5p+0 } },
{ { .fl= 0x0.3p+20 }, { .fl= 0x0.4p+20 }, { .fl= 0x0.5p+20 } },
{ { .fl= 0x0.3p+40 }, { .fl= 0x0.4p+40 }, { .fl= 0x0.5p+40 } },
{ { .fl= 0x0.3p+60 }, { .fl= 0x0.4p+60 }, { .fl= 0x0.5p+60 } },
{ { .fl= 0x0.3p+80 }, { .fl= 0x0.4p+80 }, { .fl= 0x0.5p+80 } },
{ { .fl= 0x0.3p+100 }, { .fl= 0x0.4p+100 }, { .fl= 0x0.5p+100 } },
{ { .fl= 0x0.3p+120 }, { .fl= 0x0.4p+120 }, { .fl= 0x0.5p+120 } },
};
void x_exit (int index)
{
#ifndef __AVR__
fprintf (stderr, "t[%d]: %#lx\n", index - 1, v.lo);
#endif
exit (index ? index : -1);
}
int main ()
{
union lofl_u x,y,z;
int i;
for (i = 0; i < (int) (sizeof(t) / sizeof(t[0])); i++) {
y.lo = pgm_read_dword (& t[i].y);
x.lo = pgm_read_dword (& t[i].x);
z.lo = pgm_read_dword (& t[i].z);
v.fl = hypot (x.fl, y.fl);
if ( v.lo == z.lo
|| (isnan(v.fl) && isnan(z.fl)) )
continue;
x_exit (i+1);
}
return 0;
}
|
