SSE moltiplicazione di 2 interi a 64 bit

Question

Come moltiplicare due numeri interi a 64 bit da altri 2 numeri interi a 64 bit? Non ho trovato alcuna istruzione che possa farlo.

Answer 1

È necessario implementare la propria routine di moltiplicazione a 64 bit utilizzando le operazioni di moltiplicazione a 32 bit. Probabilmente non sarà più efficiente del semplice fare questo con il codice scalare, soprattutto perché ci sarà un sacco di mischiare i vettori per ottenere tutte le operazioni necessarie.

Answer 2

So che questa è una domanda vecchia ma in realtà stavo cercando esattamente questo. Dato che non ci sono ancora istruzioni per farlo, ho implementato il 64 bit per moltiplicare me stesso con il pmuldq come menzionato da Paul R. Questo è quello che è venuto fuori con

__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)); 
}