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
|
#ifndef fma_main
#define fma_main
#if DEBUG
#include <stdio.h>
#endif
TYPE m1[32] = {
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
};
TYPE m2[32] = {
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
};
TYPE m3[32] = {
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
};
TYPE m4[32];
int test_fails = 0;
void
compare_result(char * title, TYPE *res)
{
int i;
int good = 1;
for (i =0; i < 32; i++)
if (m4[i] != res[i])
{
if (good)
{
#if DEBUG
printf ("!!!! %s miscompare\n", title);
#endif
good = 0;
}
#if DEBUG
printf ("res[%d] = %d, must be %d\n", i, (int)res[i], (int) m4[i]);
#endif
}
if (!good)
test_fails = 1;
}
static void fma_test ()
{
int i;
for (i=0; i <32; i++)
m4[i] = test_noneg_add_noneg_add (m1[i], m2[i], m3[i]);
compare_result ("test0000", res_test0000);
for (i=0; i <32; i++)
m4[i] = test_noneg_add_noneg_sub (m1[i], m2[i], m3[i]);
compare_result ("test0001", res_test0001);
for (i=0; i <32; i++)
m4[i] = test_noneg_add_neg_add (m1[i], m2[i], m3[i]);
compare_result ("test0010", res_test0010);
for (i=0; i <32; i++)
m4[i] = test_noneg_add_neg_sub (m1[i], m2[i], m3[i]);
compare_result ("test0011", res_test0011);
for (i=0; i <32; i++)
m4[i] = test_noneg_sub_noneg_add (m1[i], m2[i], m3[i]);
compare_result ("test0100", res_test0100);
for (i=0; i <32; i++)
m4[i] = test_noneg_sub_noneg_sub (m1[i], m2[i], m3[i]);
compare_result ("test0101", res_test0101);
for (i=0; i <32; i++)
m4[i] = test_noneg_sub_neg_add (m1[i], m2[i], m3[i]);
compare_result ("test0110", res_test0110);
for (i=0; i <32; i++)
m4[i] = test_noneg_sub_neg_sub (m1[i], m2[i], m3[i]);
compare_result ("test0111", res_test0111);
for (i=0; i <32; i++)
m4[i] = test_neg_add_noneg_add (m1[i], m2[i], m3[i]);
compare_result ("test1000", res_test1000);
for (i=0; i <32; i++)
m4[i] = test_neg_add_noneg_sub (m1[i], m2[i], m3[i]);
compare_result ("test1001", res_test1001);
for (i=0; i <32; i++)
m4[i] = test_neg_add_neg_add (m1[i], m2[i], m3[i]);
compare_result ("test1010", res_test1010);
for (i=0; i <32; i++)
m4[i] = test_neg_add_neg_sub (m1[i], m2[i], m3[i]);
compare_result ("test1011", res_test1011);
for (i=0; i <32; i++)
m4[i] = test_neg_sub_noneg_add (m1[i], m2[i], m3[i]);
compare_result ("test1100", res_test1100);
for (i=0; i <32; i++)
m4[i] = test_neg_sub_noneg_sub (m1[i], m2[i], m3[i]);
compare_result ("test1101", res_test1101);
for (i=0; i <32; i++)
m4[i] = test_neg_sub_neg_add (m1[i], m2[i], m3[i]);
compare_result ("test1110", res_test1110);
for (i=0; i <32; i++)
m4[i] = test_neg_sub_neg_sub (m1[i], m2[i], m3[i]);
compare_result ("test1111", res_test1111);
if (test_fails) abort ();
}
#endif
|
