[x265] [PATCH] sse.inc: Replace sse_ss8 vector class function with intrinsic
dnyaneshwar at multicorewareinc.com
dnyaneshwar at multicorewareinc.com
Wed Oct 9 12:51:40 CEST 2013
# HG changeset patch
# User Dnyaneshwar Gorade <dnyaneshwar at multicorewareinc.com>
# Date 1381315864 -19800
# Wed Oct 09 16:21:04 2013 +0530
# Node ID 0e93fd14e753842677c1b899efff16cc5b52a65f
# Parent 118741a00fca214989335c483e4ddc8c6826fc0f
sse.inc: Replace sse_ss8 vector class function with intrinsic.
diff -r 118741a00fca -r 0e93fd14e753 source/common/vec/sse.inc
--- a/source/common/vec/sse.inc Wed Oct 09 15:57:33 2013 +0530
+++ b/source/common/vec/sse.inc Wed Oct 09 16:21:04 2013 +0530
@@ -78,26 +78,34 @@
int sse_ss8(short* fenc, intptr_t strideFenc, short* fref, intptr_t strideFref)
{
int rows = ly;
+ __m128i sum = _mm_setzero_si128();
+ __m128i m1, n1, diff;
- Vec4i diff(0);
- Vec8s m1, n1;
- Vec4i sum(0);
for (; rows != 0; rows--)
{
- m1.load(fenc);
- n1.load(fref);
- diff = extend_low(m1) - extend_low(n1);
- diff = diff * diff;
- sum += diff;
+ m1 = _mm_loadu_si128((__m128i const*)fenc);
+ n1 = _mm_loadu_si128((__m128i const*)fref);
- diff = extend_high(m1) - extend_high(n1);
- diff = diff * diff;
- sum += diff;
+ __m128i sign1 = _mm_srai_epi16(m1, 15); // sign bit
+ __m128i tmp1 = _mm_unpacklo_epi16(m1, sign1); // interleave with sign extensions
+ __m128i sign2 = _mm_srai_epi16(n1, 15);
+ __m128i tmp2 = _mm_unpacklo_epi16(n1, sign2);
+ diff = _mm_sub_epi32(tmp1, tmp2);
+ diff = _mm_mullo_epi32(diff, diff);
+ sum = _mm_add_epi32(sum, diff);
+
+ m1 = _mm_unpackhi_epi16(m1, sign1);
+ n1 = _mm_unpackhi_epi16(n1, sign2);
+
+ diff = _mm_sub_epi32(m1, n1);
+ diff = _mm_mullo_epi32(diff, diff);
+ sum = _mm_add_epi32(sum, diff);
+
fenc += strideFenc;
fref += strideFref;
}
-
- return horizontal_add(sum);
+ __m128i sum1 = _mm_hadd_epi32(sum, sum); // horizontally add 4 elements
+ return _mm_cvtsi128_si32(_mm_hadd_epi32(sum1, sum1));
}
template<int ly>
More information about the x265-devel
mailing list