设为首页收藏本站

爱吱声

 找回密码
 注册
搜索
查看: 5506|回复: 25
打印 上一主题 下一主题

[信息技术] C++ 提速的新发现

[复制链接]
  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    跳转到指定楼层
    楼主
     楼主| 发表于 2022-9-24 22:54:26 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
    C++ 比 Octave 慢好多,怎么破?
    7 b0 t; B0 ^2 y5 a
    1 K$ G+ M, t; Y2 _: n自相关两层循环,内层循环涉及浮点数计算,试验了一下把内层循环内部全都 comment out 只留个壳子,  但空的内层循环本身就把速度拉下来了,看来问题并不在浮点计算。2 _: m7 S4 K* n; o( o

    8 R) D8 u, C3 H9 K" `速度优化问题真的很有意思啊。& l2 F) F, d8 n' `2 W) T

    ) _9 {  m" ?) ^8 {' o- e欢迎大家继续讨论

    评分

    参与人数 2爱元 +8 收起 理由
    helloworld + 4
    住在乡下 + 4 涨姿势

    查看全部评分

  • TA的每日心情
    开心
    10 小时前
  • 签到天数: 1925 天

    [LV.Master]无

    沙发
    发表于 2022-9-24 23:04:16 | 只看该作者
    拉下来?拉多少?; d8 N* m/ e! U& s. [
    把代码贴上来看看?5 V% ~, h8 N& p$ S3 G3 D/ T
    7 [/ K/ r5 {& s" ?, ~) T' ]
    难道分支预测不准破坏流水线执行?不该啊。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    板凳
    发表于 2022-9-24 23:15:41 | 只看该作者
    会不会代码本身的缺陷阻止了自动优化?另外,硬件配置和开发环境可能也有关系。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    2025-6-14 10:13
  • 签到天数: 124 天

    [LV.7]分神

    地板
    发表于 2022-9-24 23:33:02 | 只看该作者
    Maybe Debug mode?
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    5#
     楼主| 发表于 2022-9-24 23:54:10 | 只看该作者
    本帖最后由 雷达 于 2022-9-24 23:57 编辑 + ]5 _, g& Q7 T5 H
    数值分析 发表于 2022-9-24 23:043 i' [5 L& \. N3 c) m
    拉下来?拉多少?
    * @- X4 y. E" z/ p" N# ]& }+ R把代码贴上来看看?

    9 g7 f6 P3 H: z. P  N
    . D, b  l( E6 w1 M% Bvoid xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB), r5 y" {8 ~; x: Z8 P/ F0 U- }
    {3 }8 C5 c; t2 C( ^! D5 w9 T
            comp temp, xtimesy;5 b$ B- r- c& z9 v# p
            xtimesy.re = 0;7 I& |  T/ t& f
            xtimesy.im = 0;
    8 h% K2 X! P5 R% U  S+ P7 R        int j0 = lenB - 1;% J, Y' \0 I& Y+ }# Q: w
            int    i, j, i1, reali;
    ' `: h; u% a3 k$ y" Q        if (lenA % 2 == 1)5 o* c2 v; {& O" \5 q
                    reali = lenA + 1;6 l7 ]0 s% G3 e6 i% F: F
            else
    ; {- E, [/ w5 x) t, }                reali = lenA;
    : T1 J/ J; n2 K" T$ V5 e        reali /= 2;
    ! X1 [% S, H2 ^2 y. T6 P
    4 H# B2 n6 t4 ]* I- L        int nconv = reali + lenB;8 T1 V6 n: O% _1 j
            //#pragma omp parallel for
    6 o- X6 z0 i) {7 M4 W        for (i = reali; i < nconv; i++)$ B3 N1 E+ L4 z8 L- J2 j5 m* M7 y
            {
    # K" c+ o+ J: X- t' q                temp.re = 0;$ Z, N9 I2 p4 W0 Z
                    temp.im = 0;# \3 ^+ ]( O: e( d
                    i1 = i;8 t/ }- o5 W7 F! @2 k
                    for (j = j0; j >= 0; j--)
    / F! }* Y) ~5 N2 A                {
    . p4 f/ S  b# V; \1 K: t/ [                        /* floating date operation */# R% c) P* n$ n6 `! j- U; E7 |
                    }
    ; q/ N" X$ d- c8 }
            }
    & O* x2 s. G- m7 e, R3 m}
    ! D' Z) ?% }. V& u# f. r
    0 R  @* ^/ `2 [, B- g5 ^xcorr函数代码如上,comp是复数struct, 做过长度为11、19两个矢量的测试,和octave结果完全一样
    4 {1 g3 I' E6 d( ?  a+ w
    3 ^' C) ]- ~0 P3 X* G# j5 I红色部分是内循环,现在其内部操作都comment out 了, j0大概是 6000。
    4 M! J4 O  W2 {8 |+ ^现在call xcorr 100次,耗时78s.+ U- ~# E. T+ A# b3 v' h: o6 p0 z# h
    3 P) i: |5 B8 w0 ?' V' ]
    如果把红色部分内循环本身完全comment out, call xcorr 1000次,耗时 <1s. / C/ j) e9 P1 ~% I, Z

      C; [; E2 C# `" M
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    6#
     楼主| 发表于 2022-9-25 00:17:17 | 只看该作者
    风雨无阻 发表于 2022-9-24 23:338 S# I/ Q( j  s+ E3 g
    Maybe Debug mode?

    % u# f1 \+ m: V) D9 }; z) ]
    7 n  E3 U4 V( o/ D3 x$ H6 `4 C不应该,看我上面的回复。9 H3 x! \% I. t& J2 x3 b' [

    % n7 o! M2 p3 B- b: e* b* F, O# G我更怀疑是 VS 社区版的问题
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    10 小时前
  • 签到天数: 1925 天

    [LV.Master]无

    7#
    发表于 2022-9-25 00:20:10 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 00:24 编辑
    3 f; W1 Y& i9 f% f3 a$ E$ h
    雷达 发表于 2022-9-24 23:54$ E4 e* P! \- h( y8 y. P
    void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
    ; E  B) z. k5 X{
    8 a' \: J" k0 t2 ~        comp temp, xtimesy;
    8 l# x/ g0 K% A
    5 t3 d& E9 |# t
    这个不是这么比的吧。。。
    # u2 A! D7 i- s! p# B) R+ O/ \) A/ |7 T" X" N7 ^% Z7 s
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    3 L" i% s2 v) j3 X6 O; h  i1 q9 R8 C$ R! l: @7 n' m
    而加上内循环,光jmp和dec指令就至少多执行了6000个,慢个几十倍不是正常的么?
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    8#
     楼主| 发表于 2022-9-25 00:46:56 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 01:09 编辑 9 I2 D$ a' g3 t* `: y) L9 W
    数值分析 发表于 2022-9-25 00:20
    9 a, g5 h) A* ]+ I这个不是这么比的吧。。。
    $ I1 g6 S/ B" `( X% A0 v2 B
    2 z6 Q( D/ t  `* L+ e您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    + R* ?' d  K. K8 _* \2 M. c
    * X5 Q( u: Q1 r: V5 P; @$ b) b7 p
    有道理。! d" P$ ]6 r- E8 T
    所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,记得 numeric 库里有算向量内积的,我回头试试。
    ; v: ?  e# K1 v+ U" T
    & y8 |, J2 k& }6 K  G' U# {我先尝试尽量用标准库,一个小程序,不想搞得太复杂。多谢了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    9#
    发表于 2022-9-25 01:27:53 | 只看该作者
    雷达 发表于 2022-9-25 00:469 Y; M: M% j8 o/ `+ l  [
    有道理。2 Q9 n% U! k' P; v6 J# [& Z
    所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...

    % ?* e3 u, C) d, |- @你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子,把循环本身翻译成机器指令loop或dec/jnz,两者速度上会差很多% O# Z/ U0 Q# I4 w% P
    Why is the loop instruction slow? Couldn't Intel have implemented it efficiently?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    10#
    发表于 2022-9-25 01:48:51 | 只看该作者
    数值分析 发表于 2022-9-25 00:20
    6 |; x2 Z8 A/ {& k这个不是这么比的吧。。。
    " H# s( k, W( s) i# h+ J9 k4 ~' k; S. s* _
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    而加上内循环,光jmp和dec指令就至少多执行了6000个

    9 K6 J7 U/ [& ]; S5 {/ V- l$ \+ W8 E" s3 n( I
    现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果循环这样跑,花不了多少时间。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    10 小时前
  • 签到天数: 1925 天

    [LV.Master]无

    11#
    发表于 2022-9-25 02:06:44 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 02:16 编辑 3 E( s- A2 ^5 E# _7 T) F
    沉宝 发表于 2022-9-25 01:48! `  R2 ^( \+ J, S* R- e+ `
    现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果 ...
    # t8 P) ~7 n  q" X, |1 ?
    " j+ N2 D3 H# Q
    是的,兄台说的对。" t# A5 j- G) @0 j& K, {0 Y+ O

    , s1 Z3 [& y8 T; D; B$ c其实我想说的是 真正数值计算部分和代码中其他不直接计算的overhead的比值这个事儿。
    $ M* L  T7 q" z) g( y4 \9 r- o% B1 u6 k& {$ [# C) B, l
    雷达兄构造测试用例的时候,屏蔽掉了所有计算的部分,使得剩下的都是overhead,这样run time比较的结果就显得好像不合理了。如果把计算加回去,计算部分的run time会dominate,结果就不那么离谱了。因为不好说,所以用指令数对比的方式试图直观地说明这一点。
    7 m! W2 u  x+ Q  E5 G2 ]- I
    2 w7 F/ A1 X7 r' [比如说,如果有计算,那么跑六千个循环相对于计算应该用不了多少时间。但是如果一边是什么都不做,另一边是六千个循环,那六千个循环比什么都不做慢几十倍了,就不是那么不合理了。
    - @- d% b5 v- W, W+ _. p% x. U
    % z6 u1 P% n1 H- x' B+ _* \0 V当然也有可能像兄台说的,是优化参数的问题,但我觉得更多地是测试用例设计的不合理。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    12#
     楼主| 发表于 2022-9-25 04:47:43 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 04:49 编辑 ; ~7 ^, ?. {6 @( ~$ u( |
    沉宝 发表于 2022-9-25 01:272 @. w7 a/ U" Z/ y
    你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子 ...
    , h5 x/ z; [% Q) N

    % h) o# Z. |: M3 K) }又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差一倍,我上面这个差的太多了。
    4 Q) J; l; Y3 w6 R, L2 H
    2 K! J3 l. W4 [$ @7 D; s& O5 l6 c+ q我已经完全懵了。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    13#
    发表于 2022-9-25 05:51:23 | 只看该作者
    雷达 发表于 2022-9-25 04:47
    & s3 W7 V" |4 Z4 [0 y又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...

    & C( H: s3 c# _6 d7 z时间差一倍的结果可以接受。3 e  ^. t. Y. j! b8 Z

    3 z/ T( J4 D/ s# g# o你还是用profile工具看看吧。现在大家都主观瞎猜。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    10 小时前
  • 签到天数: 1925 天

    [LV.Master]无

    14#
    发表于 2022-9-25 14:58:59 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 15:38 编辑
    . F' u! X2 P; S9 w, f
    雷达 发表于 2022-9-25 04:47
    0 L* [; o* m5 X' e; V又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...

    & D' B( X( X# O- l  |( o$ d
    0 x- @) j1 ?) D  H# |% _; q1 x; p: k: r- Q
    ; K0 ?0 k8 I) V( n8 x% A
    能不能把这个也贴上来,看看和上一个有什么不同?
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    15#
     楼主| 发表于 2022-9-26 01:30:03 | 只看该作者
    本帖最后由 雷达 于 2022-9-27 01:17 编辑 / d4 W( X( ~& @% P3 R! d
    数值分析 发表于 2022-9-25 14:58+ S' V( Z& s: A& ?7 a, f: d
    能不能把这个也贴上来,看看和上一个有什么不同?

    7 o* [' C" N* r3 N1 P理了理思路,重新做了一个测试。5 A$ o% ~, N9 \2 U6 h7 S: n5 ?8 c1 y
    做了两个 vector 和 两个 float *, 都长 100000; g8 @6 {3 V% E9 a
    外循环 6000,里面先做随机数生成,模拟真实环境,避免数据的 cache.) ?' r& J! X; |

    ) r* `0 v; F2 ~( i6 @内循环试了4种方法,
    - Y# E, |- h- j7 N6 g1. 直接调用 vector inner_product 247s ) J' m- T  ^1 n9 |" l6 A
    2. vector 循环点乘累加 237s6 [2 ?6 q5 H5 l2 M7 c7 {
    3. float * 循环点乘累加 204s& g: ]6 t8 U( n1 h' E' C
    4. 空循环 100000 次 202s
    - l* [& e& M2 U. `: S! C
    % w+ k! a6 J, H0 o( R7 R! r不做内循环 200s
    + e1 o! q( [' J! q) s$ d% ~0 A0 e2 }; ]1 w
    你昨天说的对,内循环本身占比是很小的,大头在其他处理。( _2 `4 o! [% z. v; a
    另外可以看到, float * 循环点乘累加 并不差,比用vector 还更快。! O) a( C9 T! c1 E) G9 R
    # t" u- R! u/ N; D& A3 ]/ c; A
    至于我那个原始程序,还有一些疑问,见5楼,其他都不变仅仅是有无空的内循环就有很大不同,这是不对的,也许有一些其他缺陷我没有看到。(也许可以改成 while 试试)0 U* G  w* [- e# I* ]7 ]/ C2 e
    5 O( U0 e& a4 |; h
    (为什么下面我贴的  b1 加 方括号里的 i , 显示出来却是 b1 ?方括号 i 消失了。 LOL . 改成  jj 好了,原来 方括号里的 i 是斜体标志  LOL)
    ; X+ l- W8 b) s0 g: F: L! V
    4 w4 d% e9 E2 @7 K
            std::vector < float > vec1(N);! v' T  ?# ~5 G3 w, ?* r  c& `
            std::vector < float > vec2(N);) n4 N& L5 G4 }  {3 h1 B/ L8 n
            float* b1 = new float[N];
    $ L8 j& {. o" Y$ t9 O( \        float* b2 = new float[N];1 z- }# L# z. [9 {. N5 L
    * x. z2 O2 f8 r1 n' y4 l
            for (int j = 0; j < 6000; j++)2 i& o- A' }% n4 b8 p2 u
            {
    $ k' I6 p& c; v1 S8 i) ?7 X1 C+ h. n                std::generate(vec1.begin(), vec1.end(), []() {
    1 v8 y( u) O% H5 q                        return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
    7 K, n/ a* V% L( c  ]6 H+ m                        });/ _9 o- b. v3 L3 P8 M

    ! u2 [: l2 H$ n- P                std::generate(vec2.begin(), vec2.end(), []() {' Q2 _0 O9 u8 [
                            return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
    / g8 B5 k* [0 [7 m! m! o! a: j/ Z                        });
    & L* g4 ?" |3 x, C2 G. X: ?/ X0 E
    6 T, t: F# C0 L+ y, z4 D6 W$ s                for (size_t jj = 0; jj < vec1.size(); jj++)# @; {7 q+ s: S( T: @9 y8 ?" l/ n
                    {: I' B: |7 ]* s8 X. }1 A, ?6 O  g
                            b1[jj] = vec1[jj];7 m, c- q, E3 ^0 o& F( B9 L5 M3 g
                    }
    $ ^1 B6 @4 v. p0 S2 r% L5 S6 i$ B9 g& @5 O' q$ P, A& ~
                    for (size_t jj = 0; jj < vec2.size(); jj++)
    ) H8 a% B8 t  l! p$ l; A/ ]; d                {
      l7 R1 {8 o1 Q/ t, F' `                        b2[jj] = vec2[jj];
    ' r4 ~: j4 C6 T! m! n- H                }$ [* H0 Z# ]0 B' c+ H1 Q* q

    : N0 F* {. u5 m! W                //Method - 1  N=100000 247s  . k, Y3 A" c# }) \* n
                    //fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);
    3 T$ a3 I# A' Y$ ~) a- X                                / n' `, k$ E& E( {
                    //Method - 2  N=100000  237s% n6 \  y% C0 B! _& Q* {- w5 Q
                    /*& V% |1 x" C% g5 i1 r; h
                    for (int jj = 0; jj < N ; jj++)4 @5 A1 a- k* o/ V
                    {
    8 _, Y# o, m; H% K                        fresult += vec1[jj] * vec2[jj];
    : _& q! s# e/ y" c$ X: b* q                }
    7 X/ M* U2 w2 t5 \- u% G                */
      N/ S- M1 N1 z+ m/ R                                # v, ]4 Q+ l1 w) A2 E2 s
                    //Method - 3  N=100000 204s4 X, M7 W6 a/ D  A
                    /*! L; D0 D; g4 b' }( \+ w
                    for (int jj = 0; jj < N; jj++); A% k9 H$ ~+ W- `; S, g- ?5 @$ C* ]
                    {4 e8 C2 I2 E; J: \
                            fresult += b1[jj] * b2[jj];$ I7 N9 E7 }# L) R
                    }9 l6 t/ u8 v' X0 T
                    */
    2 O1 o1 g: U+ c4 l
    * {, t, C8 u5 i2 X6 }& w0 i                //Method - 4   202s
    9 w9 b0 n# y1 [! x# S                /*
    % F9 Z( q! A3 f  v% }# X                for (int jj = 0; jj < N; jj++)' T9 {. P3 |" M$ K8 }" R$ D$ T
                    {) V& R! ^% H- d6 F1 s. g
                            
    1 [- x& Q3 p  a) E( }& ]9 J5 s" L                }
    3 M7 a6 B; n9 q6 }/ F                */# ~0 G$ j7 d4 y0 l) G1 v
                    //comment out all methods, N=100000  202s                2 C# p' e7 r5 _+ Z5 j. Z
            }4 x& D$ P# B9 w9 N

    9 t* G9 f1 p1 U$ ^( I  D4 T" @        delete []b1;3 J7 o. u0 F3 S5 s1 B: j
            delete []b2;
    & X; _1 B& |7 \* M
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    奋斗
    2020-2-8 10:08
  • 签到天数: 2 天

    [LV.1]炼气

    16#
    发表于 2022-9-27 00:15:38 | 只看该作者
    瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?% z; a3 b! f" s
    8 J2 a; p+ I, D1 K# l- u+ I
    你第二个试验里面的j在循环里面又重新定义了啊,你确定真的跑了6000次?! r9 c# O5 t% F/ N, h) R5 c
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    擦汗
    2024-12-25 23:22
  • 签到天数: 1182 天

    [LV.10]大乘

    17#
     楼主| 发表于 2022-9-27 01:16:03 | 只看该作者
    机器猫 发表于 2022-9-27 00:15- D% W* l* U% L0 [/ G" S; x
    瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?
    ' A) ^6 r' N% W% C5 y4 k0 s0 L
    " B: `8 D( L; }8 _. ~你第二个试验里面的j在循环里面又重新定义 ...

    6 `6 d, J! _) V6 G- S内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL8 w" Z' x3 ^& w6 g( s3 m# Z
    ' ?$ a' S# S8 }5 U" l8 K
    不和它较劲了,瞎耽误工夫,我已经转到 ubuntu, 也准备顺便试试 avx2 向量化。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    奋斗
    2020-2-8 10:08
  • 签到天数: 2 天

    [LV.1]炼气

    18#
    发表于 2022-9-27 02:06:58 | 只看该作者
    雷达 发表于 2022-9-27 01:16. N+ S; a5 i; o% f2 K$ i. o- ]! U6 I# [
    内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
    / p* H- C$ |- @( @+ z5 V, n3 \$ z
    + Y1 c: _6 g! S不和它 ...
    : J/ b" b# {- `; W6 {5 m8 x/ X
    ; I2 O7 w4 q" U4 |3 u. z( n( L
    不过可以试试我说的register变量。前一个试验j是混在一堆其它变量里一起定义的,很有可能是在stack上,这样内存读写会更多,要是再碰上每次都需要加载cache就更慢了。
    2 l7 ?, ?$ j( y1 h6 ^+ x后面一个是在循环那里定义的,说不定编译器就把它优化成register变量了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    19#
    发表于 2022-9-27 07:25:57 | 只看该作者
    一个无关问题,为什么爱坛的帖子里在我这里有好些奇怪的东东在里面,是防拷贝措施吗?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    20#
    发表于 2022-9-27 20:29:50 | 只看该作者
    雷达 发表于 2022-9-24 23:543 N2 c1 B# {7 D% a+ y# O% L8 d
    void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
    ! v1 M0 X- d. a. r0 I7 f{: r3 F. Z% W  N( M
            comp temp, xtimesy;
    $ V6 O' g) C2 G7 Q% [* h1 c# \
    这个code里面如果Openmp没有被注释掉的话,那么temp那个变量应该是定义在循环里面,否则线程之间会存在争夺写入那个temp的风险。
    ' c# ~5 Q9 M" ?  ~内层for循环如果没有内部操作的话,编译时应该被优化掉了,和你完全注册掉整个循环是一回事。可能你的编译设置没有打开优化?$ I+ m6 z" i( a3 H
    VS社区版没有问题,我工作用的就是社区版,设置正常的话不会比商业版差。以前游说头头用Intel Compiler,他说不想花钱,而且差不了多少,就一直用到现在。
    回复 支持 反对

    使用道具 举报

    手机版|小黑屋|Archiver|网站错误报告|爱吱声   

    GMT+8, 2025-7-3 12:57 , Processed in 0.048621 second(s), 19 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

    快速回复 返回顶部 返回列表