Intel x86_64 汇编，如何在 x87 和 SSE2 之间移动？（计算双精度数的反正切）

Question

我尝试计算双精度浮点值的反正切，该值保存在 xmm 寄存器中。对于正常的浮点，可以使用旧的 x87 指令 FPATAN，但是如何使用 double 来做到这一点？

Answer 1

您仍然可以将数据从 xmm 复制到 x87 以使用类似 的指令fpatan，但通常您应该调用数学库函数。（fpatan速度太慢，用许多简单的指令替换它仍然很好。） Wikipedia 建议查看 Netlib 以获取可自由重新分发的 C 实现。（显然，最简单的方法是在您使用的任何系统上调用 libm 中的函数。）

如果您打算这样做，请不要使用静态存储来存储您跳转的内存；在堆栈上使用临时的。

另请注意，fpatan需要 2 个输入，因为它实现了atan2库函数，根据两个输入的符号在适当的象限中给出结果。

; assuming you did  sub  rsp, 24   or something earlier in your function

movsd   [rsp], xmm1
fld     qword [rsp]   ; st0 = xmm1
movsd   [rsp], xmm0
fld     qword [rsp]   ; st0 = xmm0,  st1 = xmm1

fpatan                ; st0 = arctan(xmm1/xmm0)

fstp    qword [rsp]     ; x87 stack is empty again
movsd   xmm0, [rsp]   ; xmm0 = arctan(xmm1/xmm0)

; and then   add rsp, 24   at some point before returning