这是另一种SSE is slower than normal code! Why?类型的问题。
我知道有很多类似的问题，但它们似乎与我的情况不符。

我正在尝试使用蒙哥马利模乘法实现Miller-Rabin 素性测试以进行快速模运算。
我尝试以标量和 SIMD 方式实现它，结果发现 SIMD 版本慢了大约 10%。
[esp+16] 或 [esp+12] 指向是否有人想知道的模倒数。

N

我真的很困惑这样一个事实，即所谓的 1 Latency 1c Throughput 1uops 指令psrldq需要超过 3 Latency 0.5c Throughput 1uops pmuludq。

下面是在 Ryzen 5 3600 上运行的 Visual Studio 的代码和运行时分析。

关于如何改进 SIMD 代码和/或为什么它比标量代码慢的任何想法表示赞赏。

PS似乎由于某种原因，一条指令关闭了运行时分析

编辑1：对图像的评论是错误的，我在下面附上了一个固定版本：

    ;------------------------------------------
    ; Calculate base^d mod x
    ;
    ; eax = 1
    ; esi = x
    ; edi = bases[eax]
    ; ebp = d
    ; while d do
    ;     if d & 1 then eax = (eax * edi) mod x
    ;     edi = (edi*edi) mod x
    ;     d >>= 1
    ; end
    ;------------------------------------------

标量代码：

LOOP_MODEXP:
    push eax

    test ebp, 1
    jz @F

    mul edi
    mov ecx, edx
    imul eax, DWORD PTR [esp+16]
    mul esi
    xor ebx, ebx
    sub ecx, edx
    cmovs ebx, esi
    add ecx, ebx
    mov DWORD PTR [esp], ecx
@@:
    mov edx, edi
    mulx ecx, edx, edi
    imul edx, DWORD PTR [esp+16]
    mulx eax, ebx, esi
    xor ebx, ebx
    sub ecx, eax
    cmovs ebx, esi
    add ecx, ebx
    mov edi, ecx

    pop eax

    shr ebp, 1
    jnz LOOP_MODEXP

SIMD 代码

    movd xmm2, DWORD PTR [esp+12]
    movd xmm3, esi
    pshufd xmm2, xmm2, 0
    pshufd xmm3, xmm3, 0
    
    movd xmm1, edi

    pshufd xmm1, xmm1, 0
    movdqa xmm0, xmm1

    pinsrd xmm0, eax, 2

LOOP_MODEXP:
    movdqa xmm4, xmm0
    pmuludq xmm0, xmm1
    movdqa xmm1, xmm0
    pmuludq xmm0, xmm2
    pmuludq xmm0, xmm3
    psubd xmm1, xmm0
    
    psrldq xmm1, 4
    pxor xmm0, xmm0
    pcmpgtd xmm0, xmm1
    blendvps xmm0, xmm3, xmm0
    paddd xmm0, xmm1

    movddup xmm1, xmm0

    test ebp, 1
    jnz @F
    blendps xmm0, xmm4, 4

@@:
    shr ebp, 1
    jnz LOOP_MODEXP

    pextrd eax, xmm0, 2

回答

mov      ebx, ebp
and      ebx, 1
sub      ebx, 1
pxor     xmm5, xmm5
pinsrd   xmm5, ebx, 2
blendvps xmm0, xmm4, xmm5

test    ebp, 1
jnz @F
blendps xmm0, xmm4, 4

xmm0 = _mm_cmplt_epi32( xmm1, _mm_setzero_si128() );
xmm0 = _mm_blendv_ps( xmm0, xmm3, xmm0 );

xmm0 = _mm_blendv_ps( _mm_setzero_si128(), xmm3, xmm1 );