我正在尝试使用 MATMPIDENSE 矩阵来解决 Ax=b 类型的系统，但我遇到了一些问题。KSP 文档说，密集矩阵的 KSPSolve 需要设置“gmres”求解器和“lu”预处理器，问题是“lu”预处理器不能并行工作，相应的“asm”不能与稠密矩阵。我该如何解决这个“循环”？

此外，我需要对这些密集矩阵进行一些操作，例如调用一个似乎不适用于 MATMPIDENSE 的 MatMatMult。是这样吗？我能做些什么 ？

我什至尝试计算我的矩阵的集合，以允许这些操作，但如果在调用 SetMatType(A,MATSEQDENSE) 和 MatGetRow(A,...) 时使用超过 1 个处理器运行，我会收到错误消息。

提前谢谢你，伊丽莎