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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
|
/*
* Copyright 2020 The Emscripten Authors. All rights reserved.
* Emscripten is available under two separate licenses, the MIT license and the
* University of Illinois/NCSA Open Source License. Both these licenses can be
* found in the LICENSE file.
*/
// This file uses SSE1 by calling different functions with different interesting inputs and prints the results.
// Use a diff tool to compare the results between platforms.
#include <xmmintrin.h>
#include "test_sse.h"
bool testNaNBits = true;
int main()
{
float *interesting_floats = get_interesting_floats();
int numInterestingFloats = sizeof(interesting_floats_)/sizeof(interesting_floats_[0]);
assert(numInterestingFloats % 4 == 0);
uint32_t *interesting_ints = get_interesting_ints();
int numInterestingInts = sizeof(interesting_ints_)/sizeof(interesting_ints_[0]);
assert(numInterestingInts % 4 == 0);
// MXCSR register state
printf("%08X\n", _mm_getcsr());
printf("%08X\n", _MM_GET_EXCEPTION_MASK());
printf("%08X\n", _MM_GET_EXCEPTION_STATE());
printf("%08X\n", _MM_GET_FLUSH_ZERO_MODE());
printf("%08X\n", _MM_GET_ROUNDING_MODE());
// SSE1 Arithmetic instructions:
Ret_M128_M128(__m128, _mm_add_ps);
Ret_M128_M128(__m128, _mm_add_ss);
Ret_M128_M128(__m128, _mm_div_ps);
Ret_M128_M128(__m128, _mm_div_ss);
Ret_M128_M128(__m128, _mm_mul_ps);
Ret_M128_M128(__m128, _mm_mul_ss);
Ret_M128_M128(__m128, _mm_sub_ps);
Ret_M128_M128(__m128, _mm_sub_ss);
// SSE1 Elementary Math functions:
Ret_M128approx(__m128, _mm_rcp_ps);
Ret_M128approx(__m128, _mm_rcp_ss);
Ret_M128approx(__m128, _mm_rsqrt_ps);
Ret_M128approx(__m128, _mm_rsqrt_ss);
Ret_M128approx(__m128, _mm_sqrt_ps);
Ret_M128approx(__m128, _mm_sqrt_ss);
// SSE1 Logical instructions:
Ret_M128_M128(__m128, _mm_and_ps);
Ret_M128_M128(__m128, _mm_andnot_ps);
Ret_M128_M128(__m128, _mm_or_ps);
Ret_M128_M128(__m128, _mm_xor_ps);
// SSE1 Compare instructions:
Ret_M128_M128(__m128, _mm_cmpeq_ps);
Ret_M128_M128(__m128, _mm_cmpeq_ss);
Ret_M128_M128(__m128, _mm_cmpge_ps);
Ret_M128_M128(__m128, _mm_cmpge_ss);
Ret_M128_M128(__m128, _mm_cmpgt_ps);
Ret_M128_M128(__m128, _mm_cmpgt_ss);
Ret_M128_M128(__m128, _mm_cmple_ps);
Ret_M128_M128(__m128, _mm_cmple_ss);
Ret_M128_M128(__m128, _mm_cmplt_ps);
Ret_M128_M128(__m128, _mm_cmplt_ss);
Ret_M128_M128(__m128, _mm_cmpneq_ps);
Ret_M128_M128(__m128, _mm_cmpneq_ss);
Ret_M128_M128(__m128, _mm_cmpnge_ps);
Ret_M128_M128(__m128, _mm_cmpnge_ss);
Ret_M128_M128(__m128, _mm_cmpngt_ps);
Ret_M128_M128(__m128, _mm_cmpngt_ss);
Ret_M128_M128(__m128, _mm_cmpnle_ps);
Ret_M128_M128(__m128, _mm_cmpnle_ss);
Ret_M128_M128(__m128, _mm_cmpnlt_ps);
Ret_M128_M128(__m128, _mm_cmpnlt_ss);
Ret_M128_M128(__m128, _mm_cmpord_ps);
Ret_M128_M128(__m128, _mm_cmpord_ss);
Ret_M128_M128(__m128, _mm_cmpunord_ps);
Ret_M128_M128(__m128, _mm_cmpunord_ss);
Ret_M128_M128(int, _mm_comieq_ss);
Ret_M128_M128(int, _mm_comige_ss);
Ret_M128_M128(int, _mm_comigt_ss);
Ret_M128_M128(int, _mm_comile_ss);
Ret_M128_M128(int, _mm_comilt_ss);
Ret_M128_M128(int, _mm_comineq_ss);
Ret_M128_M128(int, _mm_ucomieq_ss);
Ret_M128_M128(int, _mm_ucomige_ss);
Ret_M128_M128(int, _mm_ucomigt_ss);
Ret_M128_M128(int, _mm_ucomile_ss);
Ret_M128_M128(int, _mm_ucomilt_ss);
Ret_M128_M128(int, _mm_ucomineq_ss);
// SSE1 Convert instructions:
Ret_M128_int(__m128, _mm_cvt_si2ss);
Ret_M128(int, _mm_cvt_ss2si);
Ret_M128_int(__m128, _mm_cvtsi32_ss);
Ret_M128(float, _mm_cvtss_f32);
Ret_M128(int, _mm_cvtss_si32);
Ret_M128(int64_t, _mm_cvtss_si64);
Ret_M128(int, _mm_cvtt_ss2si);
Ret_M128(int, _mm_cvttss_si32);
Ret_M128(int64_t, _mm_cvttss_si64);
// SSE1 Load functions:
Ret_FloatPtr(__m128, _mm_load_ps, 4, 4);
Ret_FloatPtr(__m128, _mm_load_ps1, 1, 1);
Ret_FloatPtr(__m128, _mm_load_ss, 1, 1);
Ret_FloatPtr(__m128, _mm_load1_ps, 1, 1);
Ret_M128_FloatPtr(__m128, _mm_loadh_pi, __m64*, 2, 1);
Ret_M128_FloatPtr(__m128, _mm_loadl_pi, __m64*, 2, 1);
Ret_FloatPtr(__m128, _mm_loadr_ps, 4, 4);
Ret_FloatPtr(__m128, _mm_loadu_ps, 4, 1);
// SSE1 Miscellaneous functions:
Ret_M128(int, _mm_movemask_ps);
// SSE1 Move functions:
Ret_M128_M128(__m128, _mm_move_ss);
Ret_M128_M128(__m128, _mm_movehl_ps);
Ret_M128_M128(__m128, _mm_movelh_ps);
// SSE1 Set functions:
Ret_Float4(__m128, _mm_set_ps, 1);
Ret_Float(__m128, _mm_set_ps1, 1);
Ret_Float(__m128, _mm_set_ss, 1);
Ret_Float(__m128, _mm_set1_ps, 1);
Ret_Float4(__m128, _mm_setr_ps, 1);
__m128 zero = _mm_setzero_ps();
char str[256]; tostr(&zero, str);
printf("_mm_setzero_ps() = %s\n", str);
// SSE1 Special Math instructions:
Ret_M128_M128(__m128, _mm_max_ps);
Ret_M128_M128(__m128, _mm_max_ss);
Ret_M128_M128(__m128, _mm_min_ps);
Ret_M128_M128(__m128, _mm_min_ss);
// SSE1 Store instructions:
void_OutFloatPtr_M128(_mm_store_ps, float*, 16, 16);
void_OutFloatPtr_M128(_mm_store_ps1, float*, 16, 16);
void_OutFloatPtr_M128(_mm_store_ss, float*, 4, 1);
void_OutFloatPtr_M128(_mm_store1_ps, float*, 16, 16);
void_OutFloatPtr_M128(_mm_storeh_pi, __m64*, 8, 1);
void_OutFloatPtr_M128(_mm_storel_pi, __m64*, 8, 1);
void_OutFloatPtr_M128(_mm_storer_ps, float*, 16, 16);
void_OutFloatPtr_M128(_mm_storeu_ps, float*, 16, 1);
void_OutIntPtr_M128i(_mm_storeu_si16, unsigned short*, 2, 1);
void_OutIntPtr_M128i(_mm_storeu_si64, __m64*, 8, 1);
void_OutFloatPtr_M128(_mm_stream_ps, float*, 16, 16);
// SSE1 Swizzle instructions:
Ret_M128_M128_Tint(__m128, _mm_shuffle_ps);
__m128 m1 = _mm_set_ps(1.f, 2.f, 3.f, 4.f);
__m128 m2 = _mm_set_ps(5.f, 6.f, 7.f, 8.f);
__m128 m3 = _mm_set_ps(9.f, 10.f, 11.f, 12.f);
__m128 m4 = _mm_set_ps(13.f, 14.f, 15.f, 16.f);
_MM_TRANSPOSE4_PS(m1, m2, m3, m4);
tostr(&m1, str); printf("_MM_TRANSPOSE4_PS: m1 = %s\n", str);
tostr(&m2, str); printf("_MM_TRANSPOSE4_PS: m2 = %s\n", str);
tostr(&m3, str); printf("_MM_TRANSPOSE4_PS: m3 = %s\n", str);
tostr(&m4, str); printf("_MM_TRANSPOSE4_PS: m4 = %s\n", str);
Ret_M128_M128(__m128, _mm_unpackhi_ps);
Ret_M128_M128(__m128, _mm_unpacklo_ps);
// _mm_malloc and _mm_free
void *ptr = _mm_malloc(32, 16);
assert(((uintptr_t)ptr & 0xF) == 0);
_mm_free(ptr);
#ifdef __EMSCRIPTEN__
_mm_undefined();
_mm_undefined_ps();
#endif
}
|
