c++ - 与 pdpotrf 参数混淆

Question

我想在分布式环境中进行 Cholesky 分解。为此，我使用 pdpotrf() .但是，我很难理解该函数所需的参数，并且他们没有提供有关如何使用它的 C 示例(如果有一个示例就太好了)。

假设我有一个要分解的 NxX 矩阵。那么，参数应该取什么值呢？ uplo, a 和 info 在我的脑海里定义的很好。剩下的呢？

n 应该等于 N 我会说。然而，desca、ia 和 ja 是让我感到困惑的。而且，desca 是全局的和本地的，这是我无法理解的。有人可以解释吗？也许有一个方法示例(不一定是代码，因为我想理解，而不仅仅是复制粘贴)。

_{我不懂任何 Fortran。}

Answer 1

从这里answer ，我们有:

IA and JA: those parameters are meant to provide the starting row and column of your global matrix inside a larger matrix. They are only relevant if you have a big matrix and only want the Cholesky decomp of a submatrix. In your case, IA and JA both have to be 1!

同样，根据这个答案，存在一个 link ，它描述了 desca: