数值分析 发表于 2022-9-24 23:04# {) T- V7 V$ J5 N2 B
拉下来？拉多少？
把代码贴上来看看？
风雨无阻 发表于 2022-9-24 23:33" O* i/ G: r5 s! p+ m
Maybe Debug mode?
雷达 发表于 2022-9-24 23:548 [: h8 u$ k3 @# b# ~+ B0 Z
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB); v+ B& w4 |: o9 p; ~. o! U) m( x- p
{
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。1 y; w- Y& L7 d
. I. f0 C; E7 H
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到100个）。
雷达 发表于 2022-9-25 00:46# g# ^( `) l: Q
有道理。- M* s% d$ K; k! o9 T
所以存在内循环速度就上不去，把内循环取消，改成两个向量直接点乘再求和应该就会好得多，这大 ...
数值分析 发表于 2022-9-25 00:20# z& v3 E. o- c6 N# }# h. J
这个不是这么比的吧。。。
- X y! e3 e3 n3 E3 z
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到100个）。
而加上内循环，光jmp和dec指令就至少多执行了6000个
沉宝 发表于 2022-9-25 01:48! T% V, E$ v$ A8 l
现在的CPU，可以把判断、jmp和dec指令全部融合进一个µOp（微操作，CPU内部流水线上的执行单位）。如果 ...
沉宝 发表于 2022-9-25 01:27) P4 T9 ]$ D3 u
你两个试验之间就差了一个空循环， call 1000次按理不会有秒级差异，可能还是编译器优化的问题。举个例子 ...
雷达 发表于 2022-9-25 04:472 M# t& c7 L) h4 z& h
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
雷达 发表于 2022-9-25 04:474 T$ l! Y) D2 T4 ~+ z, B
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
数值分析 发表于 2022-9-25 14:58) P. E/ B3 M4 d' v- h' c; Y. V
能不能把这个也贴上来，看看和上一个有什么不同？
std::vector < float > vec1(N);
std::vector < float > vec2(N);4 z- u5 j, C4 @0 w' e2 H
float* b1 = new float[N];
float* b2 = new float[N];
for (int j = 0; j < 6000; j++). T# _3 M/ D$ B- X" H9 d0 g& }
{% `' V/ ]+ x4 S; r/ v1 V
std::generate(vec1.begin(), vec1.end(), []() {' O3 Q, n) y3 x7 t% R
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
});
std::generate(vec2.begin(), vec2.end(), []() {
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;
});
" ]* Y% }5 A \/ W3 Q) p$ a }
for (size_t jj = 0; jj < vec1.size(); jj++)' J( f% ~: z. k) @
{4 d6 B! s. l( C6 g4 s; S/ o
b1[jj] = vec1[jj];- V5 L& J+ |' H' c: o, I
}
for (size_t jj = 0; jj < vec2.size(); jj++)7 R* D: Q* L; R3 a
{2 I4 a6 W* s% N" ?! y: m L3 w
b2[jj] = vec2[jj];/ I/ v' w8 O/ S2 T3 Q. ~$ A/ N
}% g' C# h3 x2 u% [
//Method - 1 N=100000 247s $ p0 _2 I6 }! f3 t
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0); K) X1 r+ X$ j. s9 D
* S8 u! }* y5 |! ^8 T
//Method - 2 N=100000 237s
/*
for (int jj = 0; jj < N ; jj++): j( p2 y* |+ `/ C, @. I
{
fresult += vec1[jj] * vec2[jj];
}' u; x y, o5 m
*/
2 ?2 T! g& K. A! T0 y- g9 a9 h$ u
//Method - 3 N=100000 204s
/*
for (int jj = 0; jj < N; jj++)
{8 F4 r0 z X# M4 ^1 C1 w
fresult += b1[jj] * b2[jj];
}* b T C2 M2 \- l$ E
*/- ]) W* \3 o' c. W+ O3 D
. e! Z; ?; N& K6 _
//Method - 4 202s
/*
for (int jj = 0; jj < N; jj++)
{4 T) k/ u5 P8 r D
}0 j2 @8 Z* s8 g/ W* x
*/
//comment out all methods, N=100000 202s 2 ]. m& L/ H7 O' U! q6 e2 N
}
) U0 a# W, H) P5 o% ~- F
delete []b1;
delete []b2;
机器猫 发表于 2022-9-27 00:15( ~3 L$ q$ h5 S
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同？
: h" a1 k* p6 c: }2 P, d0 C
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:16
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
. k- ` f: N; Y) Z6 g7 F& k0 F
不和它 ...
2 N5 [0 g2 ]% E, @雷达 发表于 2022-9-24 23:548 { J- S R D; o
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{# a+ ^9 L3 ?; m2 `5 A K
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:30$ i. J3 G+ o, m$ o
理了理思路，重新做了一个测试。
做了两个 vector 和 两个 float *, 都长 1000003 }: R2 I3 z7 X1 G3 J4 S
外循环 6000，里面先做随 ...
雷声 发表于 2022-9-27 20:39
这个时间是从哪里开始算的？" O4 a& x; s# x$ N1 C1 R- [
我怀疑这个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& Q2 r- Y* c5 n
我有些迷糊，这样的code，难道不就应该时间差很多吗？也做了个简单的实验，你看看我做的有错吗
/ `, e' }: O$ k% M( Q
...
。
8 s% |0 E2 n; o( F$ m| 欢迎光临 爱吱声 (http://129.226.69.186/bbs/) | Powered by Discuz! X3.2 |