一个OpenMP的学习程序

贴一个学习OpenMP过程中的程序

  1#include <stdio.h>
  2#include <tchar.h>
  3#include <omp.h>
  4int _tmain(int argc, _TCHAR* argv[])
  5{    
  6    Timer tm;
  7//    omp_set_num_threads(2),设定活动线程数为2
  8    omp_set_num_threads(3);
  9//    omp_get_num_procs()
 10    int procs=omp_get_num_procs();
 11    printf("There are %d processor(s) available!\n",procs);
 12//  omp_get_thread_num(),返回当前线程号;
 13    int Threadid=omp_get_thread_num();
 14    printf("当前线程号为:%d\n",Threadid);
 15//    omp_get_num_threads(),返回当前线程总数;
 16    int Threadnum = omp_get_num_threads();
 17    printf("当前活动线程数为:%d\n",Threadnum);
 18
 19// 并行化编译指导语句
 20    printf("下面开始测试并行化编译指导语句!!\n\n");
 21// 1.parallel
 22    printf("------------------华丽无比的分隔线--------------------\n");
 23    printf("#pragma omp parallel\n\n");
 24
 25#pragma omp parallel
 26
 27    printf("This is Thread %d say hello to you!\n",omp_get_thread_num());
 28
 29// 2.parallel for
 30    printf("------------------华丽无比的分隔线--------------------\n");
 31    printf("#pragma omp parallel for\n\n");
 32
 33#pragma omp parallel for
 34    for (int i=0;i<4;i++)
 35    {
 36        printf("Thread %d is active when i = %d\n",omp_get_thread_num(),i);
 37    }
 38    int i,j=0;
 39// 3.parallel for private
 40    printf("------------------华丽无比的分隔线--------------------\n");
 41    printf("#pragma omp parallel for private(j)\n\n");
 42
 43    printf("使用private子句\n");
 44#pragma omp parallel for private(j)
 45    for (i=0;i<2;i++)
 46        for(j=5;j<7;j++)
 47            printf("Thread %d is active when i = %d,and j= %d\n",omp_get_thread_num(),i,j);
 48    printf("不使用private子句\n");
 49#pragma omp parallel for 
 50    for (i=0;i<2;i++)
 51//#pragma omp parallel for
 52        for(j=5;j<7;j++)
 53            printf("Thread %d is active when i = %d,and j= %d\n",omp_get_thread_num(),i,j);
 54// 4.parallel for firstprivate() and lastprivate()
 55    printf("------------------华丽无比的分隔线--------------------\n");
 56    printf("#pragma omp parallel for firstprivate and lastprivate\n\n");
 57
 58    int ai=10;
 59    int bi=0;
 60    printf("Before \"parallel for \", ai = 10 ,and bi = 0 !\n\n");
 61#pragma omp parallel for firstprivate(ai) lastprivate(bi) 
 62        for (int i=1;i<4;i++)
 63        {        
 64            bi=ai;
 65            ai=0;
 66        }
 67    printf("After \"parallel for \" bi equal to %d !\n\n",bi);
 68    printf("After \"parallel for \" ai still equal to %d !\n\n",ai);
 69// 5.critical 临界区
 70    printf("------------------华丽无比的分隔线--------------------\n");
 71    printf("#pragma omp critical \n\n");
 72
 73    int n=100000;
 74    int ver = 0;
 75    tm.start();
 76#pragma omp parallel for 
 77    for (int i=1;i<=n;i++)
 78    {
 79        //for (j=0;j<1000000;j+=2){;}
 80#pragma omp critical
 81        ver+=i;
 82        //printf("Thread %d is active when i = %d\n",omp_get_thread_num(),i);
 83    }
 84    printf("ver = %d\n",ver);
 85    tm.stop();
 86    printf("%f\n",tm.elapsedSeconds());
 87
 88// 6.reduction(<op>:<variable>)
 89    printf("------------------华丽无比的分隔线--------------------\n");
 90    printf("#pragma omp reduction \n\n");
 91    
 92    int ver2 = 0;
 93    tm.start();
 94#pragma omp parallel for reduction(+:ver2)
 95    for(int i=0;i<=n;i++)
 96    {
 97        //for (j=0;j<1000000;j+=2){;}
 98        ver2+=i;
 99    }

    printf("ver2 = %d\n",ver2);
    tm.stop();
    printf("%f\n",tm.elapsedSeconds());

// 7.#pragma omp for
    printf("------------------华丽无比的分隔线--------------------\n");
    printf("#pragma omp for \n\n");
//将内层for循环并行化
#pragma omp parallel 
    for (int i=0;i<2;i++)
    {
#pragma omp for
        for (int j=0;j<2;j++)
        {
            printf("Thread %d is active when i=%d,and j=%d\n",omp_get_thread_num(),i,j);
        }
    }
// 8.threadprivate and copyin
    printf("------------------华丽无比的分隔线--------------------\n");
    printf("#pragma omp threadprivate and copyin \n\n");

    int global=10;
//#pragma omp threadprivate(global)
//#pragma omp parallel
    for (int i=0;i<100;i++)
    {
        global++;
    }

    printf("Global is equal to %d\n",global);


// 9.功能并行#prgama omp parallel sections 
    printf("------------------华丽无比的分隔线--------------------\n");
    printf("#pragma omp parallel sections \n\n");

#pragma omp parallel sections
    {
#pragma omp section
        printf("Section 1 's thread is NO.%d\n",omp_get_thread_num());
#pragma omp section
        printf("Section 2 's thread is NO.%d\n",omp_get_thread_num());
#pragma omp section
        printf("Section 3 's thread is NO.%d\n",omp_get_thread_num());
    }

// 10.同步编译指导语句 
//    1).同步屏障
#pragma omp parallel
    {
        int initialization=0;
#pragma omp barrier
        int process=1;
    }

//    2.互斥锁机制(库函数)
#ifdef LOCKLOCK
#define LOCKLOCK
    omp_lock_t lock;
//普通锁
    omp_init_lock(&lock);
    omp_destroy_lock(&lock);
    omp_set_lock(&lock);
    omp_unset_lock(&lock);
    omp_test_lock(&lock);
//嵌套锁
    omp_init_nest_lock(&lock);
    omp_destroy_nest_lock(&lock);
    omp_set_nest_lock(&lock);
    omp_unset_nest_lock(&lock);
    omp_test_nest_lock(&lock);
#endif
//其他的OpenMP编译指导语句请参考OpenMP手册;



    return 0;
}