1

我有一个非常(非常非常)大的二维数组 - 大约一千列,但有几百万行(足以让它不适合我的 32GB 机器的内存)。我想计算每一千列的方差。一个有帮助的关键事实:我的数据是 8 位无符号整数。

这就是我打算如何处理这个问题。我将首先构造一个名为counts的新二维数组,形状为 (1000, 256),其思想是counts[i,:] == np.bincount(bigarray[:,i]). 一旦我有了这个数组,计算方差就很简单了。

麻烦的是,我不确定如何有效地计算它(这个计算必须实时运行,我希望带宽受限于我的 SSD 可以多快返回数据)。这是一些有效的方法,但速度非常慢:

counts = np.array((1000,256))
for row in iterator_over_bigaray_rows():
    for i,val in enumerate(row):
        counts[i,val] += 1

有没有办法写这个运行得更快?像这样的东西:

counts = np.array((1000,256))
for row in iterator_over_bigaray_rows():
    counts[i,:] = // magic np one-liner to do what I want
4

1 回答 1

1

我认为这就是你想要的:

counts[np.arange(1000), row] += 1

但是,如果您的数组有数百万行,您仍然需要遍历数百万行。以下技巧在我的系统上提供了接近 5 倍的加速:

chunk = np.random.randint(256, size=(1000, 1000))

def count_chunk(chunk):
    rows, cols = chunk.shape
    col_idx = np.arange(cols) * 256
    counts = np.bincount((col_idx[None, :] + chunk).ravel(),
                         minlength=256*cols)
    return counts.reshape(-1, 256)

def count_chunk_by_rows(chunk):
    counts = np.zeros(chunk.shape[1:]+(256,), dtype=np.int)
    indices = np.arange(chunk.shape[-1])
    for row in chunk:
        counts[indices, row] += 1
    return counts

现在:

In [2]: c = count_chunk_by_rows(chunk)

In [3]: d = count_chunk(chunk)

In [4]: np.all(c == d)
Out[4]: True

In [5]: %timeit count_chunk_by_rows(chunk)
10 loops, best of 3: 80.5 ms per loop

In [6]: %timeit count_chunk(chunk)
100 loops, best of 3: 13.8 ms per loop
于 2013-03-27T20:54:23.433 回答