我有一个函数void dynamics (A a, std::vector<double> &, std::vector<double> &, std::vector<double> )，我从openmp创建的线程中调用它。该函数的输入是每个线程（在并行块中创建）的私有输入。

#include <iostream>
#include <vector>
#include <chrono>
using namespace std;

class A {
 some code
};

int main(void)
{
vector<double> a (12,0.0);
vector<double> b (12,0.0);

#pragma omp parallel for shared(a,b)
for(int id = 0; id < 6; id++) {

vector<double> a_private (2,0.0);
vector<double> b_private (2,0.0);
vector<double> c_private (2,(double)id);
A d;

start_time for each thread - chrono
dynamics(d,a_private,b_private,c_private);
end_time for each thread - chrono
calculate_time for each thread

# pragma omp critical
{
for(int i = 0; i < 2; i++) a[i+(2*id)] = a_private[i];
for(int i = 0; i < 2; i++) b[i+(2*id)] = b_private[i];
}

}

print(a);
print(b);

return 0;

}

这里，为了避免争用条件，我已经将a_private和B_private的赋值放到临界区内的a和b中。
当我为每个线程计算上述代码的时间时，它比我将动态函数放在临界区中的时间要多。

# pragma omp critical
{
start_time for each thread - chrono
dynamics(d,a_private,b_private,c_private);
end_time for each thread - chrono
calculate_time for each thread
for(int i = 0; i < 2; i++) a[i+(2*id)] = a_private[i];
for(int i = 0; i < 2; i++) b[i+(2*id)] = b_private[i];
}

最后的输出（a和B）在这两种情况下是相同的（多次运行代码给予相同的结果）。因此，我相信动态是线程安全的（难道它不是线程安全的吗？）
动态的输入是在并行区域中创建的，因此，它们应该是每个线程的私有的（是吗？）。
为什么线程在一起工作时计算动态的速度比一个接一个工作时（在临界段内）慢？
我相信创建和管理线程的开销不会是一个问题，因为我比较了总是创建线程的时间（在我的上述两种情况下）。
并行化动态后的总时间低于串行版本（实现了加速），但为什么线程所用时间会有显著差异（在临界范围内，而不是：计算线程时间）。
我能想到的解释是，运行动态会产生竞态条件，即使它的输入和输出对每个线程都是私有的。（会是这样吗？）
此外，我没有使用omp get num threads和omp get thread num。
这里有什么问题？

When running dynamics in parallel
ID = 3, Dynamics Time = 410233
ID = 2, Dynamics Time = 447835
ID = 5, Dynamics Time = 532967
ID = 1, Dynamics Time = 545017
ID = 4, Dynamics Time = 576783
ID = 0, Dynamics Time = 624855

When running dynamics in critical section
ID = 0, Dynamics Time = 331579
ID = 2, Dynamics Time = 303294
ID = 5, Dynamics Time = 307622
ID = 1, Dynamics Time = 340489
ID = 3, Dynamics Time = 303066
ID = 4, Dynamics Time = 293090

（将不能提供最小的动力学再现，因为它是我的教授的专有）

• 谢谢-谢谢