@无名啊，顺便一提，在安卓上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; }
@无名啊，如果特别担心，就用联合吧。
#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; }代码实际上变简单了。
不过有趣的是，使用联合的版本给出的都是负值（虽然也是正确解，平方根有两个解），不知道符号位的处理和非联合版本有什么不同。
@无名啊，我还是要说，错误行为不是未定义行为。
解引用指向float值的long指针具有明确的定义，因为float的内存表示在IEEE754定义，long的内存表示在C中定义。在特定的实现中，两者的长度可能相同，也可能不同，但当两者长度不同时，错误一定会以规定好的方式发生：float及其后不属于它的4字节会被访问。这只是编程错误，不是未定义行为。
@无名啊，这是这个函数的PHP版本，有助于理解为什么没有未定义行为：
<?php function Q_rsqrt(float $number) { $threehalfs = 1.5; $x2 = $number * 0.5; $y = $number; $i = unpack("l", pack("f", $y))[1]; $i = 0x5f3759df - ($i >> 1); $y = unpack("f", pack("l", $i))[1]; $y = $y * ( $threehalfs - ($x2 * $y * $y) ); $y = $y * ( $threehalfs - ($x2 * $y * $y) ); return $y; } 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));在给定的定义域和有效数字范围内，它和C版本的结果一致。如果继续增加输出的位数，结果就开始不一致了，因为PHP在内部使用64位整数和双精度浮点数，而非C代码的32位整数和单精度浮点数，只在
pack和unpack时才转换为32位单精度，所以两者会有精度差异。此外32位和64位在处理符号位上可能也有差异，所以C版给出负数解的情况下PHP给出的是正数解。当然两者都是正确的解，因为负数的平方也是正数。
C版本：
#include <stdio.h> float Q_rsqrt( float number ) { long i; float x2, y; const float threehalfs = 1.5F; x2 = number * 0.5F; y = number; i = * ( long * ) &y; // evil floating point bit level hacking i = 0x5f3759df - ( i >> 1 ); // what the fuck? y = * ( float * ) &i; 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; }
老虎会游泳
