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
|
/* Float tests.
operation: ADD, SUB, MUL, DIV, REVDIV
*/
#if 1
// we are in the regression tests
#include <testfwk.h>
#define DEBUG(x)
#define {operation} 1
#else
// we are standalone
#include <stdio.h>
#define DEBUG(x) x
#define ASSERT(x)
#define ADD 1
#define SUB 1
#define MUL 1
#define DIV 1
#define REVDIV 1
#endif
#ifdef SDCC_mcs51
# define STORAGE xdata
# define XDATA xdata
#elif SDCC_pic16
# define STORAGE code
# define XDATA
#else
# define STORAGE
# define XDATA
#endif
XDATA volatile float left, right, result;
struct {
float left, right, add, sub, mul, div, revdiv;
} STORAGE cases[]={
// left right add sub mul div revdiv
{ 12.8, 25.6, 38.4, -12.8, 327.68, 0.5, 2},
{ 12.8, -25.6, -12.8, 38.4, -327.68, -0.5, -2},
{ -12.8, 25.6, 12.8, -38.4, -327.68, -0.5, -2},
{ -12.8, -25.6, -38.4, 12.8, 327.68, 0.5, 2},
{ 100.0, 10.0, 110.0, 90.0, 1000.00, 10.0, 0.1},
{ 1000.0, 10.0, 1010.0, 990.0, 10000.00, 100.0, 0.01},
{ 10000.0, 10.0, 10010.0, 9990.0, 100000.00, 1000.0, 0.001},
{ 100000.0, 10.0, 100010.0, 99990.0, 1000000.00, 10000.0, 0.0001},
{ 1000000.0, 10.0, 1000010.0, 999990.0, 10000000.00, 100000.0, 0.00001},
{10000000.0, 10.0,10000010.0, 9999990.0,100000000.00, 1000000.0, 0.000001},
{ 0x100, 0x10, 0x110, 0xf0, 0x1000, 0x10, 0.0625},
{ 0x1000, 0x10, 0x1010, 0xff0, 0x10000, 0x100, 0.00390625},
{ 0x10000, 0x10, 0x10010, 0xfff0, 0x100000, 0x1000, 0.00024414},
{ 0x100000, 0x10, 0x100010, 0xffff0, 0x1000000, 0x10000, 0 /* ignore */},
{ 0x1000000, 0x10, 0x1000010, 0xfffff0, 0x10000000, 0x100000, 0 /* ignore */},
{0x10000000, 0x10,0x10000010, 0xffffff0, (float)0x10000000*0x10,
0x1000000, 0 /* ignore */},
};
XDATA int tests=0, errors=0;
char compare (float is, float should) {
float diff = should ? is/should : 0;
tests++;
DEBUG(printf (" %1.3f (%f %f) ", is, should, diff));
if (should==0) {
DEBUG(printf ("IGNORED!\n"));
return 0;
}
// skip the fp roundoff errors
if (diff>0.999999 && diff<1.00001) {
DEBUG(printf ("OK!\n"));
ASSERT(1);
return 0;
} else {
errors++;
DEBUG(printf ("FAIL!\n"));
ASSERT(0);
return 1;
}
}
void
testFloatMath(void)
{
int i;
int t = sizeof(cases)/sizeof(cases[0]);
float result;
for (i=0; i<t; i++) {
#ifdef ADD
// add
result=cases[i].left+cases[i].right;
DEBUG(printf ("%1.3f + %1.3f =", cases[i].left, cases[i].right));
compare(result, cases[i].add);
#endif
#ifdef SUB
// sub
result=cases[i].left-cases[i].right;
DEBUG(printf ("%1.3f - %1.3f =", cases[i].left, cases[i].right));
compare(result, cases[i].sub);
#endif
#ifdef MUL
// mul
result=cases[i].left*cases[i].right;
DEBUG(printf ("%1.3f * %1.3f =", cases[i].left, cases[i].right));
compare(result, cases[i].mul);
#endif
#ifdef DIV
// div
result=cases[i].left/cases[i].right;
DEBUG(printf ("%1.3f / %1.3f =", cases[i].left, cases[i].right));
compare(result, cases[i].div);
#endif
#ifdef REVDIV
// revdiv
result=cases[i].right/cases[i].left;
DEBUG(printf ("%1.3f / %1.3f =", cases[i].right, cases[i].left));
compare(result, cases[i].revdiv);
#endif
}
DEBUG(printf ("%d tests, %d errors\n", tests, errors));
}
#ifndef REENTRANT
void main(void) {
testFloatMath();
}
#endif
|
