1 #define vector(elcount, type) \
2 __attribute__((vector_size((elcount)*sizeof(type)))) type
3
4 #define vidx(type, vec, idx) (*((type *) &(vec) + idx))
5
6 #define operl(a, b, op) (a op b)
7 #define operr(a, b, op) (b op a)
8
9 #define check(type, count, vec0, vec1, num, op, lr) \
10 do {\
11 int __i; \
12 for (__i = 0; __i < count; __i++) {\
13 if (vidx (type, vec1, __i) != oper##lr (num, vidx (type, vec0, __i), op)) \
14 __builtin_abort (); \
15 }\
16 } while (0)
17
18 #define veccompare(type, count, v0, v1) \
19 do {\
20 int __i; \
21 for (__i = 0; __i < count; __i++) { \
22 if (vidx (type, v0, __i) != vidx (type, v1, __i)) \
23 __builtin_abort (); \
24 } \
25 } while (0)
26
27 volatile int one = 1;
28
main(int argc,char * argv[])29 int main (int argc, char *argv[]) {
30 #define fvec_2 (vector(4, float)){2., 2., 2., 2.}
31 #define dvec_2 (vector(2, double)){2., 2.}
32
33
34 vector(8, short) v0 = {one, 1, 2, 3, 4, 5, 6, 7};
35 vector(8, short) v1;
36
37 vector(4, float) f0 = {1., 2., 3., 4.};
38 vector(4, float) f1, f2;
39
40 vector(2, double) d0 = {1., 2.};
41 vector(2, double) d1, d2;
42
43
44
45 v1 = 2 + v0; check (short, 8, v0, v1, 2, +, l);
46 v1 = 2 - v0; check (short, 8, v0, v1, 2, -, l);
47 v1 = 2 * v0; check (short, 8, v0, v1, 2, *, l);
48 v1 = 2 / v0; check (short, 8, v0, v1, 2, /, l);
49 v1 = 2 % v0; check (short, 8, v0, v1, 2, %, l);
50 v1 = 2 ^ v0; check (short, 8, v0, v1, 2, ^, l);
51 v1 = 2 & v0; check (short, 8, v0, v1, 2, &, l);
52 v1 = 2 | v0; check (short, 8, v0, v1, 2, |, l);
53 v1 = 2 << v0; check (short, 8, v0, v1, 2, <<, l);
54 v1 = 2 >> v0; check (short, 8, v0, v1, 2, >>, l);
55
56 v1 = v0 + 2; check (short, 8, v0, v1, 2, +, r);
57 v1 = v0 - 2; check (short, 8, v0, v1, 2, -, r);
58 v1 = v0 * 2; check (short, 8, v0, v1, 2, *, r);
59 v1 = v0 / 2; check (short, 8, v0, v1, 2, /, r);
60 v1 = v0 % 2; check (short, 8, v0, v1, 2, %, r);
61 v1 = v0 ^ 2; check (short, 8, v0, v1, 2, ^, r);
62 v1 = v0 & 2; check (short, 8, v0, v1, 2, &, r);
63 v1 = v0 | 2; check (short, 8, v0, v1, 2, |, r);
64
65 f1 = 2. + f0; f2 = fvec_2 + f0; veccompare (float, 4, f1, f2);
66 f1 = 2. - f0; f2 = fvec_2 - f0; veccompare (float, 4, f1, f2);
67 f1 = 2. * f0; f2 = fvec_2 * f0; veccompare (float, 4, f1, f2);
68 f1 = 2. / f0; f2 = fvec_2 / f0; veccompare (float, 4, f1, f2);
69
70 f1 = f0 + 2.; f2 = f0 + fvec_2; veccompare (float, 4, f1, f2);
71 f1 = f0 - 2.; f2 = f0 - fvec_2; veccompare (float, 4, f1, f2);
72 f1 = f0 * 2.; f2 = f0 * fvec_2; veccompare (float, 4, f1, f2);
73 f1 = f0 / 2.; f2 = f0 / fvec_2; veccompare (float, 4, f1, f2);
74
75 d1 = 2. + d0; d2 = dvec_2 + d0; veccompare (double, 2, d1, d2);
76 d1 = 2. - d0; d2 = dvec_2 - d0; veccompare (double, 2, d1, d2);
77 d1 = 2. * d0; d2 = dvec_2 * d0; veccompare (double, 2, d1, d2);
78 d1 = 2. / d0; d2 = dvec_2 / d0; veccompare (double, 2, d1, d2);
79
80 d1 = d0 + 2.; d2 = d0 + dvec_2; veccompare (double, 2, d1, d2);
81 d1 = d0 - 2.; d2 = d0 - dvec_2; veccompare (double, 2, d1, d2);
82 d1 = d0 * 2.; d2 = d0 * dvec_2; veccompare (double, 2, d1, d2);
83 d1 = d0 / 2.; d2 = d0 / dvec_2; veccompare (double, 2, d1, d2);
84
85 return 0;
86 }
87