数值分析 发表于 2022-9-24 23:042 P, f- N4 d9 G% i! P
拉下来?拉多少?$ ^8 S6 b% S, }& W& J
把代码贴上来看看?
风雨无阻 发表于 2022-9-24 23:33! B% y* E+ h+ Y9 {* k
Maybe Debug mode?
雷达 发表于 2022-9-24 23:544 V6 c7 d4 R; f
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。
! b! o6 S0 ^/ G$ r3 e, Z7 T
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到100个)。
雷达 发表于 2022-9-25 00:46. h1 X' W8 ~ b( Q4 {. ~
有道理。
所以存在内循环速度就上不去,把内循环取消,改成两个向量直接点乘再求和应该就会好得多,这大 ...
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。$ k. E% U$ o8 e2 g
您这个函数,不带内循环的话,汇编完总共操作也没几个(不到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:47
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
雷达 发表于 2022-9-25 04:47
又写了个小实验,没有调用子函数,双层循环,外层6千次,内循环30万次空转,有或没有空转内循环,时间差 ...
数值分析 发表于 2022-9-25 14:58
能不能把这个也贴上来,看看和上一个有什么不同?
std::vector < float > vec1(N);
std::vector < float > vec2(N);# X- f5 R- k" k" y: J
float* b1 = new float[N];' E" k1 W& D" @0 f! z
float* b2 = new float[N];
for (int j = 0; j < 6000; j++)
{. w5 E( K* M& G5 Z
std::generate(vec1.begin(), vec1.end(), []() {/ x. I# ]$ _/ N
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;- b9 Z; K: }: q, O0 C
});+ w, D, b% I# d- T) J& n
std::generate(vec2.begin(), vec2.end(), []() {) _; h# F7 X2 y! [; Y; C0 `
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;9 ^- [; C8 L& e" e: L( _
});
2 J s4 u" z* R0 y
for (size_t jj = 0; jj < vec1.size(); jj++)
{" ]7 m% B) X) I- Y$ h, d% B
b1[jj] = vec1[jj];3 v8 a" k0 r! @# L; o. M4 f M
}
for (size_t jj = 0; jj < vec2.size(); jj++)
{) f4 k7 A# Y/ x6 M8 x
b2[jj] = vec2[jj];
}
0 a# s$ |; Y" z3 O7 y- V
//Method - 1 N=100000 247s
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);
//Method - 2 N=100000 237s
/*: F& p9 i' E2 K6 M8 w3 j
for (int jj = 0; jj < N ; jj++)( g% v2 X% o, L- L0 k% E
{, }% r" _' w: y/ H3 l; m- R, x
fresult += vec1[jj] * vec2[jj];0 V4 D7 t& J/ R2 }5 J+ Z8 ~
}
*/- Q( ^4 G0 I7 p1 A2 s4 S/ o
6 ^3 ?, L7 s5 X- L5 U+ }3 s
//Method - 3 N=100000 204s
/*
for (int jj = 0; jj < N; jj++)' h0 z7 K- X7 O! K
{5 s& x k0 u( y: S8 N6 U1 P
fresult += b1[jj] * b2[jj];, v) b7 [% G+ N
}
*/
* W0 l$ |$ W& I# b
//Method - 4 202s: c! F. F) K, f. e' f J8 V! f M4 Z
/*
for (int jj = 0; jj < N; jj++)% c4 A- T$ L& s, L ^
{
}* q7 V. ]7 m: {" J3 a+ V, G; D- k
*/
//comment out all methods, N=100000 202s ; Q( j1 s" j; n# ]1 ]
}' h3 Y* z9 |, k( D" a
4 K; y6 b; N) O) R6 `
delete []b1;
delete []b2;
机器猫 发表于 2022-9-27 00:15
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同?
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:16# O8 N0 N& p- K$ g& }% n% _
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL- e6 X X I3 _6 h
& x! S& F8 \: M w& ^: `6 \- t
不和它 ...

雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:30; O5 E! Q7 g/ y' v
理了理思路,重新做了一个测试。4 F3 [8 l! [- Z
做了两个 vector 和 两个 float *, 都长 100000
外循环 6000,里面先做随 ...
雷声 发表于 2022-9-27 20:39
这个时间是从哪里开始算的?9 u- B6 L9 q5 o% Y/ g `. @
我怀疑这个200多秒里面有200秒花在产生随机数上了,真正计算大概只用了2秒, ...
雷达 发表于 2022-9-24 23:54$ E4 F7 s$ Q- A# ]8 ` u% W# o
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)$ _% `. p- }/ W5 V
{ {! ^6 O8 C* w
comp temp, xtimesy;
opensrc 发表于 2022-9-28 00:29
我有些迷糊,这样的code,难道不就应该时间差很多吗?也做了个简单的实验,你看看我做的有错吗
...
。

! t3 z5 T7 _1 v7 _2 g# e| 欢迎光临 爱吱声 (http://www.aswetalk.net/bbs/) | Powered by Discuz! X3.2 |