[英]In C89, how can I truncate and split a double precision floating point into two 32-bit words?
以C99为目标时,我可以使用以下代码截断并将double
值拆分为两个32位整数:
#include <stdint.h>
void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word)
{
uint64_t i = (uint64_t)d;
*most_significant_word = i >> 32;
*least_significant_word = i & 0xffffffff;
}
但是C89似乎没有定义64位整数类型,因此我不能使用编译器执行截断。 即使不需要截断(该值已经表示一个整数),我也无法使用&
或>>
类的位运算符&
因为这些运算符不适用于double
值。
那么,如何在纯C89中实现上面的split()
函数(从而不依赖于64位整数),返回组成存储在double
值中的53位整数的21/32位字?
忽略符号(与原始代码一样)应该像下面这样简单:
void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word)
{
*most_significant_word = d/4294967296.; // d>>32 in double
*least_significant_word = fmod(d, 4294967296.);
}
要考虑负数,请对绝对值进行运算,然后计算两者的补数:
void split(double d, unsigned long *most_significant_word, unsigned long *least_significant_word)
{
double dabs = d < 0 ? -d : d;
*most_significant_word = dabs/4294967296.;
*least_significant_word = fmod(dabs, 4294967296.);
if (d < 0) {
*most_significant_word = ~*most_significant_word;
*least_significant_word = ~*least_significant_word + 1;
if (!*least_significant_word) *most_significant_word += 1;
}
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.