在openMP中使用锁转储核心

Question

我想并行化函数S并锁定每个节点，但是我一直得到核心转储。我试图在图的每个节点上使用锁。如果我在我的节点上使用一个锁，它就能工作。

for (l = 0; l < n; l++)
omp_init_lock(&(lock[l]));

#pragma omp parallel for num_threads(16)default(none)  private(v)shared(n,Xof,lock)

for(v = 0; v < n; v++) {
    omp_set_lock(&(lock[v]));
    if(Xof[v] == NYC) 
    {
        S(v);   
        }
    omp_unset_lock(&(lock[v]));
}

Answer 1

似乎最有可能的原因是，在S中存在一个数据竞争(可能带有某种嵌入式状态)。

测试这一点的一种方法是使用n个锁，但只使用一个线程。如果这不会导致核心转储，那么S很可能有一个数据竞争。

如果它确实导致了核心转储，那么您的初始假设是由锁引起的，这可能是正确的。从代码可见，我没有看到任何明显的错误。我可能会删除共享子句，因为它是多余的，但我不认为这会产生影响。