为什么这个浮点运算编译得如此奇怪(没有优化)?

时间:2015-04-02 01:23:32

标签: c++ gcc optimization floating-point g++

我使用了0级和3级优化编译了以下函数 g++版本4.7.2 20120921

double function1(double a, double b)
{
  return (a+b)*(a+b);
}

反汇编0级优化版提供:

0000000000000000 <_Z9function1dd>:
   0:   55                      push   rbp
   1:   48 89 e5                mov    rbp,rsp
   4:   f2 0f 11 45 f8          movsd  QWORD PTR [rbp-0x8],xmm0
   9:   f2 0f 11 4d f0          movsd  QWORD PTR [rbp-0x10],xmm1
   e:   f2 0f 10 45 f8          movsd  xmm0,QWORD PTR [rbp-0x8]
  13:   66 0f 28 c8             movapd xmm1,xmm0
  17:   f2 0f 58 4d f0          addsd  xmm1,QWORD PTR [rbp-0x10]
  1c:   f2 0f 10 45 f8          movsd  xmm0,QWORD PTR [rbp-0x8]
  21:   f2 0f 58 45 f0          addsd  xmm0,QWORD PTR [rbp-0x10]
  26:   f2 0f 59 c1             mulsd  xmm0,xmm1
  2a:   f2 0f 11 45 e8          movsd  QWORD PTR [rbp-0x18],xmm0
  2f:   48 8b 45 e8             mov    rax,QWORD PTR [rbp-0x18]
  33:   48 89 45 e8             mov    QWORD PTR [rbp-0x18],rax
  37:   f2 0f 10 45 e8          movsd  xmm0,QWORD PTR [rbp-0x18]
  3c:   5d                      pop    rbp
  3d:   c3                      ret    

3级优化提供:

0000000000000000 <_Z9function1dd>:
   0:   f2 0f 58 c1             addsd  xmm0,xmm1
   4:   f2 0f 59 c0             mulsd  xmm0,xmm0
   8:   c3                      ret    

在未经优化的版本中,为什么代码会做那么多额外的工作?具体来说,是什么导致mulsd之后的4条指令?他们所做的就是将xmm0移至内存,从内存移至rax,然后返回内存,然后返回xmm0

1 个答案:

答案 0 :(得分:1)

这是编译器认为它正在做什么的一种可能的视图。 非优化编译器采用非常本地的视图。它没有展望下一步该做什么。它可能使用一组非常有限的操作。例如,它似乎将一些工作转移限制在堆栈和xmm0之间,或堆栈和rax之间。

“一刀切”方法的另一个方面是,如果某些事情需要做某些事情,它往往会一直完成。特别是,某些功能需要足够的寄存器来要求将参数和中间结果保存到堆栈中。优化编译器仅在必要时才执行此操作。非优化编译器无条件地执行它。

0000000000000000 <_Z9function1dd>:
// Push the stack
   0:   55                      push   rbp
   1:   48 89 e5                mov    rbp,rsp
// Save the parameters to stack temporaries
   4:   f2 0f 11 45 f8          movsd  QWORD PTR [rbp-0x8],xmm0
   9:   f2 0f 11 4d f0          movsd  QWORD PTR [rbp-0x10],xmm1
// Load the temporary representing a into register xmm1, via xmm0
   e:   f2 0f 10 45 f8          movsd  xmm0,QWORD PTR [rbp-0x8]
  13:   66 0f 28 c8             movapd xmm1,xmm0
// Add the temporary representing b leaving (a+b) in xmm1
  17:   f2 0f 58 4d f0          addsd  xmm1,QWORD PTR [rbp-0x10]
// Load the temporary representing a into xmm0
  1c:   f2 0f 10 45 f8          movsd  xmm0,QWORD PTR [rbp-0x8]
// Add the temporary representing b, leaving (a+b) in xmm0
  21:   f2 0f 58 45 f0          addsd  xmm0,QWORD PTR [rbp-0x10]
// Multiply (a+b)*(a+b)
  26:   f2 0f 59 c1             mulsd  xmm0,xmm1
// Store the multiply result in a stack temporary
  2a:   f2 0f 11 45 e8          movsd  QWORD PTR [rbp-0x18],xmm0
// Load the return value into rax
  2f:   48 8b 45 e8             mov    rax,QWORD PTR [rbp-0x18]
// Move the return value to xmm0 via a stack temporary
  33:   48 89 45 e8             mov    QWORD PTR [rbp-0x18],rax
  37:   f2 0f 10 45 e8          movsd  xmm0,QWORD PTR [rbp-0x18]
// and return
  3c:   5d                      pop    rbp
  3d:   c3                      ret