[英]MPI Cartesian Topology: MPI_Neighbor_alltoall wrong data received
[英]Accessing data in neighbor cells
我有一个大小为mxn
的矩阵M
,它保存为长度为m * n
的一维数组N
该数组的每个单元格都包含一些整数变量,这些变量是数据点的ID。 每个单元中的整数变量的数量随时间变化。
N[0] = {1,4,5,7}
N[1] = {2,9,3,1,7,4}
N[2] = {7,1,3,9,8}
N[3] = {6,4,2}
...
我使用返回的索引函数访问这些元素
idx = x + y * n
给定一些索引idx
我想使用邻居单元和具有索引idx
的中心单元的所有整数变量来访问大小为s
的数据点D
的数组。 大小s
可能非常大。
我的观点很清楚:不是在所有数据点上都这样循环
for(int i=0; i<s; i++)
// Do something with D[i]
我想要类似的东西(但更紧凑)
// Access central cell
idx = x + y*n;
num_Elements = Number_of_Elements_Cell(x,y);
for(int i=0; i<num_Elements; i++)
// Do something with D[N[idx][i]]
// Access right cell
idx = (x+1) + y*n;
num_Elements = Number_of_Elements_Cell(x+1,y);
for(int i=0; i<num_Elements; i++)
// Do something with D[N[idx][i]]
// Access left cell
idx = (x-1) + y*n;
num_Elements = Number_of_Elements_Cell(x-1,y);
for(int i=0; i<num_Elements; i++)
// Do something with D[N[idx][i]]
等等。 对于所有单元格,我必须执行9次。
我的问题:给定结构N
是否有更好的方法呢?
我不确定我是否很好地理解了您的问题....但是您可以尝试:
for (int i=-1 ; i <= 1 ; i++){
for (int j = -1 ; j <=1 ; j++){
idx = (x+i) + (y+j)*n;
// Check if idx is not out of bounds
num_Elements = Number_of_Elements_Cell(x+i,y+j);
for(int k=0; k<num_Elements; k++)
// Do something with D[N[idx][k]]
}
}
请注意,您的索引很可能超出了这种方法的范围,因此您必须实施测试以防止这种情况的发生。
这就是使用double for
循环简单地对单元及其8个邻居进行迭代的方法。 如果不是您所期望的,请告诉我,我将进行编辑/删除。
我不确定,但也许您正在寻找这样的东西:
var distinctDataPoints = new List<int>();
for(int z = x - 1, z <= x + 1, z++)
{
if(z < 0 || z > m)
continue;
for(int t = y-1, t <= y + 1, t++)
{
if(t < 0 || t > n)
continue;
idx = z + t * n;
for(int i = 0; i < num_Elements; i++)
{
if(!distinctDataPoints.Contains(N[idx][i])
distinctDataPoints.Add(N[idx][i])
}
}
}
for(int dpIdx = 0; dpIdx < distinctDataPoints.Count; dpIdx++)
{
//Do something with D[dpIdx]
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.