一个OpenMP的学习程序
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 }
100
101 printf("ver2 = %d\n",ver2);
102 tm.stop();
103 printf("%f\n",tm.elapsedSeconds());
104
105// 7.#pragma omp for
106 printf("------------------华丽无比的分隔线--------------------\n");
107 printf("#pragma omp for \n\n");
108//将内层for循环并行化
109#pragma omp parallel
110 for (int i=0;i<2;i++)
111 {
112#pragma omp for
113 for (int j=0;j<2;j++)
114 {
115 printf("Thread %d is active when i=%d,and j=%d\n",omp_get_thread_num(),i,j);
116 }
117 }
118// 8.threadprivate and copyin
119 printf("------------------华丽无比的分隔线--------------------\n");
120 printf("#pragma omp threadprivate and copyin \n\n");
121
122 int global=10;
123//#pragma omp threadprivate(global)
124//#pragma omp parallel
125 for (int i=0;i<100;i++)
126 {
127 global++;
128 }
129
130 printf("Global is equal to %d\n",global);
131
132
133// 9.功能并行#prgama omp parallel sections
134 printf("------------------华丽无比的分隔线--------------------\n");
135 printf("#pragma omp parallel sections \n\n");
136
137#pragma omp parallel sections
138 {
139#pragma omp section
140 printf("Section 1 's thread is NO.%d\n",omp_get_thread_num());
141#pragma omp section
142 printf("Section 2 's thread is NO.%d\n",omp_get_thread_num());
143#pragma omp section
144 printf("Section 3 's thread is NO.%d\n",omp_get_thread_num());
145 }
146
147// 10.同步编译指导语句
148// 1).同步屏障
149#pragma omp parallel
150 {
151 int initialization=0;
152#pragma omp barrier
153 int process=1;
154 }
155
156// 2.互斥锁机制(库函数)
157#ifdef LOCKLOCK
158#define LOCKLOCK
159 omp_lock_t lock;
160//普通锁
161 omp_init_lock(&lock);
162 omp_destroy_lock(&lock);
163 omp_set_lock(&lock);
164 omp_unset_lock(&lock);
165 omp_test_lock(&lock);
166//嵌套锁
167 omp_init_nest_lock(&lock);
168 omp_destroy_nest_lock(&lock);
169 omp_set_nest_lock(&lock);
170 omp_unset_nest_lock(&lock);
171 omp_test_nest_lock(&lock);
172#endif
173//其他的OpenMP编译指导语句请参考OpenMP手册;
174
175
176
177 return 0;
178}
转载于:https://www.cnblogs.com/lyq105/archive/2008/12/19/1358727.html
一个OpenMP的学习程序相关推荐
- 一个openMP编程处理图像的示例
一个openMP编程处理图像的示例: 从硬盘读入两幅图像,对这两幅图像分别提取特征点,特征点匹配,最后将图像与匹配特征点画出来.理解该例子需要一些图像处理的基本知识,我不在此详细介绍.另外,编译该例需 ...
- OpenMP用法大全
OpenMP基本概念 OpenMP是一种用于共享内存并行系统的多线程程序设计方案,支持的编程语言包括C.C++和Fortran.OpenMP提供了对并行算法的高层抽象描述,特别适合在多核CPU机器上的 ...
- OpenMP知识点汇总
1. OpenMP(Open Multi-Processing)官网:http://openmp.org/wp/ 2. OpenMP最新版本4.0,2013年7月发布.Visual Studio 20 ...
- 《OpenMP编译原理及实现技术》摘录
内容摘自<OpenMP编译原理及实现技术>第2章 代码测试环境:Windows7 64bit, VS2010, 4核机. 可以说OpenMP制导指令将C语言扩展为一个并行语言,但OpenM ...
- eclipse openmp mpi并行编程例子
openmp 首先,下载这个版本的Eclipse:http://www.eclipse.org/downloads/packages/eclipse-parallel-application-deve ...
- OpenMP基本概念
OpenMP是一种用于共享内存并行系统的多线程程序设计方案,支持的编程语言包括C.C++和Fortran.OpenMP提供了对并行算法的高层抽象描述,特别适合在多核CPU机器上的并行程序设计.编译器根 ...
- 版本为2.5的OpenMP规范中的编译制导指令汇总
OpenMP的所有编译制导指令是以#pragma omp开始 #pragmaomp指令 [子句[,子句]„] 版本为2.5的OpenMP规范中的指令有以下这些: parallel:用在一个结构块之前, ...
- keil流水灯c语言程序两个一起亮,Keil单片机点亮一个灯及循环流水灯三种实现方法详解...
实验名称:keil工程建立,点亮一个led灯 实验目的:学会keil软件安装,熟悉keil界面并学习如何新建一个工程 实验器材:安装有keil的电脑一台 预习内容及原理: Keil C51已集成到一个 ...
- Linux openmp教程,OpenMP中文教程
1.摘要 OpenMP 是一个应用程序接口(API),由一组主要的计算机硬件和软件供应商联合定义.OpenMP 为共享内存并行应用程序的开发人员提供了一个可移植的.可伸缩的模型.该API在多种体系结构 ...
- openmp多线程简单编程
OpenMP支持的编程语言包括C.C++和Fortran,简单的说,就是一种API,来编写多线程应用程序.通过使用简单的指令#pragma omp -就可以对程序进行多线程并行.OpenMP使得程序员 ...
最新文章
- python读取大数据量xml_[C#]_[使用微软OpenXmlSDK (OpenXmlReader)读取xlsx表格] 读取大数据量100万条数据Excel文件解决方案...
- 开放大学MySQL形考_95至尊考试网-国家开放大学《数据库基础与应用》形考任务1...
- 快速下单!简化EcStore的购物结算流程
- 1043 输出PATest (20分)_23行代码满分
- MFC载入JPG图片
- ant design vue中通知提醒框Notification的使用
- c++ array容器 传参_C++ 顺序容器基础知识总结
- 更改主机名后mysql无法启动_主机名变更导致MySql启动失败
- 软件平台与中间技术复习
- 组件(component)技术介绍
- PNG图片怎么转成ICO?分享两种思路
- Oracle RAC原理
- java 栈队列区别是什么意思_Java实现数据结构栈stack和队列Queue是什么?
- html中样式里面有符号,css样式 + 特殊符号
- Go: panic / recover 简介与实践
- 时下火热的 NFT 究竟有什么用?
- 基于容器的后端服务架构
- 【密码学】古代、古典密码
- 03-docker系列-docker容器的基本操作
- php保存文件快捷键,word保存快捷键是ctrl加什么