c - 处理大型矩阵(即 10000x10000 及更大)的最快方法？

Question

关闭。这个问题需要更多focused .它目前不接受答案。












想改善这个问题吗？更新问题，使其仅关注一个问题 editing this post .

7年前关闭。



Improve this question




我目前正在使用 Cython，但我的代码涉及在循环中多次取每行/列组合的范数，因此即使是 1000x1000 矩阵也非常慢。我正在考虑切换到 C+ LAPACK .. 处理巨大矩阵的标准是什么？

看起来最昂贵的步骤包括取所有行/列组合的点积并遍历矩阵的上三角形。

矩阵不是稀疏的。

Answer 1

试试 Numpy，1000x1000 听起来有点小。
它是用 C 实现的，可以利用线性代数库，因此速度很快。阅读 this .
如果条目不多，请使用稀疏矩阵。

import numpy as np
A = np.random.rand(1000,1000)
np.linalg.norm(A, axis=1)    # takes me under 5 ms

如果您需要更具体的帮助，您应该提供更多详细信息。