也许我的问题的解决方案非常明显。我想使用 openMP(进入英特尔至强融核(加速以下代码,但我不能。
int c=0, d=0, e=0, i;
#pragma opm parallel for private(c, d, e)
for(i=0; i < columns; i++)
{
if((left_side[rule*columns + i] > 0) || (right_side[rule*columns + i] > 0))
{
for(c=0; c < rows; c++)
{
if(left_side[i + c * columns] > 0)
{
if(flag[c] == 0)
{
r_prob[c] = c_vect[c];
flag[c] = 1;
for(d=0; d < columns; d++)
{
switch(left_side[c * columns + d])
{
case 0:
break;
case 1:
r_prob[c] *= M_in[d] * 1.0;
break;
case 2:
r_prob[c] *= (M_in[d] * (M_in[d] - 1)) * .5;
break;
default:
for(e=1; e <= left_side[c * columns + d]; e++)
r_prob[c] *= M_in[d] * 1.0 / (e * 1.0);
break;
}
}
}
}
}
}
}
//where r_prob, M_in, left_side, right_side, c_vect and flag are array that are in input.
此代码无法正常运行,r_prob中的值是错误的。如何并行化此代码?
您的代码存在错误,因为第一个 FOR 创建多个线程,并且这些线程同时写入数组 r_prob, 标志。解决方案可能是:
int c=0, d=0, e=0, i;
#pragma opm parallel for private(c, d, e)
for(i=0; i < columns; i++)
{
if((left_side[rule*columns + i] > 0) || (right_side[rule*columns + i] > 0))
{
for(c=0; c < rows; c++)
{
if(left_side[i + c * columns] > 0)
{
if(flag[c] == 0)
{
#pragma omp critical //Only one thread at a time can enter
{
r_prob[c] = c_vect[c];
flag[c] = 1;
}
for(d=0; d < columns; d++)
{
switch(left_side[c * columns + d])
{
case 0:
break;
case 1:
#pragma omp critical //Only one thread at a time can enter
{
r_prob[c] *= M_in[d] * 1.0;
}
break;
case 2:
#pragma omp critical //Only one thread at a time can enter
{
r_prob[c] *= (M_in[d] * (M_in[d] - 1)) * .5;
}
break;
default:
for(e=1; e <= left_side[c * columns + d]; e++)
#pragma omp critical //Only one thread at a time can enter
{
r_prob[c] *= M_in[d] * 1.0 / (e * 1.0);
}
break;
}
}
}
}
}
}
}//where r_prob, M_in, left_side, right_side, c_vect and flag are array that are in input.
问题是不同的线程写入和读取同一个r_prob[c]
变量。使用此算法,您只能并行化c
上的循环。即写:
for(int i=0; i < columns; i++)
{
if((left_side[rule*columns + i] > 0) || (right_side[rule*columns + i] > 0))
{
#pragma omp parallel for
for(int c=0; c < rows; c++){
//....
或者您可以使用以下语法:
#pragma omp critical
{
r_prob[c] *= ...
}
一次只允许一个线程访问数组的语法。
哪一个更快取决于数组的大小和线程数。
正如 Haatschii 所评论的那样,最好在 for 循环中声明变量。