1

我正在尝试使用本地指针来访问当前线程具有亲和力的内存。

不幸的是，我的本地指针似乎没有指向我认为应该指向的位置。

有人知道出了什么问题吗？

编辑：我忘了提到下面的输出是使用四个线程运行此代码生成的，即THREADS = 4.

我的代码：

#include <upc.h>
#include <stdio.h>
#include <stdlib.h>

int main(){

    shared int * T = (shared int *) upc_all_alloc(12, sizeof(int));
    if(!T)
        upc_global_exit(-1);

int i;
upc_forall(i=0; i<12; i++; &T[i]) T[i] = i;
upc_barrier;

if(MYTHREAD == 0)
    for(i=0; i<12; i++) printf("thread %d, T[%d] = %d\n", MYTHREAD, i, T[i]);
upc_barrier;

int my_start = (12/THREADS + 1)*MYTHREAD;
int my_end = (12/THREADS + 1)*(MYTHREAD+1) - 1;

int* T_local = (int*)&T[my_start];

for(i=my_start; i<=my_end; i++)
    printf("thread %d, T_local[%d] = %d, T[%d] = %d\n", MYTHREAD, 
            i-my_start, T_local[i-my_start], i, T[i]);
upc_barrier;

return 0;
}

输出（THREADS = 4）：

thread 0, T[0] = 0
thread 0, T[1] = 1
thread 0, T[2] = 2
thread 0, T[3] = 3
thread 0, T[4] = 4
thread 0, T[5] = 5
thread 0, T[6] = 6
thread 0, T[7] = 7
thread 0, T[8] = 8
thread 0, T[9] = 9
thread 0, T[10] = 10
thread 0, T[11] = 11
thread 0, T_local[0] = 0, T[0] = 0
thread 0, T_local[1] = 4, T[1] = 1
thread 0, T_local[2] = 8, T[2] = 2
thread 0, T_local[3] = 0, T[3] = 3
thread 1, T_local[0] = 4, T[4] = 4
thread 1, T_local[1] = 8, T[5] = 5
thread 1, T_local[2] = 0, T[6] = 6
thread 2, T_local[0] = 8, T[8] = 8
thread 2, T_local[1] = 0, T[9] = 9
thread 2, T_local[2] = 0, T[10] = 10
thread 2, T_local[3] = 0, T[11] = 11
thread 3, T_local[0] = 0, T[12] = 0
thread 3, T_local[1] = 0, T[13] = 0
thread 3, T_local[2] = 0, T[14] = 0
thread 3, T_local[3] = 0, T[15] = 0
thread 1, T_local[3] = 0, T[7] = 7
4

1 回答 1

1

您的数组 T 是使用循环布局分配和声明的（即 blocksize == 1）。这意味着与 MYTHREAD 有关联的第一个元素就是 T[MYTHREAD]。因此，您可能应该按如下方式初始化指向本地的指针：

int* T_local = (int*)&T[MYTHREAD];

在循环布局中，共享元素循环传递给线程，这意味着每个线程都有一个不连续的分布式数组元素块。因此，例如有 4 个线程，线程 0 将与 T[0]、T[4] 和 T[8] 具有亲和力。线程 0 上正确初始化的 T_local 指向本地的指针将访问共享数组的本地切片中的这些元素（分别为 T_local[0]、T_local[1] 和 T_local[2]）。

您对 my_start 和 my_end 的计算似乎假设了与 T 实际使用的不同（更大）的阻塞因子，这可能是您困惑的根源。

于 2014-04-13T08:01:58.040 回答