Pytrilinos Amesos 求解器参数 - UseTranspose

编程问答 2022-05-31

问题描述

我正在尝试使用 Pytrilinos Amesos 求解器（特别是 Lapack 或 Klu）求解 30000x1000 稀疏矩阵。我知道我需要使用转置来解决这个问题，并且我已将 Amesos UseTranspose 参数设置为“True”：

problem = Epetra.LinearProblem(Epmat,LHS,b)

solver = Amesos.Lapack(problem)

# Set the solver parameters
amesosList = {"PrintStatus": True,"PrintTiming": True,"UseTranspose": True
              }
solver.SetParameters(amesosList)

solver.SymbolicFactorization()
solver.NumericFactorization()
solver.solve()

但是它似乎不起作用，相反，我需要在解决它们之前手动将矩阵 A 和 b 乘以转置等等......我是否设置了错误的参数还是我遗漏了什么？

我是稀疏矩阵和这​​个 Trilinos 包的初学者，没有太多关于这方面的文档，所以我有点迷茫 - 任何帮助将不胜感激。

解决方法

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

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

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

lapackmatrixmatrixpytrilinossparse-matrixtrilinos