我想在分布式环境中进行Cholesky分解.为此目的,我使用pdpotrf().但是,我正在努力理解函数所需的参数,并且它们没有提供关于如何使用它的C示例(并且一个例子真的很棒).
假设我有一个我要分解的NxX矩阵.那么,参数应该具有什么值? uplo,a和info在我的脑海中有明确的定义.其余的怎么样?
我应该说n应该等于N.然而,desca,ia和ja让我很困惑. desca是全球性和本地化的,我无法理解.谁能解释一下?也许有一个方法示例(不一定是代码,因为我想了解,而不仅仅是复制粘贴).
我不知道任何Fortran.
解决方法:
从这answer,我们有:
IA
andJA
: 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
andJA
both have to be1
!
同样,从这个答案中,有一个描述desca的link: