ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

c – 如何实现没有未定义行为的快速逆sqrt?

2019-10-06 08:16:58  阅读:233  来源: 互联网

标签:strict-aliasing c undefined-behavior


参见英文答案 > What’s a proper way of type-punning a float to an int and vice-versa?                                    7个
根据我对strict aliasing rule的理解,fast inverse square root的此代码将导致C中的未定义行为:

float Q_rsqrt( float number )
{
    long i;
    float x2, y;
    const float threehalfs = 1.5F;

    x2 = number * 0.5F;
    y  = number;
    i  = * ( long * ) &y; // type punning
    i  = 0x5f3759df - ( i >> 1 );
    y  = * ( float * ) &i;
    y  = y * ( threehalfs - ( x2 * y * y ) );

    return y;
}

这段代码确实会导致UB吗?如果是,如何以符合标准的方式重新实现?如果没有,为什么不呢?

假设:在调用此函数之前,我们已经以某种方式检查浮点数是IEEE 754 32位格式,sizeof(long)== sizeof(float),平台是little-endian.

解决方法:

你应该使用memcpy. AFAIK这是唯一符合标准的方式,编译器足够聪明,可以用单个单词移动指令替换调用.有关这些说法背后的推理,请参阅this question.

标签:strict-aliasing,c,undefined-behavior
来源: https://codeday.me/bug/20191006/1859490.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有