SSE Multiplikation von 4 32-bit-Ganzzahlen

Wenn Sie brauchen unterzeichnet 32x32 bit integer-Multiplikation, dann dem folgenden Beispiel an software.intel.com wie es aussieht, soll das tun, was Sie wollen:

static inline __m128i muly(const __m128i &a, const __m128i &b)
{
    __m128i tmp1 = _mm_mul_epu32(a,b); /* mul 2,0*/
    __m128i tmp2 = _mm_mul_epu32( _mm_srli_si128(a,4), _mm_srli_si128(b,4)); /* mul 3,1 */
    return _mm_unpacklo_epi32(_mm_shuffle_epi32(tmp1, _MM_SHUFFLE (0,0,2,0)), _mm_shuffle_epi32(tmp2, _MM_SHUFFLE (0,0,2,0))); /* shuffle results to [63..0] and pack */
}

Möchten Sie vielleicht haben die beiden builds - eine für die alten CPUs und einen für die aktuellen CPUs, in dem Fall können Sie das folgende tun:

static inline __m128i muly(const __m128i &a, const __m128i &b)
{
#ifdef __SSE4_1__  //modern CPU - use SSE 4.1
    return _mm_mullo_epi32(a, b);
#else               //old CPU - use SSE 2
    __m128i tmp1 = _mm_mul_epu32(a,b); /* mul 2,0*/
    __m128i tmp2 = _mm_mul_epu32( _mm_srli_si128(a,4), _mm_srli_si128(b,4)); /* mul 3,1 */
    return _mm_unpacklo_epi32(_mm_shuffle_epi32(tmp1, _MM_SHUFFLE (0,0,2,0)), _mm_shuffle_epi32(tmp2, _MM_SHUFFLE (0,0,2,0))); /* shuffle results to [63..0] and pack */
#endif
}

InformationsquelleAutor der Antwort Paul R