数值分析 发表于 2022-9-24 23:040 b+ u8 M6 {, k* S0 l. q1 J
拉下来？拉多少？
把代码贴上来看看？
风雨无阻 发表于 2022-9-24 23:33
Maybe Debug mode?
雷达 发表于 2022-9-24 23:54& O7 i9 N( p% O, C; j1 I" P& X
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。 b; `4 y; I3 n- z4 v
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到100个）。
雷达 发表于 2022-9-25 00:46
有道理。$ w2 S, K; \8 G, T2 q4 ~
所以存在内循环速度就上不去，把内循环取消，改成两个向量直接点乘再求和应该就会好得多，这大 ...
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。2 O# Y9 \7 {# d6 {1 G
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到100个）。
而加上内循环，光jmp和dec指令就至少多执行了6000个
沉宝 发表于 2022-9-25 01:48" Y8 H9 Z0 D. A* V
现在的CPU，可以把判断、jmp和dec指令全部融合进一个µOp（微操作，CPU内部流水线上的执行单位）。如果 ...
沉宝 发表于 2022-9-25 01:27" Y1 P% Q: }; k) g
你两个试验之间就差了一个空循环， call 1000次按理不会有秒级差异，可能还是编译器优化的问题。举个例子 ...
雷达 发表于 2022-9-25 04:47, d) b5 J, T7 J
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
雷达 发表于 2022-9-25 04:47* e, [3 v; A) l3 |
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
数值分析 发表于 2022-9-25 14:58. _- p6 g: g1 R; G2 }1 f6 s
能不能把这个也贴上来，看看和上一个有什么不同？
std::vector < float > vec1(N);
std::vector < float > vec2(N);# C$ N6 h) v: E" I* b( h
float* b1 = new float[N];
float* b2 = new float[N];# f$ J$ U2 c! }: W" e) R9 ]
for (int j = 0; j < 6000; j++)
{
std::generate(vec1.begin(), vec1.end(), []() {
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
});
std::generate(vec2.begin(), vec2.end(), []() {8 g, ]- _5 Z+ B$ x
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
});
for (size_t jj = 0; jj < vec1.size(); jj++)
{$ d& h9 J( {0 w5 L" s6 [& y! l
b1[jj] = vec1[jj];
}
+ ?+ N2 Z) B; E) c5 A
for (size_t jj = 0; jj < vec2.size(); jj++)
{) ^& C4 w; g2 z: X2 L7 O: p+ B
b2[jj] = vec2[jj];
}. O- ^ u, {9 [. w) }- X
//Method - 1 N=100000 247s 7 {4 I7 k( j3 E7 `$ V
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);! ?% v# h% p# z
8 z$ _; k" M1 L# f g
//Method - 2 N=100000 237s; c1 D( g! e1 R1 K9 T0 D
/*$ ]# [- |: `; @' q
for (int jj = 0; jj < N ; jj++)
{
fresult += vec1[jj] * vec2[jj];! n3 o7 i( K; R
}
*/ X ^! f! W1 ]: B. S4 Y2 S
//Method - 3 N=100000 204s, ?" S0 D5 M) S9 [5 w. G+ S
/* ?6 z2 B B0 S5 D( s* K; T
for (int jj = 0; jj < N; jj++)0 ^4 A" U& ?1 Q: x g$ o6 ?* p
{
fresult += b1[jj] * b2[jj];
}
*/
) [- H0 L2 X4 J* [# {# u
//Method - 4 202s
/*2 S6 A+ r% b3 q! [. p* v
for (int jj = 0; jj < N; jj++)1 d$ @/ z) {' d9 o6 Y A% d
{2 x+ c' L" _, s& S( B! H: s
}% U( j/ W' L: u) u3 g
*/
//comment out all methods, N=100000 202s * O8 ]* x* {/ t0 O
}) c5 n) s2 N+ k( n5 }5 x
delete []b1;
delete []b2;
机器猫 发表于 2022-9-27 00:15
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同？
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:167 \" ~# y, B7 O3 U" Y* i
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
不和它 ...
! e, J( a7 L$ i2 r! m雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)8 J6 u8 }: b' g! C
{0 o2 t5 k$ C) \5 a* k
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:30/ ^/ W3 K5 _5 U/ c
理了理思路，重新做了一个测试。
做了两个 vector 和 两个 float *, 都长 100000* S0 l. V" \; V c" L" A" `: u
外循环 6000，里面先做随 ...
雷声 发表于 2022-9-27 20:39; m2 l6 `4 i7 m3 R$ ]+ b* `
这个时间是从哪里开始算的？1 n! c- Q9 K0 D
我怀疑这个200多秒里面有200秒花在产生随机数上了，真正计算大概只用了2秒， ...
雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)8 y% l8 P9 R; ^: F- C9 H7 P- t
{
comp temp, xtimesy;
opensrc 发表于 2022-9-28 00:29 V2 R- k8 F5 l0 l. M# x
我有些迷糊，这样的code，难道不就应该时间差很多吗？也做了个简单的实验，你看看我做的有错吗
...
。1 Y/ |) j! S; _6 @/ [" {" |5 _
2 T* [' y4 D! S| 欢迎光临 爱吱声 (http://www.aswetalk.net/bbs/) | Powered by Discuz! X3.2 |