关键部分openMP

Question

我想知道我们需要在哪里设置关键部分？

如果有多个线程具有一个共享数组，并且每个线程都想写在不同的位置，即使每个线程都写到数组中的不同位置，它也都需要放在关键部分吗？

可以说我有二维数组M [3] [3]，initial_array [3]和一些双精度变量，我想计算一些东西并将其存储在M的第一列。我可以使用for循环，但是我想与openMP一起使用，所以我做了：

 omp_set_num_threads(3);
 #pragma omp parallel shared(M,variable)
 {  
     int id = omp_get_thread_num();
     double init = initial_array[id]*variable;
     M[id][0] = init;
 }   

它工作正常，但我知道它可能导致死锁或运行时间不好。 我的意思是，如果我有更多的线程甚至更大的线程，该怎么办。设置关键部分的正确方法是什么？ 我想问的另一件事是有关initial_array，是否也需要共享？

Answer 1

这是安全代码。

数组中的随机访问不会导致对数组中其他元素的竞争。 只要您继续并发地读写数组中未共享的元素，您就永远不会遇到竞争状况。

请记住，取决于元素的类型和大小，读取可以与写入竞争。 您的示例显示为double ，如果您在同一元素上同时进行读取操作和写入操作，我将很担心。 在写入过程中可能会发生上下文切换，但这取决于您的体系结构/平台。 无论如何，您没有这样做，但是值得一提。

Answer 2

由于您正在访问内存的不同部分（数组的不同索引），因此我认为并发性没有任何问题，但是如果您的内核具有专用的L1高速缓存，我看到的唯一问题是性能下降。

在这种情况下，由于高速缓存一致性会导致性能下降，其中一个更新索引，使其他索引无效，执行回写等操作。对于少量线程/核心，这不是问题，但对于在大量核心上运行的线程，它可以确保它。 由于线程所运行的数据并不是真正独立的 ，因此它们将作为高速缓存中的数据块读取（如果您正在访问M [0] [0]，则不仅M [0] [0]会被读入缓存，但M [0] [0]至M [n] [col]，其中n取决于缓存块大小。 如果块很大，则可能包含更多共享数据。