Es una aplicación bastante sencilla Aquí (advertencia: código no probado):
int32_t sum_array(const int32_t a[], const int n)
{
__m128i vsum = _mm_set1_epi32(0); // initialise vector of four partial 32 bit sums
int32_t sum;
int i;
for (i = 0; i < n; i += 4)
{
__m128i v = _mm_load_si128(&a[i]); // load vector of 4 x 32 bit values
vsum = _mm_add_epi32(vsum, v); // accumulate to 32 bit partial sum vector
}
// horizontal add of four 32 bit partial sums and return result
vsum = _mm_add_epi32(vsum, _mm_srli_si128(vsum, 8));
vsum = _mm_add_epi32(vsum, _mm_srli_si128(vsum, 4));
sum = _mm_cvtsi128_si32(vsum);
return sum;
}
Tenga en cuenta que la matriz de entrada, a[]
, es necesario que estar alineado 16 bytes, y n
debe ser un múltiplo de 4.
¿Qué SIMD tenía en mente? SSE2? – harold
SSE Las siguientes características intrínsecas podrían utilizarse .__ _mm_setzero_si128 m128i() __m128i _mm_loadu_si128 (__m128i * p) __m128i _mm_add_epi32 (__m128i a, __m128i b) (a0 + B0, a1 + b1, a2 + b2, a3 + b3) void _mm_storeu_si128 (__m128i * p, __m128i a) – user1585869
Bien, entonces SSE2. ¿Qué has intentado? – harold