c++ - 对 pdpotrf 参数感到困惑

Question

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

假设我有一个想要分解的 NxX 矩阵。那么，参数应该有什么值呢？uplo，a并且info在我的脑海中定义得很好。其余的呢？

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

_{我不知道任何Fortran。}

Answer 1

从这个答案，我们有：

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

同样，从这个答案中，存在一个链接，描述desca：