设为首页收藏本站

爱吱声

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

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

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

    [LV.10]大乘

    跳转到指定楼层
    楼主
     楼主| 发表于 2022-9-24 22:54:26 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
    C++ 比 Octave 慢好多,怎么破?
    - u- ]$ E3 @+ q7 z4 U! E# s- }+ t3 c4 T  x9 M2 p3 i) J( B; {) y, p
    自相关两层循环,内层循环涉及浮点数计算,试验了一下把内层循环内部全都 comment out 只留个壳子,  但空的内层循环本身就把速度拉下来了,看来问题并不在浮点计算。4 l; o; R4 R. E3 M- V2 Y

    + ^# q1 H) ?0 q/ W: Y8 F3 k; B+ m速度优化问题真的很有意思啊。
    4 v. u! \1 K# N# h( E" _- l# P, E. D; e  f3 h7 w7 [( r% D
    欢迎大家继续讨论

    评分

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

    查看全部评分

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

    [LV.Master]无

    沙发
    发表于 2022-9-24 23:04:16 | 只看该作者
    拉下来?拉多少?3 \' a2 F- V# U
    把代码贴上来看看?2 u$ h' L8 l3 Y/ B( l

    6 |0 {2 y# l( W. q8 `# }难道分支预测不准破坏流水线执行?不该啊。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    板凳
    发表于 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 编辑
    : ]' J$ w1 o; Y& M  L7 ~: c5 _5 m
    数值分析 发表于 2022-9-24 23:04
    ) a0 @) Y$ y' E拉下来?拉多少?# C2 ~9 X3 O) l+ r" a; @% W/ E7 Y& R
    把代码贴上来看看?

    , h+ ?$ I, E1 l- t8 D8 x
      D/ b  o  U8 Z0 {& Fvoid xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)" N2 R6 P9 j5 b9 t# s, a( N/ R4 t
    {
    . U  E0 s  g" A, R3 d; c; ^5 W& s        comp temp, xtimesy;5 e; b. T  m+ G* l8 Q" l
            xtimesy.re = 0;% o# U$ q( S& c/ J2 \& K6 O
            xtimesy.im = 0;
    : R8 r! ~1 k) M: X5 t* {) R/ O, s        int j0 = lenB - 1;$ X4 [, V. W6 O
            int    i, j, i1, reali;' p0 Z3 W! A+ c/ m+ D/ k: s" D
            if (lenA % 2 == 1): Z  b0 E0 D3 C. Z+ F
                    reali = lenA + 1;
    ! V& N0 j% Q  S* x+ O        else, D/ B  E. Q+ B9 L* X( @5 q
                    reali = lenA;# {6 }( A" p; ]: Y
            reali /= 2;$ z0 c' i" g  o) e& U4 V: ~4 {

    + j' f% a/ S$ \' _8 E5 V; \" P        int nconv = reali + lenB;
    , Y6 ~% {2 A# b3 i; {7 t        //#pragma omp parallel for7 D8 q/ V3 c6 Z
            for (i = reali; i < nconv; i++)
    4 v+ L' H7 C) G# A, D) B        {3 q6 q% h1 K/ d9 v
                    temp.re = 0;
    $ ]1 u& |6 x0 `2 h. d' m4 Y* N                temp.im = 0;
    . y1 v. m0 `( D                i1 = i;6 _, C( n" r* ]1 K8 R
                    for (j = j0; j >= 0; j--)
    - G; @" L0 h& q8 R; _                {' R6 x  e1 M1 Y9 L# a; H
                            /* floating date operation */
    ) g! V8 O3 ^7 P& Z0 ?                }

    , @2 G5 ~$ J  _. w* u& E( U8 ]        }
    4 T1 q8 T' o% r$ l  u. Q}
    6 F+ L/ e. T) e- r4 B3 F; P% x
    ; ?, n$ T- O5 ?, ^xcorr函数代码如上,comp是复数struct, 做过长度为11、19两个矢量的测试,和octave结果完全一样, [4 i/ s* f; O/ L
    ) z" N7 ^: m; z! d0 ]' E
    红色部分是内循环,现在其内部操作都comment out 了, j0大概是 6000。
    ' b+ i6 c* H/ m- Q2 d; d7 |现在call xcorr 100次,耗时78s.
    ( I0 b4 y: u* _) t. h$ o& m' t9 N. R7 G
    如果把红色部分内循环本身完全comment out, call xcorr 1000次,耗时 <1s.
    1 y% O. ~( k+ j0 F  y  q5 b: i1 F2 A
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    6#
     楼主| 发表于 2022-9-25 00:17:17 | 只看该作者
    风雨无阻 发表于 2022-9-24 23:33' Z) [7 f- D$ N7 `; U
    Maybe Debug mode?

    / `# [# `* J& c5 r/ H2 T! }
    ( ~9 E7 C7 v5 s$ p+ Q; c不应该,看我上面的回复。1 i7 A. r, K# [0 F! K% K: c+ B

    ( S7 X" @. n8 ]8 n3 o0 _* x5 j& P- j我更怀疑是 VS 社区版的问题
    回复 支持 反对

    使用道具 举报

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

    [LV.Master]无

    7#
    发表于 2022-9-25 00:20:10 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 00:24 编辑
    " m9 q4 Y- Y3 z" q: p/ G5 J) S9 \
    雷达 发表于 2022-9-24 23:54! ~+ q- T9 k5 o+ H7 ^9 r+ g
    void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)' ]5 F# F: \9 q) Z
    {
    3 s; c6 ?9 \; i        comp temp, xtimesy;
    6 X6 a3 e0 v7 i* ~/ {$ W$ A
    1 @' w; d9 x+ ~& U2 U5 z
    这个不是这么比的吧。。。
    , L  M, A0 t* J1 `0 r2 n
    4 F" `9 j8 M* P6 a# g9 u( ^  C您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。) p/ l2 K: g9 u5 F

    0 A5 v" w4 _) j5 \$ W而加上内循环,光jmp和dec指令就至少多执行了6000个,慢个几十倍不是正常的么?
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    8#
     楼主| 发表于 2022-9-25 00:46:56 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 01:09 编辑
      X) N  c+ |- K" b
    数值分析 发表于 2022-9-25 00:206 ^+ X1 Y, d7 a5 X7 H: W
    这个不是这么比的吧。。。
    0 S6 a" @- f- l9 t/ g3 X+ l0 x8 \' I; q, s. a0 X, r3 j
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。

    7 }4 P1 o) s4 w$ H5 Q! f4 N4 a9 M/ A' P; y6 {6 l. F
    有道理。4 v- D+ u1 @  i+ W  o5 l8 _
    所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,记得 numeric 库里有算向量内积的,我回头试试。) p0 t- L- ~7 |9 Y& \
    8 j3 S& b# c+ ^! C6 c; g
    我先尝试尽量用标准库,一个小程序,不想搞得太复杂。多谢了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    9#
    发表于 2022-9-25 01:27:53 | 只看该作者
    雷达 发表于 2022-9-25 00:46
    / F% N9 Y& k/ b8 B- q- [; p有道理。/ U0 F' p& t; _" z) U5 G( |
    所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...
    " L, _. ~6 g% e  t
    你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子,把循环本身翻译成机器指令loop或dec/jnz,两者速度上会差很多; |: @( d5 N" U2 f! N% l( b
    Why is the loop instruction slow? Couldn't Intel have implemented it efficiently?
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    10#
    发表于 2022-9-25 01:48:51 | 只看该作者
    数值分析 发表于 2022-9-25 00:201 j: r5 U! X8 `
    这个不是这么比的吧。。。7 q. G5 w. }- p+ @
    * j5 _" S' ]! Q# Q) C
    您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
    而加上内循环,光jmp和dec指令就至少多执行了6000个

    0 B/ c0 G7 l% ]% \! w- Z5 }8 w$ W9 u( O) C4 P3 @
    现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果循环这样跑,花不了多少时间。
    回复 支持 反对

    使用道具 举报

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

    [LV.Master]无

    11#
    发表于 2022-9-25 02:06:44 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 02:16 编辑 * f' N% O) D9 e- j. t; S; P
    沉宝 发表于 2022-9-25 01:480 R6 J" L3 W! q2 y% p+ ]
    现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果 ...
    # [% S- ^; E1 a7 K% S: X; u# T, B

    & Z9 @5 Z8 ?5 w/ j; d是的,兄台说的对。" p; o. C/ ?! W' x7 T4 i' x
    / i6 B4 K. t2 ]0 z
    其实我想说的是 真正数值计算部分和代码中其他不直接计算的overhead的比值这个事儿。; E# T9 k/ w& i! R  S! W, K5 O

    ' S' ]8 l, r( A& c$ h6 @/ e6 D雷达兄构造测试用例的时候,屏蔽掉了所有计算的部分,使得剩下的都是overhead,这样run time比较的结果就显得好像不合理了。如果把计算加回去,计算部分的run time会dominate,结果就不那么离谱了。因为不好说,所以用指令数对比的方式试图直观地说明这一点。
    ( v  X, B$ j  w% ?1 V/ C5 o# W6 U; u9 U+ j
    比如说,如果有计算,那么跑六千个循环相对于计算应该用不了多少时间。但是如果一边是什么都不做,另一边是六千个循环,那六千个循环比什么都不做慢几十倍了,就不是那么不合理了。
    * i1 |' b$ s' y9 h- l  e3 C7 M  j
    5 I  p! w+ U1 A2 y5 P$ ?5 G当然也有可能像兄台说的,是优化参数的问题,但我觉得更多地是测试用例设计的不合理。
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    12#
     楼主| 发表于 2022-9-25 04:47:43 | 只看该作者
    本帖最后由 雷达 于 2022-9-25 04:49 编辑 6 d% T. g) Q9 L/ S6 E
    沉宝 发表于 2022-9-25 01:27
    ( D4 p3 j( R+ g4 w: W' ]; w- \你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子 ...

    ; m/ M( H' Y1 L, u2 h' O: n
    : }5 v: A. C% x. U2 l8 C5 D! p又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差一倍,我上面这个差的太多了。
    * a( A% F& B6 g& d8 [; G8 ?) H$ B) ~! T, a; v" S3 e
    我已经完全懵了。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    13#
    发表于 2022-9-25 05:51:23 | 只看该作者
    雷达 发表于 2022-9-25 04:47
    8 p, y8 S) `3 Z% `$ r" R2 x6 }又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...

    ) \, w6 {2 z; q6 A时间差一倍的结果可以接受。
    5 J% d! h/ V8 W$ j: C) x. F* P8 D) K
    你还是用profile工具看看吧。现在大家都主观瞎猜。
    回复 支持 反对

    使用道具 举报

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

    [LV.Master]无

    14#
    发表于 2022-9-25 14:58:59 | 只看该作者
    本帖最后由 数值分析 于 2022-9-25 15:38 编辑
    ' e8 y9 X& ^6 o3 x; L+ S0 P. u  N
    雷达 发表于 2022-9-25 04:47
    + W& I( |' _% X( F3 j" V2 o又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
    ) U7 Y6 N  S0 _$ b7 D

    + d4 k  `; d4 W" S; I$ p
    5 ]* s. L! F* g5 S* f
    ) s( o# \" t' k. z, E能不能把这个也贴上来,看看和上一个有什么不同?
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    15#
     楼主| 发表于 2022-9-26 01:30:03 | 只看该作者
    本帖最后由 雷达 于 2022-9-27 01:17 编辑
    ) E5 [4 ^4 A; p; }6 ^$ m
    数值分析 发表于 2022-9-25 14:58$ ]) `6 o% F1 d2 K
    能不能把这个也贴上来,看看和上一个有什么不同?
    ) l3 l: y+ I- e! x3 B9 o# k
    理了理思路,重新做了一个测试。
    9 B# T, T5 P' }# D; c) Y* q做了两个 vector 和 两个 float *, 都长 100000* X, d  @% d$ t  ?- b) o* j3 C' l1 {
    外循环 6000,里面先做随机数生成,模拟真实环境,避免数据的 cache.
    $ R  I" Q' |/ [8 \. i' a$ h
    " c( c0 R2 `( I内循环试了4种方法,
    , \3 A3 `  ?" _1 X/ \& D1. 直接调用 vector inner_product 247s / a2 B& N# x9 I5 `2 v
    2. vector 循环点乘累加 237s5 H) b1 E5 W+ x9 N9 o
    3. float * 循环点乘累加 204s3 ^) D- [, }/ L/ l9 X: S% h
    4. 空循环 100000 次 202s. c2 x4 h- W. m. r
    9 U" ~6 ]) T9 N& a2 |9 x/ \) U  j
    不做内循环 200s; o+ P9 q3 c: F  i5 d6 T% l3 ^
    9 }- d$ `6 E, B0 W; Q% S
    你昨天说的对,内循环本身占比是很小的,大头在其他处理。
    % b6 H2 L! b" M' Y另外可以看到, float * 循环点乘累加 并不差,比用vector 还更快。' p  G# R4 V: G) r
    ' a- L# d+ I9 c8 Z
    至于我那个原始程序,还有一些疑问,见5楼,其他都不变仅仅是有无空的内循环就有很大不同,这是不对的,也许有一些其他缺陷我没有看到。(也许可以改成 while 试试)4 c5 H4 p+ M. C1 y) G

    $ E1 B2 R" D3 `0 @5 P(为什么下面我贴的  b1 加 方括号里的 i , 显示出来却是 b1 ?方括号 i 消失了。 LOL . 改成  jj 好了,原来 方括号里的 i 是斜体标志  LOL)
    2 }( H4 `$ s4 m0 V# P9 i: j8 e+ \8 ~  S- w
            std::vector < float > vec1(N);7 o' L) _+ Q  l8 {0 h8 S2 I: `
            std::vector < float > vec2(N);
    ! t7 ?3 k1 |& }- \! A  H        float* b1 = new float[N];
    . E  E8 g; `( e3 ~" ~* t: t* n        float* b2 = new float[N];
    ! W1 Y# w, J6 C7 `/ ~) R. n* E  X6 V
            for (int j = 0; j < 6000; j++), n+ h2 `: v" E9 i5 p
            {# ]$ K) J9 N8 @* x* r- E% N
                    std::generate(vec1.begin(), vec1.end(), []() {
    2 V" U+ {% p2 n0 A                        return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;$ m' _% X2 e1 o4 e
                            });% r' K' Y7 {  s) T

    ( n# ^9 k; @! P% h4 O+ O1 w6 B% M$ [                std::generate(vec2.begin(), vec2.end(), []() {! a! a* o. M( j* p1 ~3 Q
                            return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
    / g7 H, j2 G* A! S9 |+ F8 p                        });  e" Z0 X  l) v4 C. A/ p( A

    " H0 o# y5 o7 y                for (size_t jj = 0; jj < vec1.size(); jj++)
      q- A" F/ ^% W% u/ @                {7 @% D. v* J3 J2 t
                            b1[jj] = vec1[jj];
    8 E5 {' L0 Q8 a+ i& J6 s+ X                }
    % M: k3 t3 o( _5 G& A: P1 ?0 P7 ?9 P! M5 I4 E
                    for (size_t jj = 0; jj < vec2.size(); jj++)
    & f7 J4 n. J, D- [; i& b( B5 _                {# X8 e$ W0 Y3 R2 v+ g9 }
                            b2[jj] = vec2[jj];% h) r5 j: y1 q& R0 ?6 d
                    }9 a8 T7 H/ }0 [% e' f0 s' v1 I
    , [; L2 U3 {) Y
                    //Method - 1  N=100000 247s  
    1 G- Z. u" |) x: U: p. V) v) P                //fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);; R: T2 r5 N' X. p2 m' ]6 I& c
                                    * B! f2 F, J% C9 ^, l% e% H
                    //Method - 2  N=100000  237s3 M8 N" R& Y4 _" }2 g8 d0 [
                    /*
    ) _( Z9 |  R( P; z, Z; H                for (int jj = 0; jj < N ; jj++)5 e% ?  R9 S2 h& H! ]1 v# S1 A
                    {2 [5 R- V- C0 u9 a) `
                            fresult += vec1[jj] * vec2[jj];9 m9 |( i3 g! z. u1 @$ Q0 V) v
                    }
    * U" b$ n! j3 g8 I/ R                */2 ?" U) E) }; y
                                    ; n) E0 [* v- E) y+ K  t, s( a
                    //Method - 3  N=100000 204s9 ^3 w; `' P. S; d6 y- {; M/ ]
                    /*, ~# r/ @- k9 \- S$ l
                    for (int jj = 0; jj < N; jj++)
    6 L4 s' s4 g! w( T                {
    8 S' }4 y, J- _5 w* b$ a' {+ E                        fresult += b1[jj] * b2[jj];
    # ~2 r( \# k7 _7 b% p( d, S5 [3 Z                }
    ) W5 G6 ]  m5 {7 ^% l                */
    8 `3 K; I8 N4 \- t4 A: p" o
    0 T, T0 E! Y: E- Y! }                //Method - 4   202s2 n" |. N8 M! q7 K% P
                    /*# T1 a8 Y# t, s) t- Z
                    for (int jj = 0; jj < N; jj++)
    % `, L! U  {( Y6 W3 L: r" g. e                {
    7 t6 F3 y" i, D5 v) y4 l. M* H                        
    0 E4 {) S3 d# t+ ]# a. V% b3 n                }" X5 f' E9 p4 K* u2 P) _1 ~
                    */
    5 s4 B" [% D* u! P                //comment out all methods, N=100000  202s               
    $ Q5 E5 Y+ x0 [& `7 [        }1 o, G0 {; H" _, B0 L& @3 [" e
    & l- W5 c5 _, n" m. d
            delete []b1;
    / P, }$ S3 T% Y) o2 b# E# }/ p        delete []b2;

    & `. ^: _. @" \
    回复 支持 反对

    使用道具 举报

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

    [LV.1]炼气

    16#
    发表于 2022-9-27 00:15:38 | 只看该作者
    瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?9 a' F  R& p8 q4 h. f" V# ~

    ! K+ S5 J3 ~* i  c- w你第二个试验里面的j在循环里面又重新定义了啊,你确定真的跑了6000次?9 \( i4 Z5 R& i. ]) p5 }
    回复 支持 反对

    使用道具 举报

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

    [LV.10]大乘

    17#
     楼主| 发表于 2022-9-27 01:16:03 | 只看该作者
    机器猫 发表于 2022-9-27 00:15
    3 G+ O* p& b  N瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?
    , V- |" J" [/ S; Y2 G1 a2 Y2 E( w4 g  X) f8 w7 C
    你第二个试验里面的j在循环里面又重新定义 ...

    / k; F8 {: G2 C$ O内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
    / o' G. u- i1 o3 c
    # ^! J" w0 g& J. r; X不和它较劲了,瞎耽误工夫,我已经转到 ubuntu, 也准备顺便试试 avx2 向量化。
    回复 支持 反对

    使用道具 举报

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

    [LV.1]炼气

    18#
    发表于 2022-9-27 02:06:58 | 只看该作者
    雷达 发表于 2022-9-27 01:16
    1 ^" i! p4 X2 j* M( C& O* e% ^0 z内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL9 M( D6 _  L8 S7 r: _0 {$ H
    & _  t% J3 f* ^* t1 s: `& L) [
    不和它 ...
      ~6 I  X/ H6 Q: t6 R# `8 C
    ! }+ d9 Q( ~% t1 }3 i
    不过可以试试我说的register变量。前一个试验j是混在一堆其它变量里一起定义的,很有可能是在stack上,这样内存读写会更多,要是再碰上每次都需要加载cache就更慢了。4 ]7 R1 q$ U7 ^! K
    后面一个是在循环那里定义的,说不定编译器就把它优化成register变量了
    回复 支持 反对

    使用道具 举报

    该用户从未签到

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

    使用道具 举报

    该用户从未签到

    20#
    发表于 2022-9-27 20:29:50 | 只看该作者
    雷达 发表于 2022-9-24 23:54% G6 r( w! q* I" B' j/ s; r' V3 B
    void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
    % `% H& F3 B/ ?" n) L{
    ' e& s5 ^* Z, P) t' c        comp temp, xtimesy;

    . D; w: `# P& b3 E5 ~+ h) h  v这个code里面如果Openmp没有被注释掉的话,那么temp那个变量应该是定义在循环里面,否则线程之间会存在争夺写入那个temp的风险。. p8 K. O# u6 }6 ]& X6 `. b& h# T
    内层for循环如果没有内部操作的话,编译时应该被优化掉了,和你完全注册掉整个循环是一回事。可能你的编译设置没有打开优化?. I' X+ ~2 K8 ~2 w; R
    VS社区版没有问题,我工作用的就是社区版,设置正常的话不会比商业版差。以前游说头头用Intel Compiler,他说不想花钱,而且差不了多少,就一直用到现在。
    回复 支持 反对

    使用道具 举报

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

    GMT+8, 2025-9-21 17:17 , Processed in 0.051180 second(s), 19 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

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