计算科学 - 有没有办法用 PETSc 检查稀疏矩阵的图形？ - 吾爱随笔录

计算科学矩阵宠物稀疏矩阵图论

2021-12-21 22:03:49

我目前正在尝试在 PETSc 框架内编写CA-CG 方法。此过程中的一个强制性步骤是为通用稀疏矩阵实施“矩阵幂核”算法。

这需要从每个进程“拥有”的本地行的索引集开始，计算探索与稀疏矩阵相关联的图以进行 s 步获得的索引集。

是否有使用 PETSc API 检索此类信息的干净方法？

1个回答

使用MatIncreaseOverlap()just like PCASM(additive Schwarz preconditioner) 获得重叠区域，然后MatGetSubMatrices()(also like PCASM) 拉出矩阵的重叠部分。

有关实施的进一步技术讨论，我们鼓励您订阅petsc-dev@mcs.anl.gov或仅发送电子邮件petsc-maint@mcs.anl.gov。

其它你可能感兴趣的问题