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
|
/* { dg-do run } */
/* { dg-options "-O2 -msse2" } */
/* { dg-require-effective-target sse2 } */
#include "sse2-check.h"
#include <mmintrin.h>
#define N 4
unsigned long long a[N], b[N], result[N];
unsigned long long check[N] =
{ 0x101010101010100full,
0x1010101010101010ull,
0x1010101010101010ull,
0x1010101010101010ull };
__m64
unsigned_add3 (const __m64 * a, const __m64 * b,
__m64 * result, unsigned int count)
{
__m64 _a, _b, one, sum, carry, onesCarry;
unsigned int i;
carry = _mm_setzero_si64 ();
one = _mm_cmpeq_pi8 (carry, carry);
one = _mm_sub_si64 (carry, one);
for (i = 0; i < count; i++)
{
_a = a[i];
_b = b[i];
sum = _mm_add_si64 (_a, _b);
sum = _mm_add_si64 (sum, carry);
result[i] = sum;
onesCarry = _mm_and_si64 (_mm_xor_si64 (_a, _b), carry);
onesCarry = _mm_or_si64 (_mm_and_si64 (_a, _b), onesCarry);
onesCarry = _mm_and_si64 (onesCarry, one);
_a = _mm_srli_si64 (_a, 1);
_b = _mm_srli_si64 (_b, 1);
carry = _mm_add_si64 (_mm_add_si64 (_a, _b), onesCarry);
carry = _mm_srli_si64 (carry, 63);
}
return carry;
}
void __attribute__((noinline))
sse2_test (void)
{
unsigned long long carry;
int i;
/* Really long numbers. */
a[3] = a[2] = a[1] = a[0] = 0xd3d3d3d3d3d3d3d3ull;
b[3] = b[2] = b[1] = b[0] = 0x3c3c3c3c3c3c3c3cull;
carry = (unsigned long long) unsigned_add3
((__m64 *)a, (__m64 *)b, (__m64 *)result, N);
_mm_empty ();
if (carry != 1)
abort ();
for (i = 0; i < N; i++)
if (result [i] != check[i])
abort ();
}
|
