openmp非迭代结构并行

程序函数依赖关系图

openmp非迭代结构并行_第1张图片

c代码


//gcc -lm -fopenmp -o ex4_sections ex4_sections.c && ./ex4_sections.c  
#include 
#include
#include
void  alpha(int n,int* rs)
{
    int i;
    *rs=0;
    for(i=0;i<n;i++)
       *rs+=i*2;

}
void beta(int n,int* rs)
{
    
    int i;
    *rs=0;
    for(i=0;i<n;i++)
       *rs+=i*3;

}
void  gamma1(int n,int* a,int*  b,int *rs)
{
    
    int i;
    *rs=0;
    for(i=0;i<n;i++)
       *rs+=i*(*a)/(*b);
}
void delta(int n,int* rs)
{
    
    int i;
    *rs=0;
    for(i=0;i<n;i++)
       *rs+=i*i;
}
int main()
{
    int n=10;
    int a,b,rs1,rs2;
#pragma omp parallel num_threads(4)
    {
    #pragma omp sections
        {
        #pragma omp section
        alpha(n,&a);
        #pragma omp section
        beta(n,&b);

        }
    
    #pragma omp sections
        {
        #pragma omp section
        gamma1(n,&a,&b,&rs1);
        #pragma omp section
        delta(n,&rs2);
        }   

    }
    printf("a=%d,b=%d,rs1=%d,rs2=%d,result=%d\n",a,b,rs1,rs2,pow(rs1,rs2));
    return 0;
}


结果

gcc -lm -fopenmp -o ex4_sections ex4_sections.c && ./ex4_sections
a=90,b=135,rs1=27,rs2=285,result=1081200640

你可能感兴趣的:(openmp,HPC)