与 pdpotrf 参数混淆
Confused with pdpotrf arguments
我想在分布式环境中进行 Cholesky 分解。为此，我使用 pdpotrf()。但是，我很难理解该函数所需的参数，并且他们没有提供有关如何使用它的 C 示例（如果有一个示例真的很棒）。
假设我有一个要分解的 NxX 矩阵。那么，参数应该取什么值呢？ uplo、a 和 info 在我的脑海中定义得很好。剩下的呢？
n 应该等于 N 我会说。但是，desca、ia 和 ja 让我感到困惑。而且， desca 是全局和本地的，我无法理解。有人可以解释吗？也许有一个方法示例（不一定是代码，因为我想理解，而不仅仅是复制粘贴）。
我不知道任何 Fortran。
由此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:
我想在分布式环境中进行 Cholesky 分解。为此，我使用 pdpotrf()。但是，我很难理解该函数所需的参数，并且他们没有提供有关如何使用它的 C 示例（如果有一个示例真的很棒）。
假设我有一个要分解的 NxX 矩阵。那么，参数应该取什么值呢？ uplo、a 和 info 在我的脑海中定义得很好。剩下的呢？
n 应该等于 N 我会说。但是，desca、ia 和 ja 让我感到困惑。而且， desca 是全局和本地的，我无法理解。有人可以解释吗？也许有一个方法示例（不一定是代码，因为我想理解，而不仅仅是复制粘贴）。
我不知道任何 Fortran。
由此answer，我们有：
IAandJA: 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,IAandJAboth have to be1!
同样，根据这个答案，存在一个 link，描述 desca: