标签：c intel-mkl linear-algebra mpi scalapack

我想在分布式环境中进行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!

同样,从这个答案中,有一个描述desca的link：

标签：c,intel-mkl,linear-algebra,mpi,scalapack
来源： https://codeday.me/bug/20191002/1844724.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。