@无名啊，顺便一提，在安卓上long是64位，这就是符号位差异的来源
还有，
memcpy看起来是最佳选择，因为它没有任何多余的操作——我们想要的就是内存复制，所以我们就应该写内存复制。改成memcpy后，代码比用联合与指针类型转换都简单。#include <stdio.h> #include <stdint.h> float Q_rsqrt( float number ) { int32_t i = 0; const float threehalfs = 1.5F; float x2 = number * 0.5F; float y = number; memcpy(&i, &y, sizeof(y)); // evil floating point bit level hacking i = 0x5f3759df - ( i >> 1 ); // what the fuck? memcpy(&y, &i, sizeof(y)); y = y * ( threehalfs - ( x2 * y * y ) ); // 1st iteration y = y * ( threehalfs - ( x2 * y * y ) ); // 2nd iteration, this can be removed return y; } int main() { printf("sizeof(long): %lu\n", sizeof(long)); printf("sizeof(int32_t): %lu\n", sizeof(int32_t)); printf("sizeof(float): %lu\n", sizeof(float)); printf("%0.7f\n", Q_rsqrt(3.14)); printf("%0.7f\n", Q_rsqrt(1024.0)); printf("%0.7f\n", Q_rsqrt(10086.0)); printf("%0.7f\n", Q_rsqrt(2147483647.0)); printf("%0.14f\n", Q_rsqrt(3.14)); printf("%0.14f\n", Q_rsqrt(1024.0)); printf("%0.14f\n", Q_rsqrt(10086.0)); printf("%0.14f\n", Q_rsqrt(2147483647.0)); return 0; }
@老虎会游泳，另外，我这篇帖子，主要是想讨论
restrict的（因为我已经认下严格别名了）。。我刚才看了看你一开始的其他回复。
我认为，根据 cppreference 对
restrict的叙述，应该是允许有其他别名的。只要保证：通过
restrict修饰的指针写入某个对象后，不再使用其他别名访问这个对象，即可？
@老虎会游泳，关于联合体，应该只有 C 能用，C++ 不能：
联合体的大小仅足以保有其最大的数据成员。其他数据成员在该最大成员的一部分相同的字节分配。分配的细节是实现定义的，且读取并非最近写入的联合体成员是未定义行为。许多编译器以非标准语言扩展实现读取联合体的不活跃成员的能力。
@老虎会游泳，比如，
(long *)&y，编译器认为结果是NULL。。。
@无名啊，如果特别担心，就用联合吧。
#include <stdio.h> float Q_rsqrt( float number ) { union { long l; float f; } i; float x2, y; const float threehalfs = 1.5F; x2 = number * 0.5F; y = number; i.f = y; // evil floating point bit level hacking i.l = 0x5f3759df - ( i.l >> 1 ); // what the fuck? y = i.f; y = y * ( threehalfs - ( x2 * y * y ) ); // 1st iteration y = y * ( threehalfs - ( x2 * y * y ) ); // 2nd iteration, this can be removed return y; } int main() { printf("%0.7f\n", Q_rsqrt(3.14)); printf("%0.7f\n", Q_rsqrt(1024.0)); printf("%0.7f\n", Q_rsqrt(10086.0)); printf("%0.7f\n", Q_rsqrt(2147483647.0)); printf("%0.14f\n", Q_rsqrt(3.14)); printf("%0.14f\n", Q_rsqrt(1024.0)); printf("%0.14f\n", Q_rsqrt(10086.0)); printf("%0.14f\n", Q_rsqrt(2147483647.0)); return 0; }代码实际上变简单了。
不过有趣的是，使用联合的版本给出的都是负值（虽然也是正确解，平方根有两个解），不知道符号位的处理和非联合版本有什么不同。
@老虎会游泳，
( long * ) &y这个是左值吗？解引用后，类型是long，不是字符类型，与float也不是兼容类型，所以应该是未定义行为。我 C/C++ 写得不多，现在看到那篇文章，就更怕编译器没能正确实现 UB 了。。
@老虎会游泳，cppreference 说，这是 UB：
@老虎会游泳，按照标准，可能根本不会发生解引用指向float值的long指针，因为可能已经被优化掉了。。也就没有你后面说的浮点定义如何如何……
@老虎会游泳，如果按 cppreference 所说，应该是未定义行为。
转换指针类型没问题，只要不访问就行。（但不访问，转了也没用。可认为不能转）
比如：
float a = 1.0; long * b = (long *)&a; *b = 1; return a;按照标准，编译器可认为，
a未被修改（因为*b不是a的兼容类型，所以修改*b不应该污染a），所以优化掉b，直接返回1.0。。。
@无名啊，我还是要说，错误行为不是未定义行为。
解引用指向float值的long指针具有明确的定义，因为float的内存表示在IEEE754定义，long的内存表示在C中定义。在特定的实现中，两者的长度可能相同，也可能不同，但当两者长度不同时，错误一定会以规定好的方式发生：float及其后不属于它的4字节会被访问。这只是编程错误，不是未定义行为。
老虎会游泳
