x86汇编实现表达式计算报错-1.#IND00，求排查与验证工具

问题背景 #

需要计算数学表达式：(k - y) / (b^z - (k + k^a))，其中b、z、y、a为实数，k取值为1到15。以下是编写的x86汇编实现代码：

finit
    mov ecx, 15
    mov eax, 1
    
    ; Calcul b^z
    fld b 
    fld z   ;st(0)=z,st(1)=b,st(2)=k/y
    fyl2x; st(0)=b*log*z
    fld1
    fscale; st(0)=2^[b*logz];st(1)=b*log*z
    fld1 
    fld st(1)
    fyl2x ; st(0)=[b*logz],st(1)=2^[b^logz],st(2)=b*logz
    fsubr st(0),st(2);st(0)=blogz-[b^logz]
    f2xm1
    fld1
    fadd
    fstp qword ptr [putere]
    finit
  
loop_k:
    mov [k], eax

    ; Conversie k la float
    fild dword ptr [k]
    fstp qword ptr [k_float]

    ; Calcul (k - y)
    fld qword ptr [k_float]   ; Load k
    fld qword ptr [y]         ; Load y
    fsub                      ; k - y
    fstp qword ptr [numarator]

    ; Calcul k^a
    fld k_float 
    fld a   ;st(0)=z,st(1)=b,st(2)=k/y
    fyl2x; st(0)=b*log*z
    fld1
    fscale; st(0)=2^[b*logz];st(1)=b*log*z
    fld1 
    fld st(1)
    fyl2x ; st(0)=[b*logz],st(1)=2^[b^logz],st(2)=b*logz
    fsubr st(0),st(2);st(0)=blogz-[b^logz]
    f2xm1
    fld1
    fadd
    

    ; Calcul k + k^a
    fld qword ptr [k_float]   ; Load k
    fadd                      ; k + k^a
    fstp qword ptr [sumk]

  
    ; Calcul b^z - (k + k^a)
    fld qword ptr [putere]
    fld qword ptr [sumk]
    fsub                      ; b^z - (k + k^a)

    ; Calcul (k - y) / (b^z - (k + k^a))
    fld qword ptr [numarator]
    fdivr                     ; (k - y) / (b^z - (k + k^a))
    
    fld qword ptr [sum]
    fadd
    fstp qword ptr [sum]


    ; Incrementare k și buclă
    inc eax
    loop loop_k  
  
    ; Afișarea rezultatului
    push dword ptr [sum+4]
    push dword ptr [sum]
    push offset result_msg
    call printf
    add esp, 12

现存问题 #

• 不确定代码中的浮点数缩放逻辑是否正确，希望推荐可验证该逻辑的工具
• 无论输入什么数值，程序输出结果均为-1.#IND00，仅单循环时能输出值但结果不正确

问题解答 #

一、浮点数缩放逻辑验证工具 #

• GDB调试：用GDB单步执行汇编代码，通过info float查看FPU寄存器实时值，x/q查看内存中浮点数变量内容，逐指令对比手动推导的数学结果，验证缩放、幂运算的每一步正确性。
• C语言对照实现：编写等价的C语言代码实现相同数学表达式（包括手动实现幂运算底层逻辑），编译后对比汇编代码的中间结果与C程序输出，快速定位逻辑差异。
• 手动数学推导：针对每一步浮点数指令对应数学公式验证：
  1. fyl2x功能为计算ST(1) * log2(ST(0))，执行后弹出ST(0)，结果存于ST(0)
  2. fscale将ST(0)乘以2^ST(1)，仅处理指数的整数部分
  3. f2xm1计算2^ST(0) - 1，要求输入范围在[-1,1]之间
     把指令序列对应成数学公式，检查是否符合x^y = 2^(y*log2(x))的推导逻辑。

二、-1.#IND00错误排查与修复 #

-1.#IND00表示无效浮点数运算（如除以0、负数开方、对数输入负数等），结合代码分析问题点：

1. 分母计算符号错误
   代码中计算b^z - (k + k^a)时，fsub实际执行ST(0) = ST(0) - ST(1)，即计算sumk - putere，与需求的putere - sumk完全相反。应改为：

   fld qword ptr [sumk]
   fld qword ptr [putere]
   fsub                      ; putere - sumk
   

   符号错误可能导致分母为0或无效值，触发异常。

2. 幂运算逻辑错误
   代码中实现b^z和k^a的逻辑未正确分离指数的整数与小数部分，正确的非整数指数幂运算步骤应为：

   ; 计算 x^y，x、y为内存中的浮点数
   fld qword ptr [x]
   fld qword ptr [y]
   fyl2x               ; ST(0) = y * log2(x)
   fld st0             ; 复制指数值到ST(0)
   frndint             ; ST(0) = 指数的整数部分E
   fsub st1, st0       ; ST(1) = 指数的小数部分F = 原指数 - E
   f2xm1               ; ST(0) = 2^F - 1
   fld1
   fadd                ; ST(0) = 2^F
   fscale              ; ST(0) = 2^F * 2^E = 2^(E+F) = x^y
   fstp st1            ; 清理栈，结果保留在ST(0)
   

   原代码逻辑混乱导致幂运算结果错误，进而分母可能出现0或无效值。

3. FPU状态被意外重置
   计算完b^z后调用的finit会重置FPU所有寄存器和状态，属于多余操作，可能引发异常，建议删除。

4. sum变量未初始化
   代码未对sum变量赋值为0.0，第一次循环加载的是内存垃圾值，导致后续计算始终无效。需在循环前添加：

   fldz
   fstp qword ptr [sum]
   

5. 注释错误干扰逻辑判断
   计算k^a部分的注释仍沿用b^z的内容（如st(0)=z,st(1)=b），需同步修改为对应变量注释，避免调试混淆。

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