OpenMP多核编程的配置与测试

现在的集成开发环境很好的支持了市面上的多数多核CPU，比如在intel双核CPU的系统上安装好vs2008后，

第一，在项目属性-c++-语言中选择支持openmp

第二，在环境变量中加入OMP_NUM_THREADS变量赋值2

可以了，打开vs2008

敲入下面的例子，相信我们都能看懂的：

/************************************************************************/
/* 利用Monte Carlo算法求圆周率pi                                                                     */
/************************************************************************/
#include <iostream>
#include <stdio.h>
#include <cmath>
#include <ctime>
#include <omp.h>
#include <time.h>

using namespace std;

#define pear 5000000//扔的“黄豆”数目
int isInCir(double x,double y)
{
    if ((x*x+y*y) <= 1)
    {
        return 1;//在1/4园内的话，返1
    }
    return 0;
}

int main()
{
    float kaishi;
    float jieshu;
    int i;
    int num = 0;//在园内的豆豆数目
    double x,y;
    double mianji;
    srand((unsigned)time(NULL));

    kaishi =clock();
    //singer
    int singerNum = 0;
    for (i=0;i<pear;i++)
    {
        x = rand()%RAND_MAX/float(RAND_MAX);
        y = rand()%RAND_MAX/float(RAND_MAX);

        if (isInCir(x,y) == 1)
        {
            singerNum++;
        }
    }
    jieshu = clock();
    cout<<"单核处理时间为："<<jieshu-kaishi<<endl;

    //双核2线程
    kaishi = clock();
    omp_set_num_threads(2);
#pragma omp parallel for reduction(+ : num)
    for (i=0;i<pear;i++)
    {
        x = rand()%RAND_MAX/float(RAND_MAX);
        y = rand()%RAND_MAX/float(RAND_MAX);

        if (isInCir(x,y) == 1)
        {
            num++;
        }
    }
    jieshu = clock();
    cout<<"双核两个线程处理时间为："<<jieshu-kaishi<<endl;

    //双核3-200线程
    int moreNum =0;
    for(int j = 3;j<200;j=j+20){
        moreNum =0;
    kaishi = clock();
    omp_set_num_threads(j);
#pragma omp parallel for reduction(+ : moreNum)
    for (i=0;i<pear;i++)
    {
        x = rand()%RAND_MAX/float(RAND_MAX);
        y = rand()%RAND_MAX/float(RAND_MAX);

        if (isInCir(x,y) == 1)
        {
            moreNum++;
        }
    }
    jieshu = clock();
    cout<<"双核"<<j<<"个线程处理时间为："<<jieshu-kaishi<<endl;
    }
    cout<<"豆豆落在园内的个数："<<num<<endl;
    cout<<"最大随机数："<<RAND_MAX<<endl;
    mianji = (num*1.0)/pear;
    //printf("Monte Carlo算法模拟出的半径为一的1/4圆的面积为%f\n",mianji);
    cout<<"Monte Carlo算法模拟出的半径为一的1/4圆的面积为"<<mianji<<endl;
    cout<<"计算出的pi值为："<<4*mianji<<endl;
    return 0;
}