数值分析 发表于 2022-9-24 23:04. k. U& |" d3 A8 _# [# L
拉下来?拉多少?
把代码贴上来看看?
风雨无阻 发表于 2022-9-24 23:33
Maybe Debug mode?
雷达 发表于 2022-9-24 23:54: b a; \7 \8 L5 R% W
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB) L' F5 | w* y J/ k k
{
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:20% x4 }- ?: ^) u$ n) |( U4 O' M
这个不是这么比的吧。。。) s4 G3 W0 C1 u7 T# k' O3 r
m; |% A; g& \1 Q3 b: n w
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
雷达 发表于 2022-9-25 00:463 M+ K. A8 v+ f b
有道理。
所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...
数值分析 发表于 2022-9-25 00:207 j* ]- g7 ]( g% G0 r% d6 k) ~6 [
这个不是这么比的吧。。。
5 K$ x. g" _ ~0 s" w
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
而加上内循环,光jmp和dec指令就至少多执行了6000个
沉宝 发表于 2022-9-25 01:48
现在的CPU,可以把判断、jmp和dec指令全部融合进一个µOp(微操作,CPU内部流水线上的执行单位)。如果 ...
沉宝 发表于 2022-9-25 01:27
你两个试验之间就差了一个空循环, call 1000次按理不会有秒级差异,可能还是编译器优化的问题。举个例子 ...
雷达 发表于 2022-9-25 04:473 U3 \' F* E5 ]0 A7 p1 B! i
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
雷达 发表于 2022-9-25 04:47
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
数值分析 发表于 2022-9-25 14:58, r- }, f8 [# ^4 I1 U
能不能把这个也贴上来,看看和上一个有什么不同?
std::vector < float > vec1(N);. N/ a' [: h0 V+ a
std::vector < float > vec2(N);
float* b1 = new float[N];! [ S8 f/ q$ ]1 ~
float* b2 = new float[N];
, p% b+ m* a2 {, h3 C$ u4 E
for (int j = 0; j < 6000; j++)
{; d) ~* V: x# _1 N/ A
std::generate(vec1.begin(), vec1.end(), []() {2 y t' U9 z# E5 j, ] f
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
});/ [1 W- c. h: T# y
std::generate(vec2.begin(), vec2.end(), []() {: L; c% x" q' z7 x
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
});5 K* d5 m1 o4 D* ?) d- l! f0 y/ i6 Z- t
% m: J, u$ N8 N( a5 W7 a" U" [
for (size_t jj = 0; jj < vec1.size(); jj++)
{
b1[jj] = vec1[jj];
}4 r! ]* b% M3 {" I. y: ]$ O4 _- W
9 f4 n6 i0 v+ @6 t: C! [/ p
for (size_t jj = 0; jj < vec2.size(); jj++)
{1 c8 z. c, p" A) B
b2[jj] = vec2[jj];
}
( k& o7 |8 q( |2 {- ?2 Z
//Method - 1 N=100000 247s 9 b( L# l$ K0 O. Y( @$ t
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);
: w; r% E1 t1 j3 d2 x! O
//Method - 2 N=100000 237s
/*2 b2 Q5 h' |! O$ [* J* K8 D
for (int jj = 0; jj < N ; jj++)
{* n0 c: a8 Y' `7 @; m
fresult += vec1[jj] * vec2[jj];, y+ v7 B+ D! i9 y4 m
}
*/( @ U( K% O( h
@: t; w: e+ h# w. A
//Method - 3 N=100000 204s
/*
for (int jj = 0; jj < N; jj++)! S. R7 R2 P" l; I! K5 v
{
fresult += b1[jj] * b2[jj];
}
*/6 s' i2 W z4 Q% c/ H8 i/ I
3 E0 r) _; H) l s4 O* a+ G
//Method - 4 202s( J0 O2 h) V4 M
/*
for (int jj = 0; jj < N; jj++)
{
2 C4 l4 Z& T0 q0 H) C: b( [
}: t* m0 \. z+ u# F, f1 ~! B
*// n) [5 E9 b- G$ }. K0 _* [
//comment out all methods, N=100000 202s
}
; z% f4 A6 s( u! F" H6 l* ~
delete []b1;$ o; C$ l. a, ^8 ]! K
delete []b2;
机器猫 发表于 2022-9-27 00:15
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:16% q d. z* w( f9 |4 l L
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
& k# K" |' H2 k+ G) E3 W
不和它 ...
* x' e8 M( O' `5 S, v雷达 发表于 2022-9-24 23:54' u- E H Y& T
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{% G2 k' ^9 z- g
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:304 V. `( {* P) Z- I2 q! V0 A- {+ ?
理了理思路,重新做了一个测试。6 @ Z4 R; j k. {4 W
做了两个 vector 和 两个 float *, 都长 100000
外循环 6000,里面先做随 ...
雷声 发表于 2022-9-27 20:39
这个时间是从哪里开始算的?0 K) ^9 b- M7 U
我怀疑这个200多秒里面有200秒花在产生随机数上了,真正计算大概只用了2秒, ...
雷达 发表于 2022-9-24 23:547 _4 Z' i1 J* b& T
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{3 n: [* A* \ m- @ @/ U+ l
comp temp, xtimesy;
opensrc 发表于 2022-9-28 00:29# t) K2 P$ G3 }6 ^
我有些迷糊,这样的code,难道不就应该时间差很多吗?也做了个简单的实验,你看看我做的有错吗
1 y7 i/ E: Y6 M
...
。


| 欢迎光临 爱吱声 (http://129.226.69.186/bbs/) | Powered by Discuz! X3.2 |