你确定？1万*1万方阵的平方？

import numpy as np

m = 10000
adjacent_matrix = np.random.random((m,m))
m_squared = adjacent_matrix.dot(adjacent_matrix)

【 在 cafitren 的大作中提到: 】
: 你得换电脑了，我笔记本8G内存1万x1万都是几秒钟。
:
--
FROM 219.237.112.*

对于二维矩阵，用dot还是用@或者matmul，都是一样一样的……
【 在 cafitren 的大作中提到: 】
: 囧，你这从哪学的矩阵相乘用dot的。一般用@或者matmul吧。
:
--
FROM 39.144.39.*

感谢

之前因为还要干Django等别的事情，没有使用anaconda，也就没有blas。才知道anaconda的numpy默认安装了blas的。

配置了mkl或者blas以后就是这个速度了

【 在 bsxfun 的大作中提到: 】
: 不过你电脑确实太差了，或者是你没有装blas
: 我在我的小米手机（Mix4）上运行 3000 * 3000矩阵的平方也就4分钟
: 在我的笔记本电脑上，运行只需要0.15秒
: ...................
--
修改:Blue FROM 219.237.112.*
FROM 219.237.112.*