Mercurial > vec

comparison src/impl/fallback.c @ 23:e26874655738

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
*: huge refactor, new major release (hahaha) I keep finding things that are broken... The problem NOW was that vec would unintentionally build some functions with extended instruction sets, which is Bad and would mean that for all intents and purposes the CPU detection was completely broken. Now vec is no longer header only either. Boohoo. However this gives a lot more flexibility to vec since we no longer want or need to care about C++ crap. The NEON and Altivec implementations have not been updated which means they won't compile hence why they're commented out in the cmake build file.
author Paper <paper@tflc.us>
date Sun, 24 Nov 2024 02:52:40 -0500
parents
children e49e70f7012f
comparison
equal deleted inserted replaced
22:fbcd3fa6f8fc 23:e26874655738
1 #include "vec/impl/fallback.h"
2
3 #include <string.h>
4
5 // Fallback implementations - this is what an implementation should use if it
6 // doesn't support a specific function *and* the actual representation in
7 // memory is unknown or yields incorrect results from the generic functions.
8 // This is *extremely* unlikely; for x86 the layout is exactly the same in
9 // memory as the generic functions (i.e. it is literally stored as an array of
10 // integers).
11 //
12 // These functions can probably be removed if every current implementation is
13 // found to have the same
14
15 #define VEC_FALLBACK_OPERATION(op, sign, csign, bits, size) \
16 do { \
17 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(varr1); \
18 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(varr2); \
19 \
20 v##sign##int##bits##x##size##_store_aligned(vec1, varr1); \
21 v##sign##int##bits##x##size##_store_aligned(vec2, varr2); \
22 \
23 for (int i = 0; i < size; i++) varr1[i] = (op); \
24 \
25 return v##sign##int##bits##x##size##_load_aligned(varr1); \
26 } while (0)
27
28 #define VEC_FALLBACK_CMP(op, sign, csign, bits, size) \
29 VEC_FALLBACK_OPERATION((varr1[i] op varr2[i]) ? VEC_UINT##bits##_MAX : 0, sign, csign, bits, size)
30
31 #define VEC_FALLBACK_SHIFT(op, sign, csign, bits, size) \
32 do { \
33 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(varr1); \
34 VUINT##bits##x##size##_ALIGNED_ARRAY(varr2); \
35 \
36 v##sign##int##bits##x##size##_store_aligned(vec1, varr1); \
37 vuint##bits##x##size##_store_aligned(vec2, varr2); \
38 \
39 for (int i = 0; i < size; i++) varr1[i] = (op); \
40 \
41 return v##sign##int##bits##x##size##_load_aligned(varr1); \
42 } while (0)
43
44 #define VEC_DEFINE_FALLBACK_OPERATIONS_SIGN(sign, csign, bits, size) \
45 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_splat(vec_##sign##int##bits x) \
46 { \
47 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(arr); \
48 for (int i = 0; i < size; i++) arr[i] = x; \
49 return v##sign##int##bits##x##size##_load_aligned(arr); \
50 } \
51 \
52 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_load(const vec_##sign##int##bits in[size]) \
53 { \
54 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(arr); \
55 memcpy(arr, in, sizeof(vec_##sign##int##bits) * size); \
56 return v##sign##int##bits##x##size##_load_aligned(arr); \
57 } \
58 \
59 void v##sign##int##bits##x##size##_fallback_store(v##sign##int##bits##x##size vec, vec_##sign##int##bits out[size]) \
60 { \
61 V##csign##INT##bits##x##size##_ALIGNED_ARRAY(arr); \
62 v##sign##int##bits##x##size##_store_aligned(vec, arr); \
63 memcpy(out, arr, sizeof(vec_##sign##int##bits) * size); \
64 } \
65 \
66 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_add(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
67 { \
68 VEC_FALLBACK_OPERATION(varr1[i] + varr2[i], sign, csign, bits, size); \
69 } \
70 \
71 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_sub(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
72 { \
73 VEC_FALLBACK_OPERATION(varr1[i] - varr2[i], sign, csign, bits, size); \
74 } \
75 \
76 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_mul(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
77 { \
78 VEC_FALLBACK_OPERATION(varr1[i] * varr2[i], sign, csign, bits, size); \
79 } \
80 \
81 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_div(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
82 { \
83 VEC_FALLBACK_OPERATION(varr2[i] ? (varr1[i] / varr2[i]) : 0, sign, csign, bits, size); \
84 } \
85 \
86 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_avg(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
87 { \
88 VEC_FALLBACK_OPERATION((varr1[i] + varr2[i] + 1) / 2, sign, csign, bits, size); \
89 } \
90 \
91 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_and(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
92 { \
93 VEC_FALLBACK_OPERATION(varr1[i] & varr2[i], sign, csign, bits, size); \
94 } \
95 \
96 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_or(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
97 { \
98 VEC_FALLBACK_OPERATION(varr1[i] | varr2[i], sign, csign, bits, size); \
99 } \
100 \
101 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_xor(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
102 { \
103 VEC_FALLBACK_OPERATION(varr1[i] ^ varr2[i], sign, csign, bits, size); \
104 } \
105 \
106 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_not(v##sign##int##bits##x##size vec) \
107 { \
108 return v##sign##int##bits##x##size##_xor(vec, v##sign##int##bits##x##size##_splat((vec_##sign##int##bits)VEC_UINT##bits##_MAX)); \
109 } \
110 \
111 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_cmplt(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
112 { \
113 VEC_FALLBACK_CMP(<, sign, csign, bits, size); \
114 } \
115 \
116 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_cmple(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
117 { \
118 VEC_FALLBACK_CMP(<=, sign, csign, bits, size); \
119 } \
120 \
121 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_cmpeq(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
122 { \
123 VEC_FALLBACK_CMP(==, sign, csign, bits, size); \
124 } \
125 \
126 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_cmpge(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
127 { \
128 VEC_FALLBACK_CMP(>=, sign, csign, bits, size); \
129 } \
130 \
131 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_cmpgt(v##sign##int##bits##x##size vec1, v##sign##int##bits##x##size vec2) \
132 { \
133 VEC_FALLBACK_CMP(>, sign, csign, bits, size); \
134 } \
135 \
136 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_lshift(v##sign##int##bits##x##size vec1, vuint##bits##x##size vec2) \
137 { \
138 VEC_FALLBACK_SHIFT(vec_##sign##lshift(varr1[i], varr2[i]), sign, csign, bits, size); \
139 } \
140 \
141 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_rshift(v##sign##int##bits##x##size vec1, vuint##bits##x##size vec2) \
142 { \
143 VEC_FALLBACK_SHIFT(vec_##sign##rshift(varr1[i], varr2[i]), sign, csign, bits, size); \
144 } \
145 \
146 v##sign##int##bits##x##size v##sign##int##bits##x##size##_fallback_lrshift(v##sign##int##bits##x##size vec1, vuint##bits##x##size vec2) \
147 { \
148 VEC_FALLBACK_SHIFT(vec_lrshift((vec_uint##bits)varr1[i], varr2[i]), sign, csign, bits, size); \
149 }
150
151 #define VEC_DEFINE_FALLBACK_OPERATIONS(bits, size) \
152 VEC_DEFINE_FALLBACK_OPERATIONS_SIGN( , , bits, size) \
153 VEC_DEFINE_FALLBACK_OPERATIONS_SIGN(u, U, bits, size)
154
155 // 16-bit
156 VEC_DEFINE_FALLBACK_OPERATIONS(8, 2)
157
158 // 32-bit
159 VEC_DEFINE_FALLBACK_OPERATIONS(8, 4)
160 VEC_DEFINE_FALLBACK_OPERATIONS(16, 2)
161
162 // 64-bit
163 VEC_DEFINE_FALLBACK_OPERATIONS(8, 8)
164 VEC_DEFINE_FALLBACK_OPERATIONS(16, 4)
165 VEC_DEFINE_FALLBACK_OPERATIONS(32, 2)
166
167 // 128-bit
168 VEC_DEFINE_FALLBACK_OPERATIONS(8, 16)
169 VEC_DEFINE_FALLBACK_OPERATIONS(16, 8)
170 VEC_DEFINE_FALLBACK_OPERATIONS(32, 4)
171 VEC_DEFINE_FALLBACK_OPERATIONS(64, 2)
172
173 // 256-bit
174 VEC_DEFINE_FALLBACK_OPERATIONS(8, 32)
175 VEC_DEFINE_FALLBACK_OPERATIONS(16, 16)
176 VEC_DEFINE_FALLBACK_OPERATIONS(32, 8)
177 VEC_DEFINE_FALLBACK_OPERATIONS(64, 4)
178
179 // 512-bit
180 VEC_DEFINE_FALLBACK_OPERATIONS(8, 64)
181 VEC_DEFINE_FALLBACK_OPERATIONS(16, 32)
182 VEC_DEFINE_FALLBACK_OPERATIONS(32, 16)
183 VEC_DEFINE_FALLBACK_OPERATIONS(64, 8)
184
185 #undef VEC_FALLBACK_OPERATION
186 #undef VEC_FALLBACK_CMP
187 #undef VEC_FALLBACK_SHIFT
188 #undef VEC_DEFINE_FALLBACK_OPERATIONS
189 #undef VEC_DEFINE_FALLBACK_OPERATIONS_SIGN