Question

我有一个程序可以将2个大数字作为char矢量提供。现在我必须确保它使用OpenMP。我有嵌套循环的问题，如果我使用它们，结果不是我所期待的。我非常感谢你的帮助。这是代码

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<string.h>
#include<omp.h>
#include<time.h>
#define MAX 100000

char * multiply(char [],char[], int manyThreads);
int main(){   
    //omp_set_dynamic(1);
    omp_set_nested(1);
    omp_set_num_threads(8);
    char a[MAX];
    char b[MAX];
    char *c;    

    printf("First number : ");
    scanf("%s",a);
    printf("Second number : ");
    scanf("%s",b);
    printf("Result : ");
    double start = omp_get_wtime( );

    c = multiply(a,b,1);
    printf("%s\n",c);
    double end = omp_get_wtime( );
    printf("Calculation time = %.16g\n", end - start);

    return 0;
}

char * multiply(char a[],char b[], int manyThreads){
    static char result[MAX];
    char tempResult[MAX];
    char temp[MAX];
    int aLength,bLength;
    int i,j,k=0,x=0,y;
    long int r=0;
    long sum = 0;
    aLength=strlen(a)-1;
    bLength=strlen(b)-1;

    #pragma omp parallel if(manyThreads == 1)
    {
        #pragma omp for schedule(dynamic) nowait
        for(i=0;i<=aLength;i++) 
        {
            a[i] = a[i] - 48;
        }

        #pragma omp for schedule(dynamic) nowait
        for(i=0;i<=bLength;i++) 
        {
            b[i] = b[i] - 48;
        }
    }

    #pragma omp parallel if(manyThreads == 1) 
    {       
        #pragma omp for schedule(dynamic) 
        for(i=bLength;i>=0;i--) 
        {
             r=0;

             #pragma omp parallel
             {
               #pragma omp for schedule(dynamic)
               for(j=aLength;j>=0;j--) 
               {
                  temp[k++] = (b[i]*a[j] + r)%10;
                 r = (b[i]*a[j]+r)/10;
               }
             }
             temp[k++] = r;
             x++;
             #pragma omp parallel
             {
               #pragma omp for schedule(dynamic)            
               for(y = 0;y<x;y++)
               {
                  temp[k++] = 0;
               }    
             }                   
        }
    }  

    k=0;
    r=0;
    #pragma omp parallel if(manyThreads == 1)
    {
        #pragma omp for schedule(dynamic) 
        for(i=0;i<aLength+bLength+2;i++)
        {
             sum =0;
             y=0;
             #pragma omp parallel
             {
               #pragma omp for schedule(dynamic)
               for(j=1;j<=bLength+1;j++)
               {
                  if(i <= aLength+j)
                  {
                      sum = sum + temp[y+i];
                  }
                  y += j + aLength + 1;
               }
             }

             tempResult[k++] = (sum+r) %10;
             r = (sum+r)/10;            
        }
    }
    tempResult[k] = r;
    j=0;

    #pragma omp parallel if(manyThreads == 1)
    {
        #pragma omp for schedule(dynamic) 
        for(i=k-1;i>=0;i--)
        {
             result[j++]=tempResult[i] + 48;
        }
    }

    result[j]='\0';

    if(result[0]==48)
    {
        result[0]=255;
    }

    return result;
}

Answer 1

我可以确认，你有一些数据争用错误没有OpenMP：

First number : 123456
Second number : 654321
Result :  08563613376
Calculation time = 0.005543371655221563

使用OpenMP：

First number : 123456
Second number : 654321
Result :  00000000825
Calculation time = 0.007188999978097854

我并没有解决整个问题，但根据我的看法，我有一些评论首先我要说你应该在

中指定for变量

#pragma omp for schedule(dynamic) nowait
        for(int i=0;i<=aLength;i++) 
        {
            a[i] = a[i] - 48;
        }

你所做的主要错误的一个例子（每次）......当你进入一个平行区域时，你必须要处理你的变量。什么是共享的，每个线程的私有内容是什么？在您的结束循环中，您不使用nested parallelism，但在每个帖子中都执行j++，因此您应该保护j，如下所示。

#pragma omp parallel if(manyThreads == 1) shared(j) 
    {
        int pos;
#pragma omp for schedule(dynamic) 
        for(int i=k-1;i>=0;i--)
        {
#pragma omp atomic capture
            pos = j++;

            result[pos]=(tempResult[i] + 48);
        }
    }

请注意，您还忘记了OpenMP for构造具有reduction子句

因此可以使用

重写计算sum += temp[y+i]的for循环

#pragma omp for schedule(dynamic) reduction(+:sum)

您必须遵循所有算法逻辑，在每个周围和嵌套shared区域使用private和parallel子句。

在OMP C中嵌套for循环

1 个答案: