内联汇编代码并存储 128 位结果

Question

我必须在 Ubuntu 上使用 GCC 编译器在我的 C 程序中编写以下语句作为内联汇编代码。

__int128 X = (__int128)F[0]*T[0]+(__int128)F[1]*T[1]+(__int128)F[2]*T[2]);

其中F为无符号64位整数数组，T为有符号64位整数数组。数组 F 作为参数通过引用传递，而 T 是本地数组。以上声明我翻译如下

__asm__("movq %0, %%rax;        imulq  %1;  movq %%rax, %%xmm0; movq %%rdx, %%xmm1;"
        ::"m"(F[0]), "m"(T[0]));
__asm__("movq %0, %%rax;        imulq  %1;  movq %%xmm0, %%rcx; addq %%rcx, %%rax;" 
        "movq %%rax, %%xmm0;                movq %%xmm1, %%rcx; adcq %%rcx, %%rdx;"
        "movq %%rdx, %%xmm1;"
        ::"m"(F[1]), "m"(T[1]));
__asm__("movq %2, %%rax;        imulq  %3;  movq %%xmm0, %%rcx; addq %%rcx, %%rax;" 
        "movq %%rax, %?;                    movq %%xmm1, %%rcx; adcq %%rcx, %%rdx;"
        "movq %%rdx, %?;"
        :"=m"(??), "=m"(??):"m"(F[2]), "m"(T[2]));

第一个也是最重要的问题是我做得对吗？如果是，那么我不知道如何将结果存储到 X 中，因为结果的低 64 位在 rax 中，高 64 位在 rdx 中。我已经检查过我是否可以替换？？通过 X，然后我得到了错误的结果。

使用xmm寄存器来存储是有原因的。因此，由于我对内联汇编还很天真，所以我认为有更好的方法可以做到这一点。我用上面的内联汇编代码检查了我的程序，没有错误。任何帮助或改进建议将不胜感激。

Answer 1

您正在对 F 进行符号扩展。由于没有带符号 * 无符号的乘法指令，因此必须显式地进行符号扩展（16 位到 32 位的示例）：

(0xFFFF0000 + S) * U
= 0xFFFF0000 * U + S * U
= (0x100000000 - 0x10000) * U + S * U
= 0x100000000 * U - 0x10000 * U + S * U
= S * U - 0x10000 * U  (don't care about high bits)

您不能依赖内联 asm 语句块之间寄存器中剩余的值；你必须使用变量。所有修改的寄存器必须声明为输出或破坏。

例如，一个 64 位无符号值 U 和一个 64 位有符号值 S 的单个乘法：

__int128 X;
uint64_t Utmp = U;
asm ("mov %1, %%rax;"
     "mul %2;"
     "test %2, %2;"
     "cmovns %3, %1;"
     "sub %1, %%rdx"
     : "=&A" (X), "+r" (Utmp) : "r" (S), "rm" (0UL));

编辑：无需零输入即可完成：

int64_t Stmp = S;
asm ("mov %1, %%rax;"
     "mul %2;"
     "sar , %1;"
     "and %2, %1;"
     "sub %1, %%rdx"
     : "=&A" (X), "+rm" (Stmp) : "r" (U));

内联汇编代码并存储 128 位结果

Inline assembly code and storing 128-bit result

c

gcc

inline-assembly