数值分析 发表于 2022-9-24 23:04
拉下来？拉多少？
把代码贴上来看看？
风雨无阻 发表于 2022-9-24 23:33
Maybe Debug mode?
雷达 发表于 2022-9-24 23:540 A# A; W3 n$ H2 I; @2 E3 A
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)
{. e) e+ s7 v. i* r( H6 N% p+ M
comp temp, xtimesy;
数值分析 发表于 2022-9-25 00:203 t, h1 S$ n3 [3 O1 ~( C4 a+ p4 n
这个不是这么比的吧。。。( d! v9 K1 {( V( Z/ M# ?
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到100个）。
雷达 发表于 2022-9-25 00:469 `. @' t# A1 e' O& g q
有道理。6 T" p3 I. `# a/ g0 X
所以存在内循环速度就上不去，把内循环取消，改成两个向量直接点乘再求和应该就会好得多，这大 ...
数值分析 发表于 2022-9-25 00:20
这个不是这么比的吧。。。
您这个函数，不带内循环的话，汇编完总共操作也没几个（不到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 a2 g( ~) _6 V" K+ D! r
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
雷达 发表于 2022-9-25 04:47
又写了个小实验，没有调用子函数，双层循环，外层6千次，内循环30万次空转，有或没有空转内循环，时间差 ...
数值分析 发表于 2022-9-25 14:58
能不能把这个也贴上来，看看和上一个有什么不同？
std::vector < float > vec1(N);
std::vector < float > vec2(N);
float* b1 = new float[N];
float* b2 = new float[N];
for (int j = 0; j < 6000; j++)
{- n" Z. x* `" X
std::generate(vec1.begin(), vec1.end(), []() {
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 23.23));;
});
- u. `1 |( i8 X2 i6 g* g( I
std::generate(vec2.begin(), vec2.end(), []() {; o# o( W( B$ ?8 z D2 d; v
return static_cast <float> (rand()) / (static_cast <float> (RAND_MAX / 24.31));;% `2 l' b$ Y+ ?" n: M/ h
});
for (size_t jj = 0; jj < vec1.size(); jj++)( ` |9 I: v9 T
{
b1[jj] = vec1[jj];7 I, |8 R7 B r! w" ]4 I
}
6 \( ?0 Q- G$ w @. l
for (size_t jj = 0; jj < vec2.size(); jj++) c: a) }8 @5 O! K5 { }( j* k$ p
{
b2[jj] = vec2[jj];. ?3 N! Q: ^6 x- N3 y2 s0 f" {
}2 Q2 R' C$ e+ j/ A
; Y. d$ \+ R) l; R/ r9 l
//Method - 1 N=100000 247s ' A1 k/ g# n, v& x
//fresult = inner_product(vec1.begin(), vec1.end(), vec2.begin(), 0);
; C4 D' e# }9 `6 x) R
//Method - 2 N=100000 237s
/*
for (int jj = 0; jj < N ; jj++)
{
fresult += vec1[jj] * vec2[jj];% ` J+ I4 E! u+ H( Z0 x/ E
}
*/) h& J% z: U s" s3 v
//Method - 3 N=100000 204s* {5 d0 d( l' A m3 _1 X
/*1 z& \; X9 a3 n+ _& A) x$ P8 F
for (int jj = 0; jj < N; jj++)5 L; T* s: ~6 o, `) c
{& p2 ^! Y1 G1 S9 S- m$ t
fresult += b1[jj] * b2[jj];% k0 h5 `- y. V
}& L) v. f$ f9 e: p/ V2 S
*/6 @ E, c0 E$ C! |
2 F% q, }8 ?) B$ t. }' {3 Q3 K+ }
//Method - 4 202s
/*
for (int jj = 0; jj < N; jj++)# @" }7 R" h/ w' o6 N7 s
{
) Y2 C+ H% r* w, i
}8 @4 m/ q* F* d' I \( E
*/- N" _! s) H0 s) @! Y1 G5 C! v
//comment out all methods, N=100000 202s
}; P- t, ^6 j" G! o$ b" b- g
delete []b1;
delete []b2;
机器猫 发表于 2022-9-27 00:15. W4 {5 K0 z& N G, R. U& J8 u0 p
瞎猜一下啊。把第一个的那个j定义成register变量会不会有不同？; H- f4 n5 I+ P6 W7 Z
# c R* A- I& e: i. H5 h
你第二个试验里面的j在循环里面又重新定义 ...
雷达 发表于 2022-9-27 01:16
内循环里面的 j 实际是 i, 为了规避爱坛显示的冲突帖子里临时改成了j, 现在是 jj 了。好累 、LOL
不和它 ...
雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB)% a- k3 _% ^# m' u; |% O
{
comp temp, xtimesy;
雷达 发表于 2022-9-26 01:30' e/ u+ C+ J0 p1 t
理了理思路，重新做了一个测试。
做了两个 vector 和 两个 float *, 都长 100000% Y y6 }8 R& W- f |4 g. X
外循环 6000，里面先做随 ...
雷声 发表于 2022-9-27 20:390 n/ J/ G H1 M( ]
这个时间是从哪里开始算的？
我怀疑这个200多秒里面有200秒花在产生随机数上了，真正计算大概只用了2秒， ...
雷达 发表于 2022-9-24 23:54
void xcorr(comp* outcomp, comp* A, int lenA, comp* B, int lenB), v) G$ n4 q( S. C0 Z4 Z
{
comp temp, xtimesy;
opensrc 发表于 2022-9-28 00:29
我有些迷糊，这样的code，难道不就应该时间差很多吗？也做了个简单的实验，你看看我做的有错吗
' V# k. `. [ z5 j; y1 R* l" `
...
。( @1 v0 r2 q; } u* S4 ^
| 欢迎光临 爱吱声 (http://aswetalk.net/bbs/) | Powered by Discuz! X3.2 |