首页 > 解决方案 > do_cos_slow.isra 是什么?

问题描述

我写了一个简单的代码来测试prof.

double bar_compute (double d) {
    double t = std::abs(d);
    t += std::sqrt(d);
    t += std::cos(d);

    return t;
}

// Do some computation n times
double foo_compute(unsigned n) {
    std::random_device rd;
    std::mt19937 mt(rd());
    std::uniform_real_distribution<double> dist(0.0, 1.0);

    double total = 0;
    for (int i=0; i<n; i++) {
        double d = dist(mt);
        total += bar_compute(d);
    }
    return total;
}

当我运行prof并查看输出时

  56.14%  runcode  libm-2.23.so       [.] __cos_avx                                                                                                                                
  27.34%  runcode  runcode            [.] _Z11foo_computej                                                                                                                         
  13.92%  runcode  runcode            [.] _Z11bar_computed                                                                                                                         
   0.86%  runcode  libm-2.23.so       [.] do_cos_slow.isra.1                                                                                                                       
   0.44%  runcode  runcode            [.] cos@plt                                                                                                                                  
   0.41%  runcode  libm-2.23.so       [.] sloww1                                                                                                                                   
   0.35%  runcode  libm-2.23.so       [.] __dubcos                                                                                                                                 
   0.17%  runcode  ld-2.23.so         [.] _dl_lookup_symbol_x             

是什么do_cos_slow.isra意思sloww1

cos我可以使用更快的版本吗?不然怎么会叫慢呢?

标签: glibcavxperftrigonometry

解决方案


do_cos_slow来自它在glibc/sysdeps/ieee754/dbl-64/s_sin.c中的声明。之所以调用它,是因为它比第164 行声明上方的注释do_cos_slow所基于的函数更精确。do_cos

.isra是因为该函数是 IPA SRA 根据以下堆栈溢出答案优化的版本,GCC 函数后缀“isra”是什么意思?

sloww1是一个根据上面的注释计算 sin(x+dx) 的函数。

关于更快的cos版本,我不确定是否有更快的版本,但是如果您将提供libm的glibc或libc实现更新到至少glibc 2.28,那么您将得到Wilco Dijkstra删除这些慢路径的结果dosincos 的功能和重构,可提高速度。

提交消息

Refactor the sincos implementation - rather than rely on odd partial inlining
of preprocessed portions from sin and cos, explicitly write out the cases.
This makes sincos much easier to maintain and provides an additional 16-20%
speedup between 0 and 2^27.  The overall speedup of sincos is 48% over this range.
Between 0 and PI it is 66% faster.

您可以尝试的其他替代方案是其他 libc 或 libm 实现,或其他 cos 实现,包括avx_mathfun或avx_mathfun 以及针对较新的 GCC或supersimd的一些修复


推荐阅读