对pdpotrf参数感到困惑

Question

我想在分布式环境中进行Cholesky分解。为此，我使用pdpotrf()。但是，我正在努力理解函数所需的参数，并且它们没有提供关于如何使用它的C示例（并且一个例子真的很棒）。

假设我有一个想要分解的NxX矩阵。那么，参数应该具有什么值？ uplo，a和info在我的脑海中有明确的定义。其余的怎么样？

我会说，

n应该等于N.但是，desca，ia和ja会让我感到困惑。此外，desca是全球性和本地性的，我无法理解。谁能解释一下？也许用一个方法示例（不一定是代码，因为我想理解，而不仅仅是复制粘贴）。

_{我不知道任何Fortran。}

Answer 1

从answer开始，我们有：

  

IA和JA：这些参数用于在更大的矩阵中提供全局矩阵的起始行和列。只有当你有一个大矩阵并且只想要一个子矩阵的Cholesky分解时，它们才是相关的。在您的情况下，IA和JA都必须为1！

同样，从这个答案中，有一个link，它描述了desca：