1
2 #ifndef fma_main
3 #define fma_main
4
5 #if DEBUG
6 #include <stdio.h>
7 #endif
8
9 TYPE m1[32] = {
10 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,
11 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32
12 };
13 TYPE m2[32] = {
14 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,
15 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33
16 };
17 TYPE m3[32] = {
18 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
19 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34
20 };
21 TYPE m4[32];
22 int test_fails = 0;
23
24 void
compare_result(char * title,TYPE * res)25 compare_result(char * title, TYPE *res)
26 {
27 int i;
28 int good = 1;
29 for (i =0; i < 32; i++)
30 if (m4[i] != res[i])
31 {
32 if (good)
33 {
34 #if DEBUG
35 printf ("!!!! %s miscompare\n", title);
36 #endif
37 good = 0;
38 }
39 #if DEBUG
40 printf ("res[%d] = %d, must be %d\n", i, (int)res[i], (int) m4[i]);
41 #endif
42 }
43 if (!good)
44 test_fails = 1;
45 }
46
fma_test()47 static void fma_test ()
48 {
49 int i;
50 for (i=0; i <32; i++)
51 m4[i] = test_noneg_add_noneg_add (m1[i], m2[i], m3[i]);
52 compare_result ("test0000", res_test0000);
53
54 for (i=0; i <32; i++)
55 m4[i] = test_noneg_add_noneg_sub (m1[i], m2[i], m3[i]);
56 compare_result ("test0001", res_test0001);
57
58 for (i=0; i <32; i++)
59 m4[i] = test_noneg_add_neg_add (m1[i], m2[i], m3[i]);
60 compare_result ("test0010", res_test0010);
61
62 for (i=0; i <32; i++)
63 m4[i] = test_noneg_add_neg_sub (m1[i], m2[i], m3[i]);
64 compare_result ("test0011", res_test0011);
65
66 for (i=0; i <32; i++)
67 m4[i] = test_noneg_sub_noneg_add (m1[i], m2[i], m3[i]);
68 compare_result ("test0100", res_test0100);
69
70 for (i=0; i <32; i++)
71 m4[i] = test_noneg_sub_noneg_sub (m1[i], m2[i], m3[i]);
72 compare_result ("test0101", res_test0101);
73
74 for (i=0; i <32; i++)
75 m4[i] = test_noneg_sub_neg_add (m1[i], m2[i], m3[i]);
76 compare_result ("test0110", res_test0110);
77
78 for (i=0; i <32; i++)
79 m4[i] = test_noneg_sub_neg_sub (m1[i], m2[i], m3[i]);
80 compare_result ("test0111", res_test0111);
81
82 for (i=0; i <32; i++)
83 m4[i] = test_neg_add_noneg_add (m1[i], m2[i], m3[i]);
84 compare_result ("test1000", res_test1000);
85
86 for (i=0; i <32; i++)
87 m4[i] = test_neg_add_noneg_sub (m1[i], m2[i], m3[i]);
88 compare_result ("test1001", res_test1001);
89
90 for (i=0; i <32; i++)
91 m4[i] = test_neg_add_neg_add (m1[i], m2[i], m3[i]);
92 compare_result ("test1010", res_test1010);
93
94 for (i=0; i <32; i++)
95 m4[i] = test_neg_add_neg_sub (m1[i], m2[i], m3[i]);
96 compare_result ("test1011", res_test1011);
97
98 for (i=0; i <32; i++)
99 m4[i] = test_neg_sub_noneg_add (m1[i], m2[i], m3[i]);
100 compare_result ("test1100", res_test1100);
101
102 for (i=0; i <32; i++)
103 m4[i] = test_neg_sub_noneg_sub (m1[i], m2[i], m3[i]);
104 compare_result ("test1101", res_test1101);
105
106 for (i=0; i <32; i++)
107 m4[i] = test_neg_sub_neg_add (m1[i], m2[i], m3[i]);
108 compare_result ("test1110", res_test1110);
109
110 for (i=0; i <32; i++)
111 m4[i] = test_neg_sub_neg_sub (m1[i], m2[i], m3[i]);
112 compare_result ("test1111", res_test1111);
113
114 if (test_fails) abort ();
115 }
116
117 #endif
118