C++ Div3、Sqrt 函数高性能实现（带汇编指令集）

liulilittle

于 2023-12-08 11:53:23 发布

阅读量1.2k

点赞数 8

分类专栏： C/C++ 文章标签： c++ 汇编服务器

本文链接：https://blog.csdn.net/liulilittle/article/details/134874910

版权

C/C++ 专栏收录该内容

262 篇文章 7 订阅

订阅专栏

文章探讨了在C++中使用特定的魔法数字（如0x5f375a86和0x5f3759df）实现浮点数平方根计算（floatSqrt）的方法，以及针对不同处理器架构（Intel和AT&T）优化的整数除以3（Div3）的代码示例。

摘要由CSDN通过智能技术生成

均采用魔法数字（Magic Number）实现，一个是经典求平方根函数所使用的魔法数字：0x5f375a86、0x5f3759df。

float Sqrt(float x) noexcept { /* 0x5f3759df */
    float xhalf = 0.5f * x;
    int32_t i = *(int32_t*)&x;
    i = 0x5f375a86 - (i >> 1);
    x = *(float*)&i;
    x = x * (1.5f - xhalf * x * x);
    return x;
}

unsigned int Div3(unsigned int i) noexcept {
    // AT&T:
    // movl    $2863311531, %edx
    // imulq   %rcx, %rdx
    // shrq    $33, %rdx

    // INTEL:
    // mov     edx, 2863311531
    // imul    rdx, rcx
    // shr     rdx, 33

    unsigned long long n = static_cast<long long>(i) * 2863311531;
    unsigned int r = static_cast<unsigned int>(n >> 33);
    return r;
}

unsigned long long Div3(unsigned long long i) noexcept {
    // INTEL:
    // mov     rax, -6148914691236517205; 
    // mul     rcx
    // shr     rdx, 1

    int64_t rax = -6148914691236517205;
    int64_t rdx = (rax * i) >> 1;
    return rdx;
}