解读单精度浮点运算x86汇编的逻辑、参数、寄存器作用及对应C实现

这段是x86-64 System V ABI编译的单精度浮点数组点积计算函数，核心逻辑是遍历两个等长单精度浮点数组，计算对应位置元素的乘积累加和，返回最终的单精度浮点结果。

逐行纠正原有理解错误 #

• 第1行unknown5:：函数入口标识，理解正确。
• 第2行pxor %xmm0, %xmm0：对xmm0做按位异或操作将所有位清0，作为累加和的初始值。这里xmm0不是传入的浮点参数，是承载最终结果的局部累加器，同时也符合System V约定中单精度浮点返回值的寄存器位置要求，理解基本正确。
• 第3行movl $0, %eax：x86-64架构下，写32位寄存器eax会自动将对应64位寄存器rax的高32位清零，因此这步等价于将64位寄存器rax置为0，不是为了设置整数返回值——这里rax是循环计数器，该写法是编译器的指令长度优化，比直接写movq $0, %rax更节省代码空间，原理解误以为这是在设置32位返回寄存器，属于错误认知。
• 第4行jmp .L13：跳转到循环判断位置，是循环结构的优化写法（先判断循环条件，再执行循环体），跳转目标理解正确。
• 第6行movss (%rdx,%rax,4), %xmm1：AT&T语法的内存寻址格式为(基址, 索引, 比例)，计算的实际访问地址是rdx + rax * 4，不是原理解的rdx + rax +4。因为单精度浮点数占4字节，用4作为比例因子实现数组索引到字节偏移的转换，这步是将rdx指向的数组中、索引为rax位置的单精度浮点数加载到xmm1，原寻址计算逻辑错误。
• 第7行mulss (%rsi,%rax,4), %xmm1：同样基于SIB寻址规则，访问rsi + rax*4地址对应的单精度浮点数，与xmm1中的值相乘后存回xmm1，即完成两个数组同索引元素的乘法，原寻址计算逻辑错误。
• 第8行addss %xmm1, %xmm0：将乘法结果累加到xmm0的累加器中，理解正确。
• 第9行addq $1, %rax：循环计数器加1，对应数组索引自增，理解正确。
• 第11行cmpq %rdi, %rax：AT&T语法下的比较操作是计算rax - rdi的差值，仅修改标志位不存结果。配合下一行的jb .L14（无符号小于则跳转），逻辑是当rax（循环计数器）无符号小于rdi（数组长度）时，跳回循环体继续执行。原理解的跳转条件逻辑结果正确，但需要明确rdi的存储内容是数组长度。
• 第13行rep ret：rep前缀在这里是旧CPU平台的性能对齐优化，等价于普通的ret返回指令，无实际重复操作含义。

核心疑问解答 #

函数参数与寄存器含义 #

按照x86-64 System V调用约定，整数/指针参数依次存入rdi、rsi、rdx寄存器，结合代码使用逻辑，各位置的含义为：

• %rdi：第1个参数，size_t n，即两个浮点数组的元素个数，无符号64位整数，是循环的终止边界
• %rsi：第2个参数，const float *x，第一个单精度浮点数组的起始指针
• %rdx：第3个参数，const float *y，第二个单精度浮点数组的起始指针
• %rax：临时循环计数器，从0开始递增，存储当前访问的数组元素索引，不是返回值载体
• %xmm0：局部累加变量，存储点积的累加和，同时是单精度浮点返回值的存储位置，不是传入的浮点参数
• %xmm1：临时局部变量，存储每次循环中两个元素的乘积结果，每次循环都会被覆盖，不是传入的参数

返回值判定 #

• 最终返回值存储在%xmm0中，为计算得到的单精度浮点累加和
• %rax/%eax不是该函数的返回值载体，仅作为循环计数器临时使用。因为x86-64 System V约定中，rax是调用者保存的临时寄存器，被调用函数可随意修改，只要返回时将浮点结果放到xmm0、整数结果放到rax即可，本函数不返回整数，因此rax的内容不影响返回逻辑。

%rax的设计逻辑 #

这里rax被用作循环计数器，是编译器的常规选择：

• rax是调用者保存寄存器，函数内可随意修改不需要恢复
• 写32位的eax即可自动完成64位寄存器的清零，初始化循环计数器的指令长度更短
• 相比其他专用寄存器，用rax做计数不会破坏其他参数寄存器的内容

内容的提问来源于stack exchange，提问作者Pharros