浮点数比较技巧:内联汇编

时间:2014-03-09 07:38:24

标签: c++ gcc assembly x86 compare

很久以前,我使用这个简单的x86汇编技巧来获得浮点数比较的0或1:

fld [value1]
fcom [value2]
fnstsw ax
mov al, ah
and eax, 1

如果比较结果仅影响从一组2个值中选择值,则此技巧允许避免分支。在Pentium时代它很快,现在可能不会那么快,但谁知道呢。

现在我主要使用C ++并使用英特尔C ++编译器或GCC C ++编译器进行编译。

有人可以帮助将此代码重写为2个内置汇编程序风格(Intel和GCC)。

所需的函数原型是: inline int compareDoublesIndexed(const double value1,const double value2) { ... }

使用SSE2操作可能会更有效率。你的观点?


我试过这个:

__asm__(
    "fcomq %2, %0\n"
    "fnstsw %ax\n"
    "fsubq %2, %0\n"
    "andq $L80, %eax\n"
    "shrq $5, %eax\n"
    "fmulq (%3,%eax), %0\n"
    : "=f" (penv)
    : "0" (penv), "F" (env), "r" (c)
    : "eax" );

但我在英特尔C ++编译器中出错:浮点输出约束必须指定一个寄存器。

1 个答案:

答案 0 :(得分:2)

如你所述,自奔腾日以来情况发生了变化:

  • SSE现在是浮点而不是x87的首选指令集,即使对于标量运算
  • 也是如此
  • 优化编译器现在非常好

因此,首先检查编译器生成的内容,您可能会感到惊喜。我在下面的代码

上尝试了使用-O3的g ++

fcmp.cpp:

int compareDoublesIndexed( const double value1, const double value2 ) {
    return value1 < value2 ? 1 : 0;
}

这是编译器生成的

0000000000400690 <_Z21compareDoublesIndexeddd>:
  400690:       31 c0                   xor    %eax,%eax
  400692:       66 0f 2e c8             ucomisd %xmm0,%xmm1
  400696:       0f 97 c0                seta   %al
  400699:       c3                      retq   

这就是它的含义

  xor     %eax,%eax        ; EAX = 0
  ucomisd %xmm0,%xmm1      ; compare value2 (in %xmm1) with value1 (in %xmm0)
  seta    %al              ; AL = value2 > value1 ? 1 : 0

因此,编译器通过使用seta指令避免了条件分支(如果结果在上面,则设置字节为'1',否则为'0'。