我不确定为什么比较切片的 numpy 阵列使用,
比][
. 例如:
start = time.time()
a = np.zeros((100,100))
for _ in range(1000000):
a[1:99][1:99] == 1
print(time.time() - start)
start = time.time()
a = np.zeros((100,100))
for _ in range(1000000):
a[1:99, 1:99] == 1
print(time.time() - start)
3.2756259441375732
11.044903039932251
差了 3 倍多。使用 的时间测量大致相同timeit
。
我正在研究递回算法(我打算这样做),这些问题使我的程序运行速度变慢,从大约 1 秒增加到 10 秒。我只想知道他们背后的原因。可能这是一个错误。我正在使用 Python 3.9.9。谢谢。
uj5u.com热心网友回复:
第一个与a[2:99]==1
. 一个 (98,100) 切片后跟一个 (97,100),然后是 == 测验。
In [177]: timeit (a[1:99][1:99]==1)
8.51 μs ± 16.3 ns per loop (mean ± std. dev. of 7 runs, 100000 loops each)
In [178]: timeit (a[1:99][1:99])
383 ns ± 5.73 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)
In [179]: timeit (a[1:99])
208 ns ± 10.4 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)
大部分时间是测验,而不是切片。
In [180]: a[1:99,1:99].shape
Out[180]: (98, 98)
In [181]: timeit a[1:99,1:99]==1
32.2 μs ± 12.9 ns per loop (mean ± std. dev. of 7 runs, 10000 loops each)
In [182]: timeit a[1:99,1:99]
301 ns ± 3.61 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)
同样,切片是计时的一小部分,但==
测验速度要慢得多。在第一种情况下,我们选择了行的子集,因此测验在资料缓冲区的连续块上进行。在第二个中,我们选择行和列的子集。通过资料缓冲区的迭代更为复杂。
我们可以通过测验一段列与一段行来简化比较:
In [183]: timeit a[:,2:99]==1
32.3 μs ± 13.8 ns per loop (mean ± std. dev. of 7 runs, 10000 loops each)
In [184]: timeit a[2:99,:]==1
8.58 μs ± 10.2 ns per loop (mean ± std. dev. of 7 runs, 100000 loops each)
作为进一步的测验,使用“F”顺序创建一个新阵列。现在“行”是慢切片
In [189]: b = np.array(a, order='F')
In [190]: timeit b[:,2:99]==1
8.83 μs ± 20.6 ns per loop (mean ± std. dev. of 7 runs, 100000 loops each)
In [191]: timeit b[2:99,:]==1
32.8 μs ± 31.2 ns per loop (mean ± std. dev. of 7 runs, 10000 loops each)
===
但是你为什么要比较这两个切片,一个是 (97,100) 阵列,另一个是 (98,98)。他们正在挑选a
.
我想知道您是否真的打算测验连续的行、列切片,而不是两行切片。
In [193]: timeit (a[1:99][:,1:99]==1)
32.6 μs ± 92.4 ns per loop (mean ± std. dev. of 7 runs, 10000 loops each)
仅比较切片,我们发现顺序切片更慢 - 只是一点点。
In [194]: timeit (a[1:99][:,1:99])
472 ns ± 3.76 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)
In [195]: timeit (a[1:99,1:99])
306 ns ± 3.19 ns per loop (mean ± std. dev. of 7 runs, 1000000 loops each)
0 评论