设为首页收藏本站

爱吱声

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

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

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

    [LV.10]大乘

    跳转到指定楼层
    楼主
     楼主| 发表于 2022-9-24 22:54:26 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
    C++ 比 Octave 慢好多,怎么破?+ ?6 V. S- l! E. f7 l8 r3 ?$ |& ~

    # g0 F4 B) o) l/ A2 R5 b自相关两层循环,内层循环涉及浮点数计算,试验了一下把内层循环内部全都 comment out 只留个壳子,  但空的内层循环本身就把速度拉下来了,看来问题并不在浮点计算。+ O- k1 e& N( w, C7 J
      K, M, f2 r0 z$ o. L
    速度优化问题真的很有意思啊。1 B' Z/ a+ e! V+ z. O2 U0 ^& T$ B

    & l3 f/ h9 v! y0 _9 p2 y欢迎大家继续讨论

    评分

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

    查看全部评分

  • TA的每日心情
    开心
    前天 04:43
  • 签到天数: 1942 天

    [LV.Master]无

    沙发
    发表于 2022-9-24 23:04:16 | 只看该作者
    拉下来?拉多少?/ U7 K  P5 F4 f3 o- @- H
    把代码贴上来看看?5 S4 G1 s- Z5 u, B4 s7 ]

    + l0 ?# K6 I  ]3 e( c" J难道分支预测不准破坏流水线执行?不该啊。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    板凳
    发表于 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 编辑 : ]2 r, g  X# m8 P9 W) @
    数值分析 发表于 2022-9-24 23:04
    3 u' O% h- Z+ Q/ G  C( f+ U拉下来?拉多少?
    - f! m! c, b, `9 T7 e# N/ p! g把代码贴上来看看?
    / C$ N, d& g! o; c4 m; X% o! R6 T

    : y( j7 Q2 r4 e0 Bvoid xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)0 P5 H" [4 O" p$ k. F3 ?
    {
    , ?5 Y% s: j4 b* Z( K9 [        comp temp, xtimesy;, d, \$ @' y7 n: e! B4 @9 p
            xtimesy.re = 0;
    2 W  T! `! ]( _) \4 z" z        xtimesy.im = 0;" s2 ~& w0 o  D5 k* l! e0 U: }9 ~" T
            int j0 = lenB - 1;
    2 J9 p7 w6 H" e+ a2 W        int    i, j, i1, reali;$ A$ g/ N( O, S5 _+ [6 z3 c* u% R
            if (lenA % 2 == 1)
    8 \/ O- p# \) s                reali = lenA + 1;0 I% [2 g( [! U- Y  N: L% {0 u
            else
    9 X. s6 W" c- Y                reali = lenA;
    " m6 {6 L5 n) J3 Z( Y; B' W- H3 `        reali /= 2;  ^: e: v, Y) t1 \
    3 T0 B& J) {1 m6 Y- {
            int nconv = reali + lenB;
    ! [- h3 o* `! p" S) V        //#pragma omp parallel for
    & P+ [1 T# `! q% _' A$ X) J. C' l: \6 f' L        for (i = reali; i < nconv; i++)
    0 O% V$ l6 U( P( j        {
    2 O  N  r/ Q; u                temp.re = 0;
    ) X. ?# A, V: O1 A, i$ Y& O; ~                temp.im = 0;
    ' k- k3 I8 k. Q) ~, \                i1 = i;
    , h* p% P  v( ^. h* c6 [) i" Y                for (j = j0; j >= 0; j--)
    & k8 X# Q9 D; x- D2 t: h                {( M/ [# d3 D2 S' P5 l
                            /* floating date operation */' R* J) _; w# H$ ~: w
                    }

    : q/ D/ a- X- k        }
    3 W) x5 N6 q% e: H$ [# w+ }$ q}+ R0 C5 ~! [" o, E6 I: i  W
    * ]( @: Q: y9 X8 X6 N) f' ], R7 {; E
    xcorr函数代码如上,comp是复数struct, 做过长度为11、19两个矢量的测试,和octave结果完全一样
    ) S1 `/ z; g0 W  a  c6 M7 W* R( D) {: W, J( y
    红色部分是内循环,现在其内部操作都comment out 了, j0大概是 6000。
    ! b  P* I/ h  W现在call xcorr 100次,耗时78s.
    & Y4 O) _; z/ X/ t
    3 J2 k" O7 ?  ~( ?- g; r+ n9 }如果把红色部分内循环本身完全comment out, call xcorr 1000次,耗时 <1s. # m9 O: H/ T7 N) B, ]# l. _+ Q

    ' W3 q  [4 J7 I# p8 r7 ?
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    6#
     楼主| 发表于 2022-9-25 00:17:17 | 只看该作者
    风雨无阻 发表于 2022-9-24 23:334 Z- N+ Y- E0 n- |. q" v  g1 A
    Maybe Debug mode?
    ( H( l& o  [3 J3 E( f7 P

    ' L3 {2 a) {6 l) P' \6 o  A: M; K6 D不应该,看我上面的回复。
    ; Q3 n$ Z8 E% ^4 E3 e" x4 C8 e; S* ]6 R5 z* K. e1 x  Z1 H
    我更怀疑是 VS 社区版的问题
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    前天 04:43
  • 签到天数: 1942 天

    [LV.Master]无

    7#
    发表于 2022-9-25 00:20:10 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 00:24 编辑 5 a, u  Q7 ~9 g6 H; r- x
    雷达 发表于 2022-9-24 23:54
    . N8 h# ^! p% Z8 r- q8 r9 H0 cvoid xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
    ' Q, h8 c) P+ p7 {( B- `. \4 u{+ j- P  _6 \% q5 p( }, x8 C
            comp temp, xtimesy;
    0 o% H( l" j" g7 E

    2 z# @: f. g" @3 ~/ o这个不是这么比的吧。。。
    + |9 W- _2 b3 n: g/ I7 J, v4 {; D1 f5 Q: a, h/ V! q
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。3 d8 y3 `# S/ K

    3 ]# m: T' d0 W( u& ]) B. I% [而加上内循环,光jmp和dec指令就至少多执行了6000个,慢个几十倍不是正常的么?
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    8#
     楼主| 发表于 2022-9-25 00:46:56 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 01:09 编辑
    6 p% c0 J* t' Z1 A8 ?  f" z: E
    数值分析 发表于 2022-9-25 00:20& Y+ a" L' y6 j5 d+ |9 J
    这个不是这么比的吧。。。4 e& S0 }( D$ z" \

    ) W5 Y- ^/ J* C  g您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    " R8 B3 b3 A- {
    1 [2 H5 [' M2 m
    有道理。
    2 s& _* g' h3 G9 o+ R- G# }; e所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,记得 numeric 库里有算向量内积的,我回头试试。
    2 W2 d0 g0 ^4 Z6 z: \5 M6 Y1 e
    2 D% r' ]" S0 e' ~0 t( t, A我先尝试尽量用标准库,一个小程序,不想搞得太复杂。多谢了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    9#
    发表于 2022-9-25 01:27:53 | 只看该作者
    雷达 发表于 2022-9-25 00:46, j! V0 C& e0 k+ G' I' \/ A
    有道理。
    % B' U' B! y' o: z所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...

    7 J3 a. ]8 j/ {3 f* t你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子,把循环本身翻译成机器指令loop或dec/jnz,两者速度上会差很多
    2 U+ k4 V" ]1 t% S/ aWhy is the loop instruction slow? Couldn't Intel have implemented it efficiently?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    10#
    发表于 2022-9-25 01:48:51 | 只看该作者
    数值分析 发表于 2022-9-25 00:20
    0 _; {  V# \* T- R8 o1 E这个不是这么比的吧。。。  u4 J: P* ?. X/ B& q
    ! X' ^$ N0 y' j6 m$ \2 D
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    而加上内循环,光jmp和dec指令就至少多执行了6000个
    $ f8 h+ N! i& @% W9 ?0 |& A

    2 y$ ?. R" p  g  ]- y现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果循环这样跑,花不了多少时间。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    前天 04:43
  • 签到天数: 1942 天

    [LV.Master]无

    11#
    发表于 2022-9-25 02:06:44 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 02:16 编辑 ' b/ i# \( b, u/ y
    沉宝 发表于 2022-9-25 01:48, U# ]3 s/ X9 O) a4 X
    现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果 ...

    % h1 V+ M9 g6 `, @7 ~9 y2 B! Q& I* t, M3 I1 p
    是的,兄台说的对。6 v( u/ ^; |4 u

    2 `0 T) c. j3 R/ z3 @$ |" y其实我想说的是 真正数值计算部分和代码中其他不直接计算的overhead的比值这个事儿。
    & c1 [+ r9 V( k0 R% m, ^! d$ \- z' Y- i( n  l' n8 u
    雷达兄构造测试用例的时候,屏蔽掉了所有计算的部分,使得剩下的都是overhead,这样run time比较的结果就显得好像不合理了。如果把计算加回去,计算部分的run time会dominate,结果就不那么离谱了。因为不好说,所以用指令数对比的方式试图直观地说明这一点。
    3 t; R% T/ ~6 m& r  r3 q
    1 N4 H5 U6 n+ e7 B比如说,如果有计算,那么跑六千个循环相对于计算应该用不了多少时间。但是如果一边是什么都不做,另一边是六千个循环,那六千个循环比什么都不做慢几十倍了,就不是那么不合理了。, e6 m* x1 ^2 D: ]; q$ [
    ( T! m2 W; P) n; m" v
    当然也有可能像兄台说的,是优化参数的问题,但我觉得更多地是测试用例设计的不合理。
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    12#
     楼主| 发表于 2022-9-25 04:47:43 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 04:49 编辑
    1 C% @# ?& S$ x7 E9 [% J0 A
    沉宝 发表于 2022-9-25 01:27
      m9 N" ]: W9 n0 m! m( T% W你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子 ...
    ; J/ y5 k' p) O. N/ `6 l3 ?

    / e+ i3 D4 v! y9 s4 }6 L+ l又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差一倍,我上面这个差的太多了。
      z8 \" O9 D0 u( Q" i! c
    0 I2 _, _2 J9 Z9 L我已经完全懵了。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    13#
    发表于 2022-9-25 05:51:23 | 只看该作者
    雷达 发表于 2022-9-25 04:47
      }& c0 u! ^' A( k% j- |又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...

    : o* W, u: r2 O8 z0 P! k  t时间差一倍的结果可以接受。
    ! N6 i8 l5 ~$ m* h; q6 u0 n# D
    0 {/ E8 f9 ^. ?; f5 C+ o, L你还是用profile工具看看吧。现在大家都主观瞎猜。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情
    开心
    前天 04:43
  • 签到天数: 1942 天

    [LV.Master]无

    14#
    发表于 2022-9-25 14:58:59 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 15:38 编辑
    ; M$ O: |, e1 P) M
    雷达 发表于 2022-9-25 04:47
    & i& g; P, M# ?2 n又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...

    / u% J( D4 g) D  q6 {
    " ^6 f: ^9 F* j9 c. k. k: M; {. ^
    & d2 f4 T( ]; g9 w) I9 C  y. ?2 D7 @% n' b2 p/ e- ]
    能不能把这个也贴上来,看看和上一个有什么不同?
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    15#
     楼主| 发表于 2022-9-26 01:30:03 | 只看该作者
    本帖最后由 雷达 于 2022-9-27 01:17 编辑 + G, M, _8 `6 \" l: \1 M8 [( S
    数值分析 发表于 2022-9-25 14:58
    7 T  \1 _# @( i) T' L$ c* W能不能把这个也贴上来,看看和上一个有什么不同?
    6 d$ y# U$ t( F* }4 N
    理了理思路,重新做了一个测试。
    1 D" |2 \$ X7 }7 m3 s2 A9 b做了两个 vector 和 两个 float *, 都长 100000
    $ F3 Q; `$ p. n, v+ e外循环 6000,里面先做随机数生成,模拟真实环境,避免数据的 cache.9 M+ s# y2 F8 u1 S% B

    ' u" H* ]/ Y. d/ f3 \- ~( a内循环试了4种方法,
    # y- P# q& o% T" |1. 直接调用 vector inner_product 247s
    . q2 x& s" O  C7 f; K2. vector 循环点乘累加 237s
    9 H) G( D: N, l$ S4 }. l3. float * 循环点乘累加 204s
    & s8 c% e. V# `* I9 h* G4. 空循环 100000 次 202s7 ~5 r, H$ `  u2 {2 [% q
    7 ?5 Q5 k% W& C4 t. Z1 [/ o
    不做内循环 200s+ H( X5 C2 g  f  E, ?0 w

    2 C& W, ^# V9 T. E4 @9 a8 S你昨天说的对,内循环本身占比是很小的,大头在其他处理。
    ' V6 P$ a( x) b! o$ J另外可以看到, float * 循环点乘累加 并不差,比用vector 还更快。
    % ?- d4 \$ \; @! @8 s4 \& a1 i- w4 L" z7 x4 U/ q; I
    至于我那个原始程序,还有一些疑问,见5楼,其他都不变仅仅是有无空的内循环就有很大不同,这是不对的,也许有一些其他缺陷我没有看到。(也许可以改成 while 试试)8 T; l1 O8 {) W6 T2 L
    ! O( S* V3 f% [5 T  V2 S1 a
    (为什么下面我贴的  b1 加 方括号里的 i , 显示出来却是 b1 ?方括号 i 消失了。 LOL . 改成  jj 好了,原来 方括号里的 i 是斜体标志  LOL)( ]) y: i: m4 Q: R7 h: u
    : J' ~3 T+ b. A; X2 F4 e' B& l
            std::vector < float > vec1(N);
    7 q. r, v$ @7 K) I" t        std::vector < float > vec2(N);6 G& w/ |: R5 L3 `- u9 W
            float* b1 = new float[N];
    ) U6 P8 J/ l! q3 B4 ^' |. j+ p8 @' k' X        float* b2 = new float[N];
    : v- Z; W+ b$ _& J( m  d8 o/ y7 Y! Q- m# L
            for (int j = 0; j < 6000; j++)
    , i+ r1 u4 ?2 @" h4 B        {3 M$ ^; ?" f# D( Y$ i
                    std::generate(vec1.begin(), vec1.end(), []() {# J5 `" t$ L: S, b: R# u' @
                            return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;4 }2 [! Y) g" M3 h) `8 d
                            });
    - R; z2 {3 l/ k% v' K2 v. n- k$ y& x4 h2 p1 G$ i6 z
                    std::generate(vec2.begin(), vec2.end(), []() {7 l# i5 y0 J" E1 @" q
                            return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;6 U# c  b0 k8 {; Q
                            });
    * z  c4 l' B) l9 x) J, b: l9 S/ A2 |4 X* e' C
                    for (size_t jj = 0; jj < vec1.size(); jj++)
    2 u6 A0 x! q' j, q5 @% X+ f; G" i                {' A  [: v: l# d. U9 J) n( n
                            b1[jj] = vec1[jj];
    : h' N3 J( j5 L/ x% l4 b                }
    ! ^! c1 B: C, H6 l, y8 Q- B0 @
      O7 A9 r( ?1 R# z- M                for (size_t jj = 0; jj < vec2.size(); jj++)
    1 a. t$ M" F& t                {
    0 ?5 h4 o0 P8 ^# Q                        b2[jj] = vec2[jj];
    6 t4 P9 Q3 ?$ K6 u: }4 k                }
    : x0 W7 k' U. F6 n( q: l8 O" Y9 _% H  Y/ [
                    //Method - 1  N=100000 247s  2 V7 G- `( L( G7 F$ O
                    //fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);% x) G( g; R5 V. j7 X8 K
                                    1 ^9 r  K. q2 D  c
                    //Method - 2  N=100000  237s* N: p! P/ ]2 i  x( w& S7 k& k+ c
                    /*) F0 E; g9 S1 o( x2 B# E: h/ {2 w
                    for (int jj = 0; jj < N ; jj++)7 r$ T& L. F6 K+ k& K/ F+ z4 V# U% S9 S
                    {
    8 p9 a7 I, O8 i                        fresult += vec1[jj] * vec2[jj];
    & ^) G8 Y  s1 V5 y" C                }
    $ }/ v6 n& ~+ F                */
    % j1 L2 T$ a% Y# j& D$ H# j1 N  e( t                                
    : w; r( j  g& f* C5 b6 P                //Method - 3  N=100000 204s5 ^- E' |2 k6 i1 C& D
                    /*8 w/ c4 f- X* {; e
                    for (int jj = 0; jj < N; jj++)2 S( H- `, a. ?" K
                    {) p8 c7 Z" A% S. o1 V: Y
                            fresult += b1[jj] * b2[jj];: F+ k3 O& {( O
                    }
    2 m( {- \2 ]8 ?2 \# j                */; W8 V$ M) Q& |% P/ o9 `

    " C0 j8 K, g& o. ]/ @; I, z                //Method - 4   202s$ N  ~) o( x% ]+ N5 ]! S0 r% [( Y
                    /*) G# X7 y) _4 p# v
                    for (int jj = 0; jj < N; jj++)
    & l" z$ z8 c+ a' F. |  ~                {0 Q* Q1 P9 Z/ J) ~9 m
                            + w; o* ^5 ~; \
                    }
    7 [, C& w5 a; j6 V  n8 n                */
    4 u" j5 p/ {0 w% ]                //comment out all methods, N=100000  202s                , U3 G6 A2 H. [' j/ I/ r3 D9 M  k' E
            }* R5 V8 h; M2 D2 n1 c

    ! ]" R# n' J' [+ v& f        delete []b1;
    ; R0 B/ p% C' V$ y9 O        delete []b2;

    3 K$ N! v7 f, ^5 ^, S2 K* R: [
    回复 支持 反对

    使用道具 举报

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

    [LV.1]炼气

    16#
    发表于 2022-9-27 00:15:38 | 只看该作者
    瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?: C3 T9 f% |, R! n9 A) ?6 H  p% C
    8 ~9 e- A# y4 H& X
    你第二个试验里面的j在循环里面又重新定义了啊,你确定真的跑了6000次?9 x! r$ y  E! s" o8 ]
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    17#
     楼主| 发表于 2022-9-27 01:16:03 | 只看该作者
    机器猫 发表于 2022-9-27 00:15
    ; t; G# P& o5 X+ _  T5 z6 x瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?
    % V5 s: n3 G) J
    2 Z  d( i! Z3 l$ g你第二个试验里面的j在循环里面又重新定义 ...
    / @* s1 A0 R9 U/ I7 H
    内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL( L% P( L- i. ~5 V- ]6 l2 \. m

    9 B! L3 }8 i! f% M6 g不和它较劲了,瞎耽误工夫,我已经转到 ubuntu, 也准备顺便试试 avx2 向量化。
    回复 支持 反对

    使用道具 举报

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

    [LV.1]炼气

    18#
    发表于 2022-9-27 02:06:58 | 只看该作者
    雷达 发表于 2022-9-27 01:16, }; [6 S2 {% V- X
    内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL) A  b) F4 e* c/ ]  K' [# r6 G
    3 `/ `9 `6 p7 H: o5 w% B
    不和它 ...

    . D8 {; N! m! n) S  C
    / ^$ u! n0 R, a& o: s7 X& D不过可以试试我说的register变量。前一个试验j是混在一堆其它变量里一起定义的,很有可能是在stack上,这样内存读写会更多,要是再碰上每次都需要加载cache就更慢了。
    : _1 y" W+ l7 A- G后面一个是在循环那里定义的,说不定编译器就把它优化成register变量了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

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

    使用道具 举报

    该用户从未签到

    20#
    发表于 2022-9-27 20:29:50 | 只看该作者
    雷达 发表于 2022-9-24 23:548 Y) O; f4 q* T. n( S
    void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
    8 M- D# u% b. [" P{
    2 x8 F) P8 T* ?' c' K" j$ G        comp temp, xtimesy;
    6 y, R: g  W- O! ?! j
    这个code里面如果Openmp没有被注释掉的话,那么temp那个变量应该是定义在循环里面,否则线程之间会存在争夺写入那个temp的风险。
    " l/ k* O' ~2 W- z7 L内层for循环如果没有内部操作的话,编译时应该被优化掉了,和你完全注册掉整个循环是一回事。可能你的编译设置没有打开优化?
      u4 s( _8 X1 p& hVS社区版没有问题,我工作用的就是社区版,设置正常的话不会比商业版差。以前游说头头用Intel Compiler,他说不想花钱,而且差不了多少,就一直用到现在。
    回复 支持 反对

    使用道具 举报

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

    GMT+8, 2025-9-20 20:53 , Processed in 0.046204 second(s), 19 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

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