用python求解大量方程

我正在写一个有限元程序，想用scipy寻根器解一个非线性方程组。例如，公式是Ax-B=0，A是（2000，2000）数组。因此，寻根器需要花费很长的时间来计算雅可比矩阵，然后求解方程。如何加快解决过程

欢迎来到Stackoverflow。你试过什么？如果你坚持提问，你会得到更好的结果。发布你的代码可以让其他人更容易地提出改进方法。这个问题太开放了，无法得到有用的答案。如何加速你的问题完全取决于你的偏微分方程系统，因为最优方法将取决于矩阵A的稀疏性和条件。也就是说，一些起点是牛顿-克雷洛夫解算器（）和gmres（）。