这段时间在学习DSP浮点运算转定点运算过程中,在网上看到了一位网友给出的一段程序.如下所示:
---------------------------------------------------------
声明:
此文为原创,欢迎转载,转载请保留如下信息
作者:聂飞(afreez) 北京-中关村
联系方式:afreez@sina.com (欢迎与作者交流)
初次发布时间:2006-11-28
不经本人同意,不得用语商业或赢利性质目的,否则,作者有权追究相关责任!
---------------------------------------------------------
/***************************************************************
// 将32位浮点数fval转换为32位整数并存储在ival中
// 小数部分将被裁剪掉
//
void TruncToInt32 (int &ival, float fval)
{
ival = *(int *)&fval;
// 提取尾数
// 注意实际的尾数前面还有一个被省略掉的1
int mantissa = (ival & 0x07fffff) | 0x800000;
// 提取指数
// 以23分界,指数大于23则左移,否则右移
// 由于指数用偏移表示,所以23+127=150
int exponent = 150 - ((ival >> 23) & 0xff);
if (exponent < 0)
ival = (mantissa << -exponent);
else
ival = (mantissa >> exponent);
// 如果小于0,则将结果取反
if ((*(int *)&fval) & 0x80000000)
ival = -ival;
}
*******************************************************/
程序流程写的很清晰.兴奋之余马上在自己的机器上测试了一把,结果发现总是不正确,大体如下:
//主函数
...
int ival;
float fval=4.23;
TruncToInt32 (&ival, fval);
//这里的输出结果不正确
printf("ival=%d/n",ival);
...
//函数定义
void TruncToInt32 (int &ival, float fval)
{
//skip
...
// 如果小于0,则将结果取反
if ((*(int *)&fval) & 0x80000000)
ival = -ival;
//发现在这里输出的结果是正确的
printf("ival=%d/n",ival);
}
仔细看了一眼,才发现,原来是函数调用时参数传递有问题。主要上理解指针和指针的指针,大家对照我的思路考虑考虑吧。
修改主程序,如下所示:
...
int ival;
float fval=4.23;
RightTruncToInt32 (&ival, fval);//第n行,假定是经过修正后的函数
ival=(int)fval;//第n+1行
printf("ival=%d/n",ival);
...
对n行和n+1行进行时间效率测试,显示,第n+1行所需的计算时间少于第n行。
对主程序进一步改善,大致如下:
...
int ival;
float fval=4.23;
//-------Begin more fast section?-----------
// 提取尾数
int ival= ((*(int *)(&fval)) & 0x07fffff) | 0x800000;
// 提取指数
int exponent = 150 - (((*(int *)(&fval)) >> 23) & 0xff);
if (exponent < 0)
ival = (ival<< -exponent);
else
ival = (ival >> exponent);
// 如果小于0,则将结果取反
if ((*(int *)&fval) & 0x80000000)
ival = -ival;
//-------End more fast section?-----------
ival=(int)fval;//第n+1行
...
结果显示,在RedHat 9.0平台上测试显示,如果编译不加优化选项,more fast section的执行时间效率优与第n+1行;加优化选项,第n+1行的执行时间效率优与more fast section。
分析结论:
由于more fast section还有一些代码行可以优化,我想,在一些特定场合下,还是可以借鉴这个思路的。