MIPS32下IEEE 754双精度数相加时的结果长度判定问题

嘿，我太懂你这种卡在浮点归一化小数点上的头疼了——毕竟IEEE 754双精度的加法流程里，最容易绕晕的就是同指数相加后的进位处理，尤其是两个1.xxxxx格式的尾数加出1x.xxxxx的时候。其实你不用死盯着小数点位置，换个思路用整数运算的逻辑来处理，问题就清晰多了：

• 直接抓进位标志，比盯小数点靠谱
  双精度的尾数是隐含最高位1的52位小数，所以同指数相加的两个尾数本质上是1.f（f是52位二进制小数）的形式。当你把它们相加时，要么结果在[1.0, 2.0)之间，要么刚好等于2.0（极端情况）。
  在MIPS里做加法时，你可以直接关注加法的进位位：

  1. 把两个尾数（连同隐含的最高位1）当成整数来加（比如1.f等价于整数(1 << 52) + f，直接用dadd指令）
  2. 如果加法产生了最高位的进位（也就是结果的第53位是1），就说明结果是1x.xxxxx的形式
  3. 这时候把结果右移1位（相当于除以2，变回1.xxxxx的归一化格式），同时把指数加1（因为右移一位等价于指数+1，保证数值不变）

• 你的长度判定思路可以落地，但要抓对“有效长度”
  你说的“判定结果长度是否超过输入数”是对的，但这里的“长度”要指包含隐含位的总有效位数。输入的每个尾数是53位（1个隐含位+52个小数位），如果相加后的结果是54位（多了进位的1），那就必须归一化。
  在代码里可以这么实现：

  • 把相加后的结果和1 << 53做比较（因为53位的最大值是(1 << 53) - 1）
  • 如果结果大于等于这个值，说明有进位，执行右移+指数加1的操作

• 别忘了极端边界情况
  比如两个全1的尾数相加（1.111...111 + 1.111...111），结果刚好是2.0，右移后变成1.000...000，指数加1，这完全符合IEEE 754规范；另外如果相加结果是1.000...000（比如一个尾数加0），那就直接跳过归一化步骤就行。

给你一段MIPS伪代码参考，一看就明白：

# 假设a0和a1存着同指数的尾数（已包含隐含最高位1），exp存着当前指数
dadd a2, a0, a1       # 执行尾数相加
srl t0, a2, 53        # 提取第53位（进位标志位）
beqz t0, skip_norm    # 无进位则跳过归一化
srl a2, a2, 1         # 尾数右移一位，回归1.xxxxx格式
addi exp, exp, 1      # 指数加1，补偿右移的倍数
skip_norm:
# 现在a2是归一化后的尾数，exp是调整后的指数

这样处理下来，根本不用费力跟踪小数点在哪，用整数运算的进位和移位逻辑就能完美搞定归一化～

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