poster
2019-12-10, 20:30
我刚刚开始学习如何使用OpenMP (http://en.wikipedia.org/wiki/OpenMP) 。我试图弄清楚为什么以下代码不能与Visual Studio 2008并行运行。它可以编译并正常运行。但是，它在我的四核计算机上仅使用一个核。这是我试图移植到MATLAB (http://en.wikipedia.org/wiki/MATLAB) mex函数的代码的一部分。任何指针表示赞赏。
#pragma omp parallel for default(shared) private(dz, t, v, ts_count) reduction(+: sum_v) for(t = 0; trandn(); v += mrdt* (tv - v) + vv_v_sqrt_dt * dz + vv_vv_v_dt*(dz*dz - 1.); sum_v += v; if(t == ts_count-1) { int_v->at_w(k++) = sum_v/(double)(t+1); ts_count += ts; } }
回答：
使用上一次迭代的v值计算v变量
for(t = 0; t
#pragma omp parallel for default(shared) private(dz, t, v, ts_count) reduction(+: sum_v) for(t = 0; trandn(); v += mrdt* (tv - v) + vv_v_sqrt_dt * dz + vv_vv_v_dt*(dz*dz - 1.); sum_v += v; if(t == ts_count-1) { int_v->at_w(k++) = sum_v/(double)(t+1); ts_count += ts; } }
回答：
使用上一次迭代的v值计算v变量
for(t = 0; t