c - MPI 三角拓扑

Question

我需要对(上)三角矩阵进行一些计算，我正在考虑为此使用 MPI。似乎为它定义一个合适的拓扑会很方便。笛卡尔拓扑朝着那个方向发展，但还不是我所需要的。我需要一个拓扑结构，它可以给出进程在严格上三角矩阵上的等级，即严格在对角线上方。

我需要的三角形拓扑类型如下所示(比如 3x3 笛卡尔网格):

i       j      rank     
0       0      NULL       
0       1       1 
0       2       2
1       0      NULL       
1       1      NULL
1       2       3
2       0      NULL       
2       1      NULL    
2       2      NULL

感谢您的帮助!

Answer 1

MPI 不提供您要求的拓扑类型。它提供笛卡尔拓扑或一般图形拓扑。对于一般的图拓扑，每个等级都有一个相邻等级的列表，但是很难将半个网格映射到它上面。

您可能想要编写自己的一组例程来管理此类拓扑。您实际上需要重新实现的是将坐标映射到等级(反之亦然)的函数，以及查找给定等级的邻居的函数。您可以预先计算映射并将其作为属性缓存在相应的通信器中，即使用 MPI_COMM_GET_ATTR 和 MPI_COMM_SET_ATTR 调用。