[x265] [PATCH] pixel8.inc: sad_x4_8, vector replaced with intrinsic
praveen at multicorewareinc.com
praveen at multicorewareinc.com
Mon Aug 26 10:48:42 CEST 2013
# HG changeset patch
# User praveentiwari
# Date 1377506902 -19800
# Node ID f2892688afb2f52de6821b1865e34e7e674bf01d
# Parent 77a45df4b06947c4d7c17fe97390f25548c43fbf
pixel8.inc: sad_x4_8, vector replaced with intrinsic
diff -r 77a45df4b069 -r f2892688afb2 source/common/vec/pixel8.inc
--- a/source/common/vec/pixel8.inc Mon Aug 26 13:09:20 2013 +0530
+++ b/source/common/vec/pixel8.inc Mon Aug 26 14:18:22 2013 +0530
@@ -5372,81 +5372,1323 @@
}
#endif /* if HAVE_MMX */
-#endif /* if INSTRSET >= X265_CPU_LEVEL_SSE41 */
-
+
+#if HAVE_MMX
template<int ly>
void sad_x4_8(pixel *fenc, pixel *fref1, pixel *fref2, pixel *fref3, pixel *fref4, intptr_t frefstride, int *res)
{
- Vec16uc m1, n1, n2, n3, n4;
-
- Vec4i sum1(0), sum2(0), sum3(0), sum4(0);
- Vec8us sad1(0), sad2(0), sad3(0), sad4(0);
- int max_iterators = (ly >> 4) << 4;
- int row;
-
- for (row = 0; row < max_iterators; row += 16)
+ assert((ly % 4) == 0);
+
+ __m64 sum0 = _mm_setzero_si64();
+ __m64 sum1 = _mm_setzero_si64();
+ __m64 sum2 = _mm_setzero_si64();
+ __m64 sum3 = _mm_setzero_si64();
+
+ __m64 T00, T01, T02, T03, T04, T05, T06, T07;
+ __m64 T0, T1, T2, T3, T4, T5, T6, T7;
+ __m64 T10, T11, T12, T13, T14, T15, T16, T17;
+ __m64 T20, T21, T22, T23, T24, T25, T26, T27;
+
+ if (4 == ly)
{
- for (int i = 0; i < 16; i++)
+ T00 = (*(__m64*)(fenc + 0 * FENC_STRIDE));
+ T01 = (*(__m64*)(fenc + 1 * FENC_STRIDE));
+ T02 = (*(__m64*)(fenc + 2 * FENC_STRIDE));
+ T03 = (*(__m64*)(fenc + 3 * FENC_STRIDE));
+
+ T10 = (*(__m64*)(fref1 + 0 * frefstride));
+ T11 = (*(__m64*)(fref1 + 1 * frefstride));
+ T12 = (*(__m64*)(fref1 + 2 * frefstride));
+ T13 = (*(__m64*)(fref1 + 3 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ res[0] = _m_to_int(sum0);
+
+ T10 = (*(__m64*)(fref2 + 0 * frefstride));
+ T11 = (*(__m64*)(fref2 + 1 * frefstride));
+ T12 = (*(__m64*)(fref2 + 2 * frefstride));
+ T13 = (*(__m64*)(fref2 + 3 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ res[1] = _m_to_int(sum1);
+
+ T10 = (*(__m64*)(fref3 + 0 * frefstride));
+ T11 = (*(__m64*)(fref3 + 1 * frefstride));
+ T12 = (*(__m64*)(fref3 + 2 * frefstride));
+ T13 = (*(__m64*)(fref3 + 3 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ res[2] = _m_to_int(sum2);
+
+ T10 = (*(__m64*)(fref4 + 0 * frefstride));
+ T11 = (*(__m64*)(fref4 + 1 * frefstride));
+ T12 = (*(__m64*)(fref4 + 2 * frefstride));
+ T13 = (*(__m64*)(fref4 + 3 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ res[3] = _m_to_int(sum3);
+ }
+ else if (8 == ly)
+ {
+ T00 = (*(__m64*)(fenc + 0 * FENC_STRIDE));
+ T01 = (*(__m64*)(fenc + 1 * FENC_STRIDE));
+ T02 = (*(__m64*)(fenc + 2 * FENC_STRIDE));
+ T03 = (*(__m64*)(fenc + 3 * FENC_STRIDE));
+ T04 = (*(__m64*)(fenc + 4 * FENC_STRIDE));
+ T05 = (*(__m64*)(fenc + 5 * FENC_STRIDE));
+ T06 = (*(__m64*)(fenc + 6 * FENC_STRIDE));
+ T07 = (*(__m64*)(fenc + 7 * FENC_STRIDE));
+
+ T10 = (*(__m64*)(fref1 + 0 * frefstride));
+ T11 = (*(__m64*)(fref1 + 1 * frefstride));
+ T12 = (*(__m64*)(fref1 + 2 * frefstride));
+ T13 = (*(__m64*)(fref1 + 3 * frefstride));
+ T14 = (*(__m64*)(fref1 + 4 * frefstride));
+ T15 = (*(__m64*)(fref1 + 5 * frefstride));
+ T16 = (*(__m64*)(fref1 + 6 * frefstride));
+ T17 = (*(__m64*)(fref1 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ sum0 = _mm_add_pi16(sum0, T24);
+ sum0 = _mm_add_pi16(sum0, T25);
+ sum0 = _mm_add_pi16(sum0, T26);
+ sum0 = _mm_add_pi16(sum0, T27);
+ res[0] = _m_to_int(sum0);
+
+ T10 = (*(__m64*)(fref2 + 0 * frefstride));
+ T11 = (*(__m64*)(fref2 + 1 * frefstride));
+ T12 = (*(__m64*)(fref2 + 2 * frefstride));
+ T13 = (*(__m64*)(fref2 + 3 * frefstride));
+ T14 = (*(__m64*)(fref2 + 4 * frefstride));
+ T15 = (*(__m64*)(fref2 + 5 * frefstride));
+ T16 = (*(__m64*)(fref2 + 6 * frefstride));
+ T17 = (*(__m64*)(fref2 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ sum1 = _mm_add_pi16(sum1, T24);
+ sum1 = _mm_add_pi16(sum1, T25);
+ sum1 = _mm_add_pi16(sum1, T26);
+ sum1 = _mm_add_pi16(sum1, T27);
+ res[1] = _m_to_int(sum1);
+
+ T10 = (*(__m64*)(fref3 + 0 * frefstride));
+ T11 = (*(__m64*)(fref3 + 1 * frefstride));
+ T12 = (*(__m64*)(fref3 + 2 * frefstride));
+ T13 = (*(__m64*)(fref3 + 3 * frefstride));
+ T14 = (*(__m64*)(fref3 + 4 * frefstride));
+ T15 = (*(__m64*)(fref3 + 5 * frefstride));
+ T16 = (*(__m64*)(fref3 + 6 * frefstride));
+ T17 = (*(__m64*)(fref3 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ sum2 = _mm_add_pi16(sum2, T24);
+ sum2 = _mm_add_pi16(sum2, T25);
+ sum2 = _mm_add_pi16(sum2, T26);
+ sum2 = _mm_add_pi16(sum2, T27);
+ res[2] = _m_to_int(sum2);
+
+ T10 = (*(__m64*)(fref4 + 0 * frefstride));
+ T11 = (*(__m64*)(fref4 + 1 * frefstride));
+ T12 = (*(__m64*)(fref4 + 2 * frefstride));
+ T13 = (*(__m64*)(fref4 + 3 * frefstride));
+ T14 = (*(__m64*)(fref4 + 4 * frefstride));
+ T15 = (*(__m64*)(fref4 + 5 * frefstride));
+ T16 = (*(__m64*)(fref4 + 6 * frefstride));
+ T17 = (*(__m64*)(fref4 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ sum3 = _mm_add_pi16(sum3, T24);
+ sum3 = _mm_add_pi16(sum3, T25);
+ sum3 = _mm_add_pi16(sum3, T26);
+ sum3 = _mm_add_pi16(sum3, T27);
+ res[3] = _m_to_int(sum3);
+ }
+ else if (16 == ly)
+ {
+ T00 = (*(__m64*)(fenc + 0 * FENC_STRIDE));
+ T01 = (*(__m64*)(fenc + 1 * FENC_STRIDE));
+ T02 = (*(__m64*)(fenc + 2 * FENC_STRIDE));
+ T03 = (*(__m64*)(fenc + 3 * FENC_STRIDE));
+ T04 = (*(__m64*)(fenc + 4 * FENC_STRIDE));
+ T05 = (*(__m64*)(fenc + 5 * FENC_STRIDE));
+ T06 = (*(__m64*)(fenc + 6 * FENC_STRIDE));
+ T07 = (*(__m64*)(fenc + 7 * FENC_STRIDE));
+ T0 = (*(__m64*)(fenc + 8 * FENC_STRIDE));
+ T1 = (*(__m64*)(fenc + 9 * FENC_STRIDE));
+ T2 = (*(__m64*)(fenc + 10 * FENC_STRIDE));
+ T3 = (*(__m64*)(fenc + 11 * FENC_STRIDE));
+ T4 = (*(__m64*)(fenc + 12 * FENC_STRIDE));
+ T5 = (*(__m64*)(fenc + 13 * FENC_STRIDE));
+ T6 = (*(__m64*)(fenc + 14 * FENC_STRIDE));
+ T7 = (*(__m64*)(fenc + 15 * FENC_STRIDE));
+
+ T10 = (*(__m64*)(fref1 + 0 * frefstride));
+ T11 = (*(__m64*)(fref1 + 1 * frefstride));
+ T12 = (*(__m64*)(fref1 + 2 * frefstride));
+ T13 = (*(__m64*)(fref1 + 3 * frefstride));
+ T14 = (*(__m64*)(fref1 + 4 * frefstride));
+ T15 = (*(__m64*)(fref1 + 5 * frefstride));
+ T16 = (*(__m64*)(fref1 + 6 * frefstride));
+ T17 = (*(__m64*)(fref1 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ sum0 = _mm_add_pi16(sum0, T24);
+ sum0 = _mm_add_pi16(sum0, T25);
+ sum0 = _mm_add_pi16(sum0, T26);
+ sum0 = _mm_add_pi16(sum0, T27);
+
+ T10 = (*(__m64*)(fref1 + 8 * frefstride));
+ T11 = (*(__m64*)(fref1 + 9 * frefstride));
+ T12 = (*(__m64*)(fref1 + 10 * frefstride));
+ T13 = (*(__m64*)(fref1 + 11 * frefstride));
+ T14 = (*(__m64*)(fref1 + 12 * frefstride));
+ T15 = (*(__m64*)(fref1 + 13 * frefstride));
+ T16 = (*(__m64*)(fref1 + 14 * frefstride));
+ T17 = (*(__m64*)(fref1 + 15 * frefstride));
+
+ T20 = _mm_sad_pu8(T0, T10);
+ T21 = _mm_sad_pu8(T1, T11);
+ T22 = _mm_sad_pu8(T2, T12);
+ T23 = _mm_sad_pu8(T3, T13);
+ T24 = _mm_sad_pu8(T4, T14);
+ T25 = _mm_sad_pu8(T5, T15);
+ T26 = _mm_sad_pu8(T6, T16);
+ T27 = _mm_sad_pu8(T7, T17);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ sum0 = _mm_add_pi16(sum0, T24);
+ sum0 = _mm_add_pi16(sum0, T25);
+ sum0 = _mm_add_pi16(sum0, T26);
+ sum0 = _mm_add_pi16(sum0, T27);
+ res[0] = _m_to_int(sum0);
+
+ T10 = (*(__m64*)(fref2 + 0 * frefstride));
+ T11 = (*(__m64*)(fref2 + 1 * frefstride));
+ T12 = (*(__m64*)(fref2 + 2 * frefstride));
+ T13 = (*(__m64*)(fref2 + 3 * frefstride));
+ T14 = (*(__m64*)(fref2 + 4 * frefstride));
+ T15 = (*(__m64*)(fref2 + 5 * frefstride));
+ T16 = (*(__m64*)(fref2 + 6 * frefstride));
+ T17 = (*(__m64*)(fref2 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ sum1 = _mm_add_pi16(sum1, T24);
+ sum1 = _mm_add_pi16(sum1, T25);
+ sum1 = _mm_add_pi16(sum1, T26);
+ sum1 = _mm_add_pi16(sum1, T27);
+
+ T10 = (*(__m64*)(fref2 + 8 * frefstride));
+ T11 = (*(__m64*)(fref2 + 9 * frefstride));
+ T12 = (*(__m64*)(fref2 + 10 * frefstride));
+ T13 = (*(__m64*)(fref2 + 11 * frefstride));
+ T14 = (*(__m64*)(fref2 + 12 * frefstride));
+ T15 = (*(__m64*)(fref2 + 13 * frefstride));
+ T16 = (*(__m64*)(fref2 + 14 * frefstride));
+ T17 = (*(__m64*)(fref2 + 15 * frefstride));
+
+ T20 = _mm_sad_pu8(T0, T10);
+ T21 = _mm_sad_pu8(T1, T11);
+ T22 = _mm_sad_pu8(T2, T12);
+ T23 = _mm_sad_pu8(T3, T13);
+ T24 = _mm_sad_pu8(T4, T14);
+ T25 = _mm_sad_pu8(T5, T15);
+ T26 = _mm_sad_pu8(T6, T16);
+ T27 = _mm_sad_pu8(T7, T17);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ sum1 = _mm_add_pi16(sum1, T24);
+ sum1 = _mm_add_pi16(sum1, T25);
+ sum1 = _mm_add_pi16(sum1, T26);
+ sum1 = _mm_add_pi16(sum1, T27);
+ res[1] = _m_to_int(sum1);
+
+ T10 = (*(__m64*)(fref3 + 0 * frefstride));
+ T11 = (*(__m64*)(fref3 + 1 * frefstride));
+ T12 = (*(__m64*)(fref3 + 2 * frefstride));
+ T13 = (*(__m64*)(fref3 + 3 * frefstride));
+ T14 = (*(__m64*)(fref3 + 4 * frefstride));
+ T15 = (*(__m64*)(fref3 + 5 * frefstride));
+ T16 = (*(__m64*)(fref3 + 6 * frefstride));
+ T17 = (*(__m64*)(fref3 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ sum2 = _mm_add_pi16(sum2, T24);
+ sum2 = _mm_add_pi16(sum2, T25);
+ sum2 = _mm_add_pi16(sum2, T26);
+ sum2 = _mm_add_pi16(sum2, T27);
+
+ T10 = (*(__m64*)(fref3 + 8 * frefstride));
+ T11 = (*(__m64*)(fref3 + 9 * frefstride));
+ T12 = (*(__m64*)(fref3 + 10 * frefstride));
+ T13 = (*(__m64*)(fref3 + 11 * frefstride));
+ T14 = (*(__m64*)(fref3 + 12 * frefstride));
+ T15 = (*(__m64*)(fref3 + 13 * frefstride));
+ T16 = (*(__m64*)(fref3 + 14 * frefstride));
+ T17 = (*(__m64*)(fref3 + 15 * frefstride));
+
+ T20 = _mm_sad_pu8(T0, T10);
+ T21 = _mm_sad_pu8(T1, T11);
+ T22 = _mm_sad_pu8(T2, T12);
+ T23 = _mm_sad_pu8(T3, T13);
+ T24 = _mm_sad_pu8(T4, T14);
+ T25 = _mm_sad_pu8(T5, T15);
+ T26 = _mm_sad_pu8(T6, T16);
+ T27 = _mm_sad_pu8(T7, T17);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ sum2 = _mm_add_pi16(sum2, T24);
+ sum2 = _mm_add_pi16(sum2, T25);
+ sum2 = _mm_add_pi16(sum2, T26);
+ sum2 = _mm_add_pi16(sum2, T27);
+ res[2] = _m_to_int(sum2);
+
+ T10 = (*(__m64*)(fref4 + 0 * frefstride));
+ T11 = (*(__m64*)(fref4 + 1 * frefstride));
+ T12 = (*(__m64*)(fref4 + 2 * frefstride));
+ T13 = (*(__m64*)(fref4 + 3 * frefstride));
+ T14 = (*(__m64*)(fref4 + 4 * frefstride));
+ T15 = (*(__m64*)(fref4 + 5 * frefstride));
+ T16 = (*(__m64*)(fref4 + 6 * frefstride));
+ T17 = (*(__m64*)(fref4 + 7 * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ sum3 = _mm_add_pi16(sum3, T24);
+ sum3 = _mm_add_pi16(sum3, T25);
+ sum3 = _mm_add_pi16(sum3, T26);
+ sum3 = _mm_add_pi16(sum3, T27);
+
+ T10 = (*(__m64*)(fref4 + 8 * frefstride));
+ T11 = (*(__m64*)(fref4 + 9 * frefstride));
+ T12 = (*(__m64*)(fref4 + 10 * frefstride));
+ T13 = (*(__m64*)(fref4 + 11 * frefstride));
+ T14 = (*(__m64*)(fref4 + 12 * frefstride));
+ T15 = (*(__m64*)(fref4 + 13 * frefstride));
+ T16 = (*(__m64*)(fref4 + 14 * frefstride));
+ T17 = (*(__m64*)(fref4 + 15 * frefstride));
+
+ T20 = _mm_sad_pu8(T0, T10);
+ T21 = _mm_sad_pu8(T1, T11);
+ T22 = _mm_sad_pu8(T2, T12);
+ T23 = _mm_sad_pu8(T3, T13);
+ T24 = _mm_sad_pu8(T4, T14);
+ T25 = _mm_sad_pu8(T5, T15);
+ T26 = _mm_sad_pu8(T6, T16);
+ T27 = _mm_sad_pu8(T7, T17);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ sum3 = _mm_add_pi16(sum3, T24);
+ sum3 = _mm_add_pi16(sum3, T25);
+ sum3 = _mm_add_pi16(sum3, T26);
+ sum3 = _mm_add_pi16(sum3, T27);
+ res[3] = _m_to_int(sum3);
+ }
+ else if (!(ly % 8))
+ {
+ for (int i = 0; i < ly; i += 8)
{
- m1.load_a(fenc);
- n1.load(fref1);
- n2.load(fref2);
- n3.load(fref3);
- n4.load(fref4);
-
- sad1.addSumAbsDiff(m1, n1);
- sad2.addSumAbsDiff(m1, n2);
- sad3.addSumAbsDiff(m1, n3);
- sad4.addSumAbsDiff(m1, n4);
-
- fenc += FENC_STRIDE;
- fref1 += frefstride;
- fref2 += frefstride;
- fref3 += frefstride;
- fref4 += frefstride;
+ T00 = (*(__m64*)(fenc + (i + 0) * FENC_STRIDE));
+ T01 = (*(__m64*)(fenc + (i + 1) * FENC_STRIDE));
+ T02 = (*(__m64*)(fenc + (i + 2) * FENC_STRIDE));
+ T03 = (*(__m64*)(fenc + (i + 3) * FENC_STRIDE));
+ T04 = (*(__m64*)(fenc + (i + 4) * FENC_STRIDE));
+ T05 = (*(__m64*)(fenc + (i + 5) * FENC_STRIDE));
+ T06 = (*(__m64*)(fenc + (i + 6) * FENC_STRIDE));
+ T07 = (*(__m64*)(fenc + (i + 7) * FENC_STRIDE));
+
+ T10 = (*(__m64*)(fref1 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref1 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref1 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref1 + (i + 3) * frefstride));
+ T14 = (*(__m64*)(fref1 + (i + 4) * frefstride));
+ T15 = (*(__m64*)(fref1 + (i + 5) * frefstride));
+ T16 = (*(__m64*)(fref1 + (i + 6) * frefstride));
+ T17 = (*(__m64*)(fref1 + (i + 7) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ sum0 = _mm_add_pi16(sum0, T24);
+ sum0 = _mm_add_pi16(sum0, T25);
+ sum0 = _mm_add_pi16(sum0, T26);
+ sum0 = _mm_add_pi16(sum0, T27);
+ res[0] = _m_to_int(sum0);
+
+ T10 = (*(__m64*)(fref2 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref2 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref2 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref2 + (i + 3) * frefstride));
+ T14 = (*(__m64*)(fref2 + (i + 4) * frefstride));
+ T15 = (*(__m64*)(fref2 + (i + 5) * frefstride));
+ T16 = (*(__m64*)(fref2 + (i + 6) * frefstride));
+ T17 = (*(__m64*)(fref2 + (i + 7) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ sum1 = _mm_add_pi16(sum1, T24);
+ sum1 = _mm_add_pi16(sum1, T25);
+ sum1 = _mm_add_pi16(sum1, T26);
+ sum1 = _mm_add_pi16(sum1, T27);
+ res[1] = _m_to_int(sum1);
+
+ T10 = (*(__m64*)(fref3 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref3 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref3 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref3 + (i + 3) * frefstride));
+ T14 = (*(__m64*)(fref3 + (i + 4) * frefstride));
+ T15 = (*(__m64*)(fref3 + (i + 5) * frefstride));
+ T16 = (*(__m64*)(fref3 + (i + 6) * frefstride));
+ T17 = (*(__m64*)(fref3 + (i + 7) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ sum2 = _mm_add_pi16(sum2, T24);
+ sum2 = _mm_add_pi16(sum2, T25);
+ sum2 = _mm_add_pi16(sum2, T26);
+ sum2 = _mm_add_pi16(sum2, T27);
+ res[2] = _m_to_int(sum2);
+
+ T10 = (*(__m64*)(fref4 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref4 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref4 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref4 + (i + 3) * frefstride));
+ T14 = (*(__m64*)(fref4 + (i + 4) * frefstride));
+ T15 = (*(__m64*)(fref4 + (i + 5) * frefstride));
+ T16 = (*(__m64*)(fref4 + (i + 6) * frefstride));
+ T17 = (*(__m64*)(fref4 + (i + 7) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+ T24 = _mm_sad_pu8(T04, T14);
+ T25 = _mm_sad_pu8(T05, T15);
+ T26 = _mm_sad_pu8(T06, T16);
+ T27 = _mm_sad_pu8(T07, T17);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ sum3 = _mm_add_pi16(sum3, T24);
+ sum3 = _mm_add_pi16(sum3, T25);
+ sum3 = _mm_add_pi16(sum3, T26);
+ sum3 = _mm_add_pi16(sum3, T27);
+ res[3] = _m_to_int(sum3);
}
-
- sum1 += extend_low(sad1);
- sum2 += extend_low(sad2);
- sum3 += extend_low(sad3);
- sum4 += extend_low(sad4);
- sad1 = 0;
- sad2 = 0;
- sad3 = 0;
- sad4 = 0;
}
-
- while (row++ < ly)
+ else
{
- m1.load_a(fenc);
- n1.load(fref1);
- n2.load(fref2);
- n3.load(fref3);
- n4.load(fref4);
-
- sad1.addSumAbsDiff(m1, n1);
- sad2.addSumAbsDiff(m1, n2);
- sad3.addSumAbsDiff(m1, n3);
- sad4.addSumAbsDiff(m1, n4);
-
- fenc += FENC_STRIDE;
- fref1 += frefstride;
- fref2 += frefstride;
- fref3 += frefstride;
- fref4 += frefstride;
+ for (int i = 0; i < ly; i += 4)
+ {
+ T00 = (*(__m64*)(fenc + (i + 0) * FENC_STRIDE));
+ T01 = (*(__m64*)(fenc + (i + 1) * FENC_STRIDE));
+ T02 = (*(__m64*)(fenc + (i + 2) * FENC_STRIDE));
+ T03 = (*(__m64*)(fenc + (i + 3) * FENC_STRIDE));
+
+ T10 = (*(__m64*)(fref1 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref1 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref1 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref1 + (i + 3) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum0 = _mm_add_pi16(sum0, T20);
+ sum0 = _mm_add_pi16(sum0, T21);
+ sum0 = _mm_add_pi16(sum0, T22);
+ sum0 = _mm_add_pi16(sum0, T23);
+ res[0] = _m_to_int(sum0);
+
+ T10 = (*(__m64*)(fref2 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref2 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref2 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref2 + (i + 3) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum1 = _mm_add_pi16(sum1, T20);
+ sum1 = _mm_add_pi16(sum1, T21);
+ sum1 = _mm_add_pi16(sum1, T22);
+ sum1 = _mm_add_pi16(sum1, T23);
+ res[1] = _m_to_int(sum1);
+
+ T10 = (*(__m64*)(fref3 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref3 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref3 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref3 + (i + 3) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum2 = _mm_add_pi16(sum2, T20);
+ sum2 = _mm_add_pi16(sum2, T21);
+ sum2 = _mm_add_pi16(sum2, T22);
+ sum2 = _mm_add_pi16(sum2, T23);
+ res[2] = _m_to_int(sum2);
+
+ T10 = (*(__m64*)(fref4 + (i + 0) * frefstride));
+ T11 = (*(__m64*)(fref4 + (i + 1) * frefstride));
+ T12 = (*(__m64*)(fref4 + (i + 2) * frefstride));
+ T13 = (*(__m64*)(fref4 + (i + 3) * frefstride));
+
+ T20 = _mm_sad_pu8(T00, T10);
+ T21 = _mm_sad_pu8(T01, T11);
+ T22 = _mm_sad_pu8(T02, T12);
+ T23 = _mm_sad_pu8(T03, T13);
+
+ sum3 = _mm_add_pi16(sum3, T20);
+ sum3 = _mm_add_pi16(sum3, T21);
+ sum3 = _mm_add_pi16(sum3, T22);
+ sum3 = _mm_add_pi16(sum3, T23);
+ res[3] = _m_to_int(sum3);
+ }
}
-
- sum1 += extend_low(sad1);
- sum2 += extend_low(sad2);
- sum3 += extend_low(sad3);
- sum4 += extend_low(sad4);
-
- res[0] = horizontal_add(sum1);
- res[1] = horizontal_add(sum2);
- res[2] = horizontal_add(sum3);
- res[3] = horizontal_add(sum4);
}
+#else /* if HAVE_MMX */
+
+template<int ly>
+void sad_x4_8(pixel *fenc, pixel *fref1, pixel *fref2, pixel *fref3, pixel *fref4, intptr_t frefstride, int *res)
+{
+ assert((ly % 4) == 0);
+ __m128i sum0 = _mm_setzero_si128();
+
+ __m128i T00, T01, T02, T03;
+ __m128i T10, T11, T12, T13;
+ __m128i T20, T21;
+
+ if (ly == 4)
+ {
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (0) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (1) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (2) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (3) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = _mm_cvtsi128_si32(sum0);
+ }
+ else if (ly == 8)
+ {
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (0) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (1) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (2) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (3) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = _mm_cvtsi128_si32(sum0);
+
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (4) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (5) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (6) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (7) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+ }
+ else if (ly == 16)
+ {
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (0) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (1) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (2) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (3) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = _mm_cvtsi128_si32(sum0);
+
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (4) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (5) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (6) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (7) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (8) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (9) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (10) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (11) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (8) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (9) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (10) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (11) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (8) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (9) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (10) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (11) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (8) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (9) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (10) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (11) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (8) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (9) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (10) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (11) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (12) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (13) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (14) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (15) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (12) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (13) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (14) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (15) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (12) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (13) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (14) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (15) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (12) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (13) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (14) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (15) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (12) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (13) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (14) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (15) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+ }
+ else if ((ly % 8) == 0)
+ {
+ res[0] = res[1] = res[2] = res[3] = 0;
+ for (int i = 0; i < ly; i += 8)
+ {
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (i + 0) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (i + 1) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (i + 2) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (i + 3) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (i + 4) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (i + 5) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (i + 6) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (i + 7) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 4) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 5) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 6) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 7) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+ }
+ }
+ else
+ {
+ res[0] = res[1] = res[2] = res[3] = 0;
+ for (int i = 0; i < ly; i += 4)
+ {
+ T00 = _mm_loadl_epi64((__m128i*)(fenc + (i + 0) * FENC_STRIDE));
+ T01 = _mm_loadl_epi64((__m128i*)(fenc + (i + 1) * FENC_STRIDE));
+ T01 = _mm_unpacklo_epi64(T00, T01);
+ T02 = _mm_loadl_epi64((__m128i*)(fenc + (i + 2) * FENC_STRIDE));
+ T03 = _mm_loadl_epi64((__m128i*)(fenc + (i + 3) * FENC_STRIDE));
+ T03 = _mm_unpacklo_epi64(T02, T03);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref1 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[0] = res[0] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref2 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[1] = res[1] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref3 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[2] = res[2] + _mm_cvtsi128_si32(sum0);
+
+ T10 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 0) * frefstride));
+ T11 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 1) * frefstride));
+ T11 = _mm_unpacklo_epi64(T10, T11);
+ T12 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 2) * frefstride));
+ T13 = _mm_loadl_epi64((__m128i*)(fref4 + (i + 3) * frefstride));
+ T13 = _mm_unpacklo_epi64(T12, T13);
+
+ T20 = _mm_sad_epu8(T01, T11);
+ T21 = _mm_sad_epu8(T03, T13);
+ T21 = _mm_add_epi32(T20, T21);
+ sum0 = _mm_shuffle_epi32(T21, 2);
+ sum0 = _mm_add_epi32(sum0, T21);
+ res[3] = res[3] + _mm_cvtsi128_si32(sum0);
+ }
+ }
+}
+
+#endif /* if HAVE_MMX */
+#endif /* if INSTRSET >= X265_CPU_LEVEL_SSE41 */
+
/* For performance - This function assumes that the *last load* can access 16 elements. */
template<int ly>
void sad_x4_12(pixel *fenc, pixel *fref1, pixel *fref2, pixel *fref3, pixel *fref4, intptr_t frefstride, int *res)
More information about the x265-devel
mailing list