OpenMP 二三

最新推荐文章于 2023-11-01 18:50:05 发布
最新推荐文章于 2023-11-01 18:50:05 发布
阅读量229 收藏 3
点赞数
分类专栏： C++ 并行 文章标签： openmp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_44076667/article/details/116665325
版权

OpenMP

OpenMP

第一部分：基本知识

1、编译制导

  编译制导法：也称为智能编译，它是隐式并行策略的体现，主要是由并行编译系统进行程序表示、控制流的分析、相关分析、优化分析和并行化划分，由相关分析得到方法库管理方案，由优化分析得到知识库管理方案，由并行化划分得到程序重构，从而形成并行程序。

  OpenMP编译制导指令以#pragma omp 开始，后边跟具体的功能指令，格式如：

#pragma omp 指令[子句],[子句] …]

注意：后面接的大括号（{ ）需要另起一行；

1.1、 指令

指令解释
parallel用在一个代码段之前，表示这段代码将被多个线程并行执行
for用于for循环之前，将循环分配到多个线程中并行执行，必须保证每次循环之间无相关性。
parallel forparallel 和 for语句的结合，也是用在一个for循环之前，表示for循环的代码将被多个线程并行执行
parallel sections用在可被并行执行的代码段之前，用于实现多个结构块语句的任务分担，可并行执行的代码段各自用section指令标出（注意区分sections和section）
critical用在一段代码临界区之前
single用在一段只被单个线程执行的代码段之前，表示后面的代码段将被单线程执行
flush保证各个OpenMP线程的数据影像的一致性
barrier用于并行区内代码的线程同步，所有线程执行到barrier时要停止，直到所有线程都执行到barrier时才继续往下执行
atomic用于指定一块内存区域被自动更新
master用于指定一段代码块由主线程执行
ordered用于指定并行区域的循环按顺序执行
task申请一个线程去执行指定区域的代码

1.2、子句

子句解释
private指定一个或多个变量在每个线程中都有它自己的私有副本
firstprivate指定一个或多个变量在每个线程都有它自己的私有副本，并且私有变量要在进入并行域或任务分担域时，继承主线程中的同名变量的值作为初值
lastprivate是用来指定将线程中的一个或多个私有变量的值在并行处理结束后复制到主线程中的同名变量中，负责拷贝的线程是for或sections任务分担中的最后一个线程
reduction用来指定一个或多个变量是私有的，并且在并行处理结束后这些变量要执行指定的归约运算，并将结果返回给主线程同名变量
nowait指出并发线程可以忽略其他制导指令暗含的路障同步
num_threads指定并行域内的线程的数目
schedule指定for任务分担中的任务分配调度类型
shared指定一个或多个变量为多个线程间的共享变量
ordered用来指定for任务分担域内指定代码段需要按照串行循环次序执行
copyprivate配合single指令，将指定线程的专有变量广播到并行域内其他线程的同名变量中
copyin用来指定一个threadprivate类型的变量需要用主线程同名变量进行初始化
default用来指定并行域内的变量的使用方式，缺省是shared

2、库函数

  通过调用OpenMP里面的API实现并行（类似平常使用的函数）。

库函数解释
omp_get_num_procs返回运行本线程的多处理机的处理器个数。
omp_get_num_threads返回当前并行区域中的活动线程个数。
omp_get_thread_num返回线程号
omp_set_num_threads设置并行执行代码时的线程个数
omp_init_lock初始化一个简单锁
omp_set_lock上锁操作
omp_unset_lock解锁操作，要和omp_set_lock函数配对使用。
omp_destroy_lockomp_init_lock函数的配对操作函数，关闭一个锁

第二部分：详细说明

1、#pragma omp parallel

  使用parallel指令只是产生了并行域，在没有其他配合指令的情况下只是让多个线程分别执行相同的任务，没有实际价值。
代码：

#define _CRT_SECURE_NO_WARNINGS 1
#include"omp.h"
#include<iostream>
#include<ctime>
#include<Windows.h>

using namespace std;

int main() {
#pragma omp parallel num_threads(5)
    {
        // 基于当前系统的当前日期/时间
        time_t now = time(0);
        // 把 now 转换为字符串形式
        char* dt = ctime(&now);

        cout << "Thread: " << omp_get_thread_num() << "\tTime: " << dt;
        Sleep(2000);
    }
    return 0;
}

结果：

Thread: 2       Time: Tue May 11 22:35:08 2021
Thread: 4       Time: Tue May 11 22:35:08 2021
Thread: 0       Time: Tue May 11 22:35:08 2021
Thread: 1       Time: Tue May 11 22:35:08 2021
Thread: 3       Time: Tue May 11 22:35:08 2021

2、#pragma omp parallel for

  该语句告知编译器，紧接着的for循环可以被并行执行，每次循环之间不能有关系，另外可以让系统默认分配线程个数，也可以使用num_threads子句指定线程个数。

  • 循环中迭代次数必须在执行循环前就可以算出
  • 循环中不能包含break,return或exit
  • 循环中不能包含前往循环外的goto语句
    代码：
#define _CRT_SECURE_NO_WARNINGS 1
#include"omp.h"
#include<iostream>
#include<time.h>
#define MAX_NUM 100000000

using namespace std;

int main() {
    clock_t start_time = clock();
    for (int i = 0; i < MAX_NUM; i++) {
        i* i;
    }
    clock_t end_time = clock();
    cout << "串行运行时间：" << static_cast<double> (end_time - start_time) / CLOCKS_PER_SEC * 1000 << "ms" << endl;

    start_time = clock();
#pragma omp parallel for
    for(int i = 0; i < MAX_NUM; i++) {
        i* i;
    }
    end_time = clock();
    cout << "并行运行时间：" << static_cast<double> (end_time - start_time) / CLOCKS_PER_SEC * 1000 << "ms" << endl;
    return 0;
}

结果：

串行运行时间：182ms
并行运行时间：30ms

3、#pragma omp sections & #pragma omp parallel sections

  #pragma omp sections 串行执行， #pragma omp parallel sections 才是并行执行

代码：

#define _CRT_SECURE_NO_WARNINGS 1
#include"omp.h"
#include<iostream>
#include<time.h>
#include<Windows.h>

using namespace std;

void thread_time() {
    time_t now = time(0);
    char* dt = ctime(&now);

    cout << "Thread: " << omp_get_thread_num() << "\tTime: " << dt;
    Sleep(2000);
}

int main() {
    omp_set_num_threads(2);

#pragma omp sections
    {
#pragma omp section
        {thread_time(); }
#pragma omp section
        {thread_time(); }
    }

    cout << endl;

#pragma omp parallel sections
    {
#pragma omp section
        {thread_time(); }
#pragma omp section
        {thread_time(); }
    }

    return 0;
}

结果：

Thread: 0       Time: Tue May 11 22:36:47 2021
Thread: 0       Time: Tue May 11 22:36:49 2021

Thread: 0       Time: Tue May 11 22:36:51 2021
Thread: 1       Time: Tue May 11 22:36:51 2021

4、#pragma omp single

  在使用了 #pragma omp single 之后，后面指定区域只有一个线程去执行，即使前面使用了#pragma omp parallel num_threads(5) 指定线程数为5；
代码：使用#pragma omp single

#define _CRT_SECURE_NO_WARNINGS 1
#include"omp.h"
#include<iostream>
#include<time.h>
#include<Windows.h>

using namespace std;

void thread_time() {
    time_t now = time(0);
    char* dt = ctime(&now);

    cout << "Thread: " << omp_get_thread_num() << "\tTime: " << dt;
}

int main() {
#pragma omp parallel num_threads(5)
    {
#pragma omp single
        {thread_time(); }
    }

    return 0;
}

结果：

Thread: 0       Time: Tue May 11 22:43:59 2021

代码：去掉#pragma omp single

#define _CRT_SECURE_NO_WARNINGS 1
#include"omp.h"
#include<iostream>
#include<time.h>
#include<Windows.h>

using namespace std;

void thread_time() {
    time_t now = time(0);
    char* dt = ctime(&now);

    cout << "Thread: " << omp_get_thread_num() << "\tTime: " << dt;
}

int main() {
#pragma omp parallel num_threads(5)
    {
        thread_time();
    }

    return 0;
}

结果：

Thread: 1       Time: Tue May 11 22:45:31 2021
Thread: 3       Time: Tue May 11 22:45:31 2021
Thread: 0       Time: Tue May 11 22:45:31 2021
Thread: 4       Time: Tue May 11 22:45:31 2021
Thread: 2       Time: Tue May 11 22:45:31 2021

5、#pragma omp task

  在递归算法中，可以利用 #pragma omp task 不断申请线程去执行递归程序；

#include"omp.h"
#include<iostream>
#include<unistd.h>
#include<ctime>

using namespace std;

void task(int parent_thread, int count) {
    if(count > 2)
        return;
    count++;

    time_t timel;
    time(&timel);
    int child_thread = omp_get_thread_num();

    cout << "Thread " << child_thread << " from Thread " << parent_thread << "\tTime: " << asctime(gmtime(&timel));
    sleep(1);

#pragma omp task
    {task(child_thread, count); }
#pragma omp task
    {task(child_thread, count); }
}

int main(int argc, char *argv[]){
    int count = 0;
#pragma omp parallel num_threads(5)
    {
#pragma omp single
        {task(omp_get_thread_num(), count); }
    }

    return 0;
}

结果：

Thread 1 from Thread 1	Time: Wed May 12 11:41:52 2021
Thread 1 from Thread 1	Time: Wed May 12 11:41:53 2021
Thread 2 from Thread 1	Time: Wed May 12 11:41:53 2021
Thread 1 from Thread 1	Time: Wed May 12 11:41:54 2021
Thread 0 from Thread 1	Time: Wed May 12 11:41:54 2021
Thread 2 from Thread 2	Time: Wed May 12 11:41:54 2021
Thread 4 from Thread 2	Time: Wed May 12 11:41:54 2021

从结果可知：第一秒，线程1执行task；第二秒，在 #pragma omp task 处申请了一个新线程，两个线程并行执行；第三秒，两个线程分别申请一个新线程，四个线程并行。

参考文献

[1]初步认识并行计算openmp – elimuzi
[2]OpenMp之sections用法 – 很厉害的名字
[3]并行与分布式计算导论（五）OpenMP基础详解 - 大选帝侯
[4]【OpenMP】parallel for使用多线程进行并行加速 - 颜丑文良777

yimik
关注
专栏目录
性能优化-OpenMP基础教程（三）
发狂的小花的博客
01-03 1958
本文主要介绍OpenMP并行编程的环境变量和实战、主要对比理解嵌套并行的效果。🎬：一个全栈工程师的升级之路！🎀CSDN主页🌄人生秘诀：学习的本质就是极致重复!
OpenMP使用详解
最新发布
程序猿老樊的博客
06-23 1152
OpenMP 提供了简单而强大的并行编程能力，适合在共享内存系统上实现高效的并行计算。通过合理使用 OpenMP 指令，可以显著提高程序的执行效率。
编译原理-语法制导翻译、后缀表达式、三元、四元
宇丁加
06-19 6002
此时的每条规则后面会对应一条语义动作，即制导翻译。 一、中间语言 中间语言就是和源程序等价的一种编码方式，复杂性也是介于源程序和机器语言中间。
openmp #pragma omp task总结
weixin_44604698的博客
11-01 794
如果产生task的数量过多且每个task处理的内容很少，会发生极其严重的负向优化；在unbloundloop结构中，如果process的的内容与前后的task有关，也是负向优化；在unbloundloop结构中，只有每个task处理的内容与前后无关，会产生正向优化；taskloop可以有效降低运行时间，但是grainsize不应该选用过小。
【openMP】single （主从模式）
xll_bit的博客
10-21 2892
omp single 指定在当前的并行域里当前代码块只被执行一次（任意线程均可） 用法： #pragma omp single[clause,clause, ...] 参数： clause Can be one or more of the following clauses: copyprivate(list) Provides a mechanism to use a privat..
并行与分布式计算导论（五）OpenMP基础详解
Kaiser_syndrom的博客
04-01 2485
并行与分布式计算：OpenMP详解（六）OpenMP编译指导语句指导句：#pragma omp parallel for指导句：#pragma omp parallel指导句：#pragma omp for指导句：#pragma omp single指导句：#pragma omp section(s)指导句：#pragma omp task指导句：#progma omp barrier 与 #pr...
vs上C/C++并行计算#pragma omp
u011089523的博客
10-16 2618
在一个vs内的工程进行并行计算，首先先修改属性内的C/C++ — 语言—OpenMP:是。 头文件不一定需要#include <omp.h>。举例常见的几个用法 section： #pragma omp parallel sections//告诉编译器有几个部分需要使用并行算法 { #pragma omp section//第一块并行部分 {//} #pragma...
eclipse 下配置OpenMP包
10-01
### 三、在Eclipse中配置OpenMP #### 3.1 准备工作 - **安装Eclipse**：首先确保已经安装了最新版本的Eclipse IDE for C/C++ Developers。 - **安装OpenMP库**：下载并安装适合您操作系统的OpenMP库。通常，OpenMP...
论文研究-基于OpenMP的三维并行Delaunay网格生成算法及实现.pdf
07-22
针对大多数并行Delaunay网格生成算法无法充分利用多核共享内存结构的问题，在原有面向共享内存的二维并行算法基础上，根据三维问题的特点提出基于OpenMP的三维并行Delaunay网格生成算法。算法采用划分求解区域为方格...
实现（串行，openmp、mpi、openmp+mpi）快速排序算法，并作出时间对比图
12-26
（5）(选做)在实验报告中通过图表说明三种并行程序使用不同的数据分配方法在各种规模的运行时间。 设计思路 步骤一： 主要采用快速排序实现（串行，openmp、mpi、openmp+mpi）排序算法，所需环境为VS2019+openmp+...
OpenMP3.0的新特性Task指令基础
。。。。
11-23 8960
从OpenMP3.0开始，OpenMP增加了task指令，这是OpenMP3.0中最激动人心的一个新特性。 本文的”术语“大多数是根据个人理解用词，不保证用词准确性。 （1）task基础 OpenMP Tutorials: 任务构造https://computing.llnl.gov/tutorials/openMP/#Task 从功能上说： 1. The TASK construc
【OpenMp】openmp库的基本语法
种一棵树最好的时间是十年前，其次是现在
08-23 1186
动态调度模式是先到先得的方式进行任务分配，不用size参数的时候，先把任务干完的线程先取下一个任务，以此类推，而不是一开始就分配固定的任务数。总之，reduction指令是OpenMP中用于对并行循环中的变量进行归约操作的指令，能够自动将每个线程的局部变量的结果累加到一个全局变量中。在并行循环中，每个线程都会拥有自己的sum的局部副本，并在循环的迭代过程中进行累加操作。当循环结束后，reduction指令将会将每个线程的sum的局部副本的结果相加，并将最终的总和存储到全局的sum变量中。
OpenMP编程（4）—sections、single指令
yiguagua的博客
07-01 3152
1. sections指令 1.1 sections指令用途 SECTIONS指令用于非迭代的多线程共享区。它指定各个section代码段分配给一组线程中部分线程。 多个独立的section指令嵌套在sections指令中，每个section由于其中一个线程执行一次。不同的section可以由不同的线程执行。当然对于一个线程来说，如果它运行足够快，是有可能执行多个section。 1.2 sections语法格式 #pragma omp sections [clause ...] newline
OPENMP学习笔记（3）——命令
qq_23858785的博客
07-20 2621
OPENMP使用 paralell parallel表示其后语句将被多个线程并行执行，“#pragma omp parallel”后面的语句（或者，语句块）被称为parallel region。 多个线程的执行顺序是不能保证的。 for 对一个计算量庞大的任务进行划分，让多个线程分别执行计算任务的某一部分，从而达到缩短计算时间的目的。关键是，每个线程执行的计算互不相同（操作的数据不同或者计算任务本...
并行编程OpenMP基础及简单示例
weixin_30832143的博客
12-24 6634
OpenMP基本概念 OpenMP是一种用于共享内存并行系统的多线程程序设计方案，支持的编程语言包括C、C++和Fortran。OpenMP提供了对并行算法的高层抽象描述，特别适合在多核CPU机器上的并行程序设计。编译器根据程序中添加的pragma指令，自动将程序并行处理，使用OpenMP降低了并行编程的难度和复杂度。当编译器不支持OpenMP时，程序会退化成...
OpenMP编译器指令
weixin_41860751的博客
02-11 2412
1. 编译器指令 OpenMP通过在串行程序中插入编译制导指令, 来实现并行化, 支持OpenMP的编译器可以识别, 处理这些指令并实现对应的功能. 所有的编译制导指令都是以#pragma omp开始, 后面跟具体的功能指令(directive)或者命令. 一般格式如下所示: #pragma omp directive [clause [[,] clause]...] structured bloc Parallel Construct(并行域结构) 为了使程序可以并行执行, 我们首先要构造一个并
for编译制导语句
koastal的博客
12-27 1221
#include #include #define N 100 #define CHUNKSIZE 10 double a[N], b[N], c[N]; int main () { int i, chunk; /* Some initializations */ for (i=0; i < N; i++) a[i]
编译原理：语义分析之语法制导的定义的简单实现
qq_42689684的博客
11-28 2143
语法制导的定义：为每一个产生式编制一个语义子程序。当一个产生式获得匹配时，调用相应的语义子程序实现语义检查和翻译。 可以类比递归下降的预测分析 要求：实现简易计算器 其实只需要在LR分析器的基础上作适当修改 代码 #语义分析，计算器 import pandas as pd data={'id':['s5',' ',' ',' ','s5',' ','s5','s5','s5',...
编译原理-语法制导翻译
WhereIsMyChair的博客
03-19 2201
在学习自制编译器过程中，进行到这一步，不知道如何进行下去。以下是一些个人总结。一、语法制导翻译的任务编译器的前端对输入流进行词法分析、语法分析之后，得到合乎词法规则和语法规则的抽象语法树，他并不能给出程序正确的答案，需要进行语义分析和翻译，才能告诉后端该程序段是具体含义，由后端完成结果计算并给出答案。因此，语法制导翻译是为了理解程序段的语义而设计的阶段，用来将语法树翻译成后端能够理解的表示形式。二...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值