数值分析 发表于 2022-9-24 23:04
拉下来?拉多少?0 u; O- g j- x6 e. }
把代码贴上来看看?
风雨无阻 发表于 2022-9-24 23:33
Maybe Debug mode?
雷达 发表于 2022-9-24 23:54* v8 D7 e1 E5 D
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)3 n8 B/ o z9 n# l& t
{6 b+ _, v. ^" I
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:20! |5 J( t$ q( }% o; n
这个不是这么比的吧。。。
: V( r. ?% ?) ^$ G4 t
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
雷达 发表于 2022-9-25 00:46
有道理。" t- G- e' ?( o
所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...
数值分析 发表于 2022-9-25 00:209 ~ I" l4 m& A3 @
这个不是这么比的吧。。。) g9 _& z, ^; e* f$ y
( p& H& Y- [4 |+ n: c
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
而加上内循环,光jmp和dec指令就至少多执行了6000个
沉宝 发表于 2022-9-25 01:481 d0 j( P0 e$ t! [4 P
现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果 ...
沉宝 发表于 2022-9-25 01:27" w/ f, t3 l6 u3 s' i7 ~1 N2 F
你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子 ...
雷达 发表于 2022-9-25 04:473 g" k X9 o6 b# o! w! ^
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
雷达 发表于 2022-9-25 04:47# J" T4 D+ m. |1 o L6 x
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
数值分析 发表于 2022-9-25 14:588 H1 z: ]* `7 U8 B
能不能把这个也贴上来,看看和上一个有什么不同?
std::vector < float > vec1(N);/ a$ b+ V& k& F5 D8 b
std::vector < float > vec2(N);
float* b1 = new float[N];# Z' E% Q: M0 j. F# }9 C
float* b2 = new float[N];
7 H. y/ c5 p8 {+ L& R
for (int j = 0; j < 6000; j++)
{3 g- k) a! n1 r7 H; H5 \
std::generate(vec1.begin(), vec1.end(), []() {
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
});, t( _( q/ {$ h, i
* X( s$ L0 [9 f, O& x& y
std::generate(vec2.begin(), vec2.end(), []() {
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
});: N$ o# L8 K( G3 q1 C2 n# @
" i) U7 P! [* F( Y# n! ^' T
for (size_t jj = 0; jj < vec1.size(); jj++)
{0 U9 o2 }* H. M
b1[jj] = vec1[jj];. w! G0 C1 _. L+ G4 }3 c
}* I) s3 ~# N8 S( J$ F
for (size_t jj = 0; jj < vec2.size(); jj++)
{ T5 _8 K( h& m: {
b2[jj] = vec2[jj];
}7 n" b; S5 `* U8 g. T
9 g9 n6 d7 s& \% i6 p! v7 y
//Method - 1 N=100000 247s
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);
//Method - 2 N=100000 237s
/*
for (int jj = 0; jj < N ; jj++)- I! H3 L7 p m1 b! P' d
{( I. x) Z! u" _ O
fresult += vec1[jj] * vec2[jj];7 w9 j5 \2 }& V* w4 F: @
}. o6 p$ o& Y+ [( p* n% L7 n
*/
//Method - 3 N=100000 204s
/*3 H J6 V) d; |" J
for (int jj = 0; jj < N; jj++); B, A' v! k8 s
{6 N. q4 z* r, G2 k
fresult += b1[jj] * b2[jj];
}0 Q0 \" _7 f/ K
*/: B$ d G" H6 B5 I- x4 N* G
//Method - 4 202s, X* [8 |; L" T7 s
/*
for (int jj = 0; jj < N; jj++)
{
}* d& |! V& X4 q4 i
*/
//comment out all methods, N=100000 202s # z" p2 z# c; ~ p
}
( R1 O: U- k2 p! E) f/ y
delete []b1;
delete []b2;
机器猫 发表于 2022-9-27 00:15/ D8 h" U5 |, |% s
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?+ `1 M Q* y6 @- P- C
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:164 O% M( c& G8 Y; ]1 c X0 N }
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
不和它 ...
/ y: J6 H; Z8 L- ]1 X* k雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB), V6 _( G0 z3 z& ]
{
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:30" ~/ a3 G; @$ E& a3 c. s
理了理思路,重新做了一个测试。6 `2 q0 F$ i7 t' T* _
做了两个 vector 和 两个 float *, 都长 100000* ^$ [) D" ^) d, b
外循环 6000,里面先做随 ...
雷声 发表于 2022-9-27 20:39) O# G# j) k# ]5 ^9 ?
这个时间是从哪里开始算的?
我怀疑这个200多秒里面有200秒花在产生随机数上了,真正计算大概只用了2秒, ...
雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{
comp temp, xtimesy;
opensrc 发表于 2022-9-28 00:29
我有些迷糊,这样的code,难道不就应该时间差很多吗?也做了个简单的实验,你看看我做的有错吗0 A8 ^# I9 {3 O7 {# Z
, V1 G/ S. W$ L5 n" U* V9 J
...
。% M# i4 g2 X9 K' c

3 r. @- y, K2 n5 r+ I6 V| 欢迎光临 爱吱声 (http://129.226.69.186/bbs/) | Powered by Discuz! X3.2 |