我有一个大小为12 GB的numpy-memmap矩阵S
。我正在尝试对每一行进行argsort。为此,我定义了另一个内存映射数组first_k
以保存结果。问题是发生内存错误。
代码如下:
first_k = np.memmap('first_k', dtype='float32', mode='w+', shape=S.shape)
first_k[:] = np.memmap.argsort(S, axis=1)
任何可能的解决方案? 我正在考虑将其切片处理...
预先感谢
答案 0 :(得分:0)
最后,我在切片中处理memmap。这是代码。
N = S.shape[0]
first_k = np.memmap('first_k', dtype='float32', mode='w+', shape=(N, N))
del first_k
step = 1000
for row in np.arange(0, N, step):
size = min(step, N-row)
first_k = np.memmap('first_k',
dtype='float32',
mode='r+',
shape=(size, N),
offset=4*N*row)
first_k[:] = np.memmap.argsort(S[row:row+size], axis=1)
del first_k