python - python array MemoryError：无法分配形状为（250000、250000）的数组

Question

出于分析反演阵列研究的目的，我需要构建一个形状为（250000、250000）的对角矩阵。

nx=500
ny=500
D=np.diag(np.full(nx*ny,1.,dtype=np.float32))
Dinv=inv(D)

不幸的是，这个数组似乎太大了（？？），我收到了这个错误信息：

MemoryError: Unable to allocate array with shape (250000, 250000) and data type float32

我该如何解决这个错误？

标签： pythonarrays