汇编浮点指令异常

使用汇编进行浮点数开发时需要注意FPU堆栈平衡，原因是FPU寄存器（ST0~ST7）是以栈的形式操作的。当你使用完FPU寄存器时需要进行释放，即FLD指令匹配FST指令。下面是我开发中遇到的问题，

 

#include <stdio.h> #include <cmath> // 求e的x次幂 double my_exp(double x) { double value; __asm{ fld x; fldl2e; fmul st, st(1); fst st(1); frndint; fxch st(1); fsub st, st(1); f2xm1; fld1; faddp st(1), st; fscale; fstp value; // ST0仍然被占用 } return value; } int _tmain(int argc, _TCHAR* argv[]) { int i=0,j=0,k=0; double vx = 0.1; for(int i = 0; i < 20; i++,vx+=0.1) printf("exp(%f)=%.8f,my_exp(%f)=%.8f/r/n", vx, exp(vx), vx, my_exp(vx)); return 0; } 

 double my_exp(double x)函数在每次执行时都会占用一个FPU寄存器，这样当循环到8次时已经没有任何FUP寄存器可用，这时对任何FPU寄存器的操作都会溢出，包括fld1简单赋值语句。然而这样的溢出错误应用程序并不捕获，给调试带来非常大的麻烦。

 

处理方法就是平衡FPU堆栈，用到哪个就释放掉或者进行清除操作。针对上述代码可以修改如下：

第一种方法：

... faddp st(1), st; fscale; fstp st(1); // <- 释放ST1寄存器 fstp value; ...

 

第二种方法：

... faddp st(1), st; fscale; fstp value; finit; // <- 初始化FPU寄存器，清除堆栈 ...

finit指令也可以放在函数入口。

 

相关参考：

关于汇编80X87FPU浮点运算

http://www.cppblog.com/saga/archive/2007/01/17/17721.html

 

 恶斗TCC(Tiny C Compiler)

http://blog.csdn.net/soloist/archive/2009/07/12/4342331.aspx