Open MP计算π的值

<pre name="code" class="cpp">#include<stdio.h>

#include<time.h>

#include<omp.h>

 

static long num_steps = 100000;

double step,pi;

void main()

{

    clock_t t1=clock();

    int i;

    double x,sum = 0.0;

    step = 1.0/(double)num_steps;

#pragma omp parallel for reduction(+:sum)

    for(i = 0;i<num_steps;i++)

    {

        x = (i+0.5)*step;

        sum = sum + 4.0/(1.0+x*x);

        printf("%d",omp_get_thread_num());

    }

    pi = step *sum;

    clock_t t2 =clock();

    printf("pi=%f,用时:%d\n",pi,t2-t1);

}
上述代码使用了规约rduction(op:list);这样每个线程会有一个私有的list,当所有线程的list计算完成之后,每个线程的私有list值将被op规约为一个list值。如果不进行归约的话可能会导致每次的运行结果都不一样
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值