SSE mnożenie z 2 64-bitowych liczb całkowitych

Question

Jak pomnożyć dwa 64-bitowe liczby całkowite przez kolejne 2 64-bitowe liczby całkowite? Nie znalazłem instrukcji, które mogą to zrobić.

Answer 1

Trzeba by zaimplementować własną 64 bitową procedurę mnożenia za pomocą 32-bitowych operacji mnożyć. Prawdopodobnie nie będzie to bardziej efektywne niż wykonanie tego za pomocą kodu skalarnego, szczególnie, że będzie dużo tasowania wektorów, aby uzyskać wszystkie wymagane operacje.

Answer 2

wiem, że to jest stare pytanie, ale ja rzeczywiście szukają dokładnie tego. Ponieważ wciąż nie ma instrukcji do tego, zaimplementowałem 64-bitowe mnożenie się z pmuldq jak wspomniał Paul R. To właśnie wymyśliłem:

__m128i Multiply64Bit(__m128i a, __m128i b) 
{ 
    auto ax0_ax1_ay0_ay1 = a; 
    auto bx0_bx1_by0_by1 = b; 

    // i means ignored 

    auto ax1_i_ay1_i = _mm_shuffle_epi32(ax0_ax1_ay0_ay1, _MM_SHUFFLE(3, 3, 1, 1)); 
    auto bx1_i_by1_i = _mm_shuffle_epi32(bx0_bx1_by0_by1, _MM_SHUFFLE(3, 3, 1, 1)); 

    auto ax0bx0_ay0by0 = _mm_mul_epi32(ax0_ax1_ay0_ay1, bx0_bx1_by0_by1); 
    auto ax0bx1_ay0by1 = _mm_mul_epi32(ax0_ax1_ay0_ay1, bx1_i_by1_i); 
    auto ax1bx0_ay1by0 = _mm_mul_epi32(ax1_i_ay1_i, bx0_bx1_by0_by1); 

    auto ax0bx1_ay0by1_32 = _mm_slli_epi64(ax0bx1_ay0by1, 32); 
    auto ax1bx0_ay1by0_32 = _mm_slli_epi64(ax1bx0_ay1by0, 32); 

    return _mm_add_epi64(ax0bx0_ay0by0, _mm_add_epi64(ax0bx1_ay0by1_32, ax1bx0_ay1by0_32)); 
}