在 Fipy 中求解的线性系统

Question

我是否通过缓存 Matrix 和 RHS 并运行求解/扫描 1 次迭代来获得线性系统？另外，是否可以手动设置正在使用的矩阵？

谢谢

Answer 1

如果您只需要矩阵和 RHS，您可以在方程式上调用 Term._prepareLinearSystem()。

虽然可以使用SparseMatrix.put() 和SparseMatrix.addAt() 构建您自己的矩阵，但FiPy 并不是真的打算以这种方式工作。如果您只想构建和求解自己的线性系统，最好直接使用底层求解器库。