在 scipy 中使用共轭梯度法求解线性方程组

问题描述

我想使用 scipy.sparse.linalg.cg 求解格式为 Ax=b 的方程组。 我使用以下代码行

from scipy.sparse.linalg import cg
x=cg(arg_stiffness,arg_force,1e-10,300 )

arg_stifness 和 arg_force 是 csr_matrix 格式的形状为 (593457,593457) 和 (593457,1) 的稀疏数组，容差为 1e-10，最大迭代次数为 300。

当我运行此代码时，出现以下错误：

ValueError: A and b have incompatible dimensions

有人可以帮我解决这个问题吗？

非常感谢！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）