@无名啊,你说得对,
( long * ) &y
是左值,* ( long * ) &y
是对它求值。不过
Q_rsqrt()
改写成符合标准似乎很容易。还有符号位的处理与union版本存在差异,所以确实可能涉及未定义行为。得到负数解可能才是这种内存操作应该有的结果。
@老虎会游泳,诶,好像只是访问时 UB 噢,转换指针类型应该没问题?
@老虎会游泳,比如,
(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 了。。
@无名啊,经过一番思考之后,我还是认为
i = * ( long * ) &y
没有问题,因为没有生成新的别名,整个表达式应该被视为一个右值。相反,把它分开的操作反而是有问题的,这违反了严格别名规则。
long i; long *p_i; float y; float *p_y; p_y = &y; p_i = (long *) p_y; i = *p_i;
但它应该也没有副作用,因为不涉及对
*p_i
的写入。
最重要的是,
i = * ( long * ) &y
的目标是读取y的值,它根本没有任何优化空间。&y
意味着y
一定得在内存,所以无论怎么优化,结果应该都是正确的。
@老虎会游泳,我没能力修改
gcc
,所以就认了。要么加-fno-strict-aliasing
,要么用memcpy
,来防止结果出错。但这还会影响到性能低下问题(特别是修改了
char *
后,编译器会认为一大堆对象有可能被修改了,所以缓存失效,需要重新读取),所以我要搞懂restrict
层主 @老虎会游泳 于 2023-01-27 16:16 删除了该楼层。
@老虎会游泳,cppreference 说,这是 UB:
@无名啊,这里没有未定义行为,因为取地址操作会阻止优化。因为
&a
,所以a
必须在内存,不能优化到寄存器。所以该代码没有未定义行为,但存在出现编程错误的风险(如果float和long长度不同)。float a = 1.0; long * b = (long *)&a; *b = 1; return a;
@老虎会游泳,按照标准,可能根本不会发生解引用指向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字节会被访问。这只是编程错误,不是未定义行为。
需要注意的是,错误行为不是未定义行为。
char c; long i; // 这个行为非常不恰当,会导致紧接着`c`后面的3个字节被访问,这3个字节不属于`c`。 // 但它只是错误行为,不是未定义行为。 // 这个行为会发生什么具有明确的定义,就是`c`所指向的内存地址及其后方3个字节一同被赋值给`i`,在所有平台上都会发生同样的事情。 // 所以,这里不含未定义行为,只含编程错误。 i = * ( long * ) &c;
@老虎会游泳,你看下 cppreference - 指针 - 注解 说的:
尽管任何指向对象的指针能被转型成指向其他类型对象的指针,解引用指向类型异于对象声明类型的指针几乎总是未定义行为。细节见严格别名使用。
@无名啊,我已经对上述问题进行了回答。
@无名啊,把这段代码拆分成多个部分,应该有助于理解为什么没有未定义行为:
long i; long *p_i; float y; float *p_y; p_y = &y; // 只是一个简单的取地址操作,不是未定义行为 p_i = (long *) p_y; // 对指针进行类型转换不是未定义行为,所有指针类型都是互相兼容的 i = *p_i; // i 和 *p_i 类型一致,没有未定义行为
操作的每一步都不含未定义行为,所以整体不含未定义行为。
@老虎会游泳,我比较怕的是未定义行为导致的结果错误,所以想弄清楚别名规则。
看到知乎那篇文章中的第三个例子,我又觉得
restrict
有助于减少性能损耗(修改char *
导致编译器认为this
可能被修改,进而没法重复利用缓存好的this
),所以顺便想弄清楚restrict