在阅读Lua's源代码时,我注意到Lua使用macro将double舍入到32位int。我提取了macro,它看起来像这样:
union i_cast {double d; int i[2]};
#define double2int(i, d, t) \
{volatile union i_cast u; u.d = (d) + 6755399441055744.0; \
(i) = (t)u.i[ENDIANLOC];}
此处ENDIANLOC定义为endianness,0定义为小端,1定义为大端。 Lua小心翼翼地处理字节序。 t代表整数类型,例如int或unsigned int。
我做了一些研究,并且有一个更简单的macro格式使用了同样的想法:
#define double2int(i, d) \
{double t = ((d) + 6755399441055744.0); i = *((int *)(&t));}
或者用C ++风格:
inline int double2int(double d)
{
d += 6755399441055744.0;
return reinterpret_cast(d);
}
这个技巧可以在任何使用IEEE 754的机器上运行(这意味着今天几乎每台机器都有)。它适用于正数和负数,并且舍入遵循Banker's Rule。 (这并不令人惊讶,因为它遵循IEEE 754。)
我写了一个小程序来测试它:
int main()
{
double d = -12345678.9;
int i;
double2int(i, d)
printf("%d\n", i);
return 0;
}
按预期输出-12345679。
我想详细介绍这个棘手的macro是如何工作的。幻数6755399441055744.0实际上是2^51 + 2^52,或1.5 * 2^52,而二进制中的1.5可以表示为1.1。当任何32位整数被添加到这个神奇数字时,好吧,我从这里迷路了。这个技巧如何运作?
P.S:这是Lua源代码Llimits.h。
更新:
正如@Mysticial所指出的,这种方法并不局限于32位int,
只要数字在,它也可以扩展为64位int
范围为2 ^ 52。 (macro需要进行一些修改。)
有些材料称此方法无法在Direct3D中使用。
使用Microsoft汇编程序for x86时,有一个偶数
写在macro中的更快assembly(这也是从Lua源中提取的):
#define double2int(i,n) __asm {__asm fld n __asm fistp i}
单精度数有一个类似的幻数:1.5 * 2 ^23
醇>