性能优化-OpenMP基础教程（一）

本文主要介绍OpenMP并行编程技术，编程模型、指令和函数的介绍、以及OpenMP实战的几个例子。希望给OpenMP并行编程者提供指导。

🎬个人简介：一个全栈工程师的升级之路！
📋个人专栏：高性能（HPC）开发基础教程
🎀CSDN主页 发狂的小花
🌄人生秘诀：学习的本质就是极致重复!

目录

一、OpenMP 简介

二、OpenMP 编程模型

1. 指令与库函数

1.1 OpenMP指令格式

1.1.1 并行区域（Parallel Region）

1.1.2 并行构造（Parallel Construct）

1.1.3 任务（Task）

1.1.4 同步（Synchronize）

1.2 OpenMP常用的指令和函数

1.3 OpenMP常用库函数

2. 并行执行

3. 线程管理

4. 同步与通信

5. 调度策略

三、OpenMP编程实战

1 Linux下编译选项

2 C语言 OpenMP 并行化程序示例（包含）

3 C++ OpenMP并行编程示例（包含宏定义#ifdef _OPENMP）

4 OpenMP 多线程性能对比

一、OpenMP 简介

        OpenMP 是一个为共享内存并行计算设计的编程接口，广泛应用于 Fortran、C 和 C++ 语言。它提供了一套编译器指令和库函数，使得开发者能够轻松地编写并行程序。OpenMP 的“fork/join”模型是其中最核心的并行执行模式，其中最初只有一个主线程在运行。当遇到需要并行计算的部分时，主线程会派生出其他线程来执行并行任务。当并行代码执行完毕，派生的线程会退出或挂起，控制权回到主线程。类似与多线程技术。

二、OpenMP 编程模型

1. 指令与库函数

        OpenMP 的基本语法是通过预处理指令 #pragma omp 来实现的。例如，#pragma omp parallel for 用于并行化 for 循环。此外，OpenMP 还提供了一系列的库函数，用于线程的创建、同步等操作。这些库函数和指令使得开发者能够更灵活地控制并行程序的执行。

1.1 OpenMP指令格式

1.1.1 并行区域（Parallel Region）

        用于指定一个代码块，该代码块将在多个线程上并行执行。

#pragma omp parallel
{// 并行执行的代码块
}
1.1.2 并行构造（Parallel Construct）

        用于创建一个新线程并执行指定的代码块

#pragma omp parallel sections
{#pragma omp section{// 线程1执行的代码块}#pragma omp section{// 线程2执行的代码块}
}
1.1.3 任务（Task）

        用于创建一个新任务并在当前线程上执行指定的代码块。

#pragma omp task firstprivate(a, b) shared(c)
{// 任务执行的代码块，使用变量a和b，以及共享变量c
}
1.1.4 同步（Synchronize）

        用于等待所有线程完成指定的任务。

#pragma omp for schedule(static, chunk_size) reduction(+:sum)
for (int i = 0; i < n; i++) {// 循环体，使用变量i和sum
}

1.2 OpenMP常用的指令和函数

  1. parallel：用于指定一个代码段，该代码段将在多个线程上并行执行。

  2. for：用于for循环之前，将循环分配到多个线程中并行执行，必须保证每次循环之间无相关性。

  3. parallel for：parallel 和 for语句的结合，也是用在一个for循环之前，表示for循环的代码将被多个线程并行执行。

  4. sections：用在可能会被并行执行的代码段之前。

  5. parallel sections：parallel和sections两个语句的结合。

  6. critical：用在一段代码临界区之前。

  7. single：用在一段只被单个线程执行的代码段之前，表示后面的代码段将被单线程执行。

  8. flush：用来保证线程的内存临时视图和实际内存保持一致，即各个线程看到的共享变量是一致的。

  9. barrier：用于并行区内代码的线程同步，所有线程执行到barrier时要停止，直到所有线程都执行到barrier时才继续往下执行。

  10. atomic：用于指定一块内存区域被制动更新。

  11. master：用于指定一段代码块由主线程执行。

  12. ordered：用于指定并行区域的循环按顺序执行。

  13. threadprivate：用于指定一个变量是线程私有的。

  14. copyprivate：配合single指令，将指定线程的专有变量广播到并行域内其他线程的同名变量中；

  15. copyin n：用来指定一个threadprivate类型的变量需要用主线程同名变量进行初始化；

  16. default：用来指定并行域内的变量的使用方式，缺省是shared。

1.3 OpenMP常用库函数

        OpenMP库函数是一组用于并行计算的函数，它们可以帮助程序员在C、C++和Fortran等编程语言中实现多线程编程。以下是一些常用的OpenMP库函数：

  1. omp_get_num_threads()：返回正在执行的线程数。
  2. omp_get_max_threads()：返回支持的最大线程数。
  3. omp_get_thread_num()：返回当前线程的编号。
  4. omp_get_num_procs()：返回正在执行的程序的处理器数。
  5. omp_set_num_threads()：设置并行区域中的线程数。
  6. omp_get_nestled()：测试当前块是否嵌套在其他并行区域内。
  7. omp_set_nestled()：设置当前块是否嵌套在其他并行区域内。
  8. omp_get_schedule()：获取指定并行区域的调度策略。
  9. omp_set_schedule()：设置指定并行区域的调度策略。
  10. omp_get_chunk_size()：获取指定并行区域的块大小。
  11. omp_set_chunk_size()：设置指定并行区域的块大小。
  12. omp_barrier()：在所有线程都到达该点时阻塞所有线程。
  13. omp_critical()：创建一个临界区，确保同一时间只有一个线程可以执行该段代码。
  14. omp_atomic()：对一个变量进行原子操作，确保多个线程对该变量的操作是有序的。
  15. omp_flush()：将缓冲区中的数据立即写入共享内存或设备。
  16. omp_lock_t：用于同步的锁类型。
  17. omp_init_lock()：初始化锁对象。
  18. omp_destroy_lock()：销毁锁对象。
  19. omp_set_lock()：对锁对象加锁。
  20. omp_unset_lock()：对锁对象解锁。

2. 并行执行

        OpenMP 提供了多种并行执行的方法，如 parallel for、parallel sections 等。这些方法使得开发者能够将代码块分配给多个线程执行，从而实现更高效的计算。通过合理地划分代码块和选择合适的并行执行方法，开发者可以显著提高程序的性能。

3. 线程管理

        OpenMP 提供了一些指令和函数，如 num_threads、thread_bind 等，用于设置和控制并行区域中的线程数量和绑定策略。这些功能使得开发者能够更好地控制并行程序的执行流程，确保程序的正确性和稳定性。

4. 同步与通信

        为了确保并行执行的正确性，OpenMP 提供了一些同步机制，如 barrier、critical、atomic 等。这些机制确保了线程之间的正确协作和数据一致性。此外，还提供了一些数据传输函数，如 reduction，用于实现线程之间的数据共享和计算结果的汇总。这些同步和通信机制是并行程序中必不可少的部分，它们确保了程序的正确性和可靠性。

5. 调度策略

        OpenMP 支持多种调度策略，如静态调度、动态调度和运行时调度。这些调度策略允许开发者根据需要选择合适的调度策略来优化程序的性能。通过合理地选择调度策略，开发者可以更好地平衡线程的负载和利用系统资源，从而提高程序的执行效率。

三、OpenMP编程实战

1 Linux下编译选项

        Linux下GCC编译器仅仅编译选项增加-fopenmp即可完成对OpenMP的支持。

2 C语言 OpenMP 并行化程序示例（包含<omp.h>）

#include <omp.h>#include <stdio.h>int main() {#pragma omp parallel forfor (int i = 0; i < 10; i++) {printf("Thread %d: %d\n", omp_get_thread_num(), i);}return 0;}

这个程序使用了 #pragma omp parallel for 指令将 for 循环进行并行化。在循环体内部，使用 omp_get_thread_num() 函数获取当前线程的编号，并打印出来。这个示例展示了 OpenMP 的基本用法和并行化效果，通过简单的修改和调整，你可以将其应用于更复杂的并行计算任务。

        运行结果：

        由于使用的电脑是八核的，因此，最多有八个线程，由上述的线程编号可以看出。

如果将上述的循环代码变成8个，如下：

#include <omp.h>#include <stdio.h>int main() {#pragma omp parallel forfor (int i = 0; i < 8; i++) {printf("Thread %d: %d\n", omp_get_thread_num(), i);}return 0;}

        运行结果：

        运行结果是八个线程，线程编号和循环编号相同。

3 C++ OpenMP并行编程示例（包含宏定义#ifdef _OPENMP）

#include <iostream>
#include <omp.h>
int main()
{#ifdef _OPENMP // 如果定义了这个宏std::cout << "Hello, OpenMP!" << std::endl;#pragma omp parallel forfor (int i = 0;i < 8;i++){printf("thread ID is %d i = %d\n",omp_get_thread_num(),i);}#elsestd::cout << "OpenMP is not enabled." << std::endl;#endifreturn 0;
}

        运行结果：

        C++ OpenMP并行编程例子。-fopenmp编译选项开启后，_OPENMP宏被打开。

4 OpenMP 多线程性能对比

#include <stdlib.h>
#include <stdio.h>
#include "omp.h"void test()
{for (int i = 0; i < 80000; i++){//执行代码}
}int main(int argc, char **argv){#ifdef _OPENMPprintf("OpenMP is Enable!\n");#elseprintf("OpenMP is Disable!\n");#endiffloat startTime = omp_get_wtime();//指定2个Thread
#pragma omp parallel for num_threads(2)for (int i = 0; i < 80000; i++){test();}float endTime = omp_get_wtime();printf("2 个Thread,latency: %f\n", endTime - startTime);startTime = endTime;//指定4个Thread
#pragma omp parallel for num_threads(4)for (int i = 0; i < 80000; i++){test();}endTime = omp_get_wtime();printf("4 个Thread,latency: %f\n", endTime - startTime);startTime = endTime;//指定8个Thread
#pragma omp parallel for num_threads(8)for (int i = 0; i < 80000; i++){test();}endTime = omp_get_wtime();printf("8 个Thread,latency: %f\n", endTime - startTime);startTime = endTime;//指定12个Thread#pragma omp parallel for num_threads(10)for (int i = 0; i < 80000; i++){test();}endTime = omp_get_wtime();printf("10 个Thread,latency: %f\n", endTime - startTime);startTime = endTime;//不使用OpenMPfor (int i = 0; i < 80000; i++){test();}endTime = omp_get_wtime();printf("不使用OpenMP Mutil Thread,latency: %f\n", endTime - startTime);startTime = endTime;return 0;
}

        运行结果：

        分析结果可知，随着线程数量的增加运行的时间减少，由于使用的电脑是八核的，因此并行只能同时有八个线程，使用十个线程的运行效率不增反减。

🌈我的分享也就到此结束啦🌈
如果我的分享也能对你有帮助，那就太好了！
若有不足，还请大家多多指正，我们一起学习交流！
📢未来的富豪们：点赞👍→收藏⭐→关注🔍，如果能评论下就太惊喜了！
感谢大家的观看和支持！最后，☺祝愿大家每天有钱赚！！！

下一节将继续开展OpenMP编程更加详细的实战。

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.hqwc.cn/news/315248.html

如若内容造成侵权/违法违规/事实不符，请联系编程知识网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

前端框架变迁史

前端框架变迁史

前端框架变迁史 过去的前端 JQuery 简单的选择器和链式操作跨浏览器兼容性强大的动画功能简单的 Ajax 支持插件系统 MVVM 模式 什么是 MVVM 什么是 MVVM &#xff1f; 现在的前端 主流技术 react VUE angular 作为现代前端框架&#xff0c;有如下 4 个共有的特点&…
阅读更多...
操作细节无遗漏：上海迅软DSE如何精准记录用户不规范行为

操作细节无遗漏：上海迅软DSE如何精准记录用户不规范行为

企业在部署加密软件后&#xff0c;虽然数据信息得到了加密保护&#xff0c;但难免会发生内部人员滥用管理身份职权&#xff0c;随意解密、外发文件等情况&#xff0c;造成明明对数据加密了&#xff0c;却还是无法完全掌控数据使用情况的现象发生&#xff0c;因此审计功能变得尤…
阅读更多...
抖店和商品橱窗有什么区别？新手应该选哪个？

抖店和商品橱窗有什么区别？新手应该选哪个？

我是电商珠珠 临近年底了&#xff0c;有的人已经开始为下一年筹谋&#xff0c;有的去抖音做账号做直播带货&#xff0c;不会直播带货的就想尝试做下抖店&#xff0c;来为以后的经济打基础。 刚想要接触却对这类有些迷糊&#xff0c;发现商品橱窗和抖店都可以卖货&#xff0c;…
阅读更多...
天融信TOPSEC安全管理系统存在远程命令执行漏洞

天融信TOPSEC安全管理系统存在远程命令执行漏洞

文章目录 产品简介漏洞概述指纹识别漏洞利用修复建议 产品简介 天融信TopSec 安全管理系统&#xff0c;是基于大数据架构&#xff0c;采用多种技术手段收集各类探针设备安全数据&#xff0c;围绕资产、漏洞、攻击、威胁等安全要素进行全面分析&#xff0c;提供统一监测告警、集…
阅读更多...
MySQL基础学习: linux系统mysql 密码插件 validate_password安装

MySQL基础学习: linux系统mysql 密码插件 validate_password安装

1、没有安装mysql密码插件&#xff0c;执行命令&#xff1a;SHOW VARIABLES LIKE ‘validate_password%’; 2、安装mysql密码插件&#xff0c;执行命令&#xff1a;install plugin validate_password soname ‘validate_password.so’; 3、再次执行&#xff1a;SHOW VARIABLE…
阅读更多...
嵌入式视频播放器（mplayer）

嵌入式视频播放器（mplayer）

1.文件准备&#xff1a; MPlayer-1.0rc2.tar.bz2 libmad-0.15.1b.tar.gz 直接Git到本地 git clone https://gitee.com/zxz_FINE/mplayer_tarball.git 2.文件夹准备&#xff1a; src存放解压后的源码文件&#xff0c;target_Mplayer存放编译安装的目标文件 mkdir src targe…
阅读更多...
通过CANopen转Modbus TCP网关的伺服系统

通过CANopen转Modbus TCP网关的伺服系统

通过开疆CANopen转Modbus TCP网关的伺服系统&#xff0c;我们能够实现CANopen协议与Modbus TCP协议之间的无缝转换&#xff0c;从而为伺服电机的通讯和控制带来更大的便利和灵活性。 首先我们可以通过CANopen协议与伺服电机进行通讯。CANopen协议是一种先进的、高效的通讯协议&…
阅读更多...
7+WGCNA+机器学习+泛癌生信思路,非肿瘤也能结合泛癌分析

7+WGCNA+机器学习+泛癌生信思路,非肿瘤也能结合泛癌分析

今天给同学们分享一篇生信文章“Analysis and Experimental Validation of Rheumatoid Arthritis Innate Immunity Gene CYFIP2 and Pan-Cancer”&#xff0c;这篇文章发表在Front Immunol期刊上&#xff0c;影响因子为7.3。 结果解读&#xff1a; DEG筛选和数据预处理 数据在…
阅读更多...
舞动指尖的代码之旅

舞动指尖的代码之旅

嗨&#xff0c;伟大的程序员们&#xff01; 随着2023年的尾声临近&#xff0c;我回顾这一年的使用 Vim 的心得体会&#xff0c;愿与你分享我的点滴成长。 初尝 Vim 之苦与乐 在开始使用 Vim 时&#xff0c;曾有过一些困惑和挫折。但正是在不断的摸索中&#xff0c;我领略到了…
阅读更多...
shell 切片参数解释

shell 切片参数解释

result"set-cookie: versionv1; Max-Age2592000" result${result:20:22}本来只想要 version 结果 v1&#xff0c;第一次使用一直输出这么长&#xff0c;仔细研究才发现第二个参数的含义是切片的长度 改为以下语句后成功输出 version 对应的值 result"set-cook…
阅读更多...
openmediavault(OMV) (26)网络(1)ddns-go

openmediavault(OMV) (26)网络(1)ddns-go

简介 "ddns-go" 是一个动态域名解析（Dynamic DNS）工具，用于更新域名的IP地址。它可以自动检测你的公共IP地址，并将其更新到指定的域名解析服务商，以确保你的域名始终与最新的IP地址相匹配。 安装 hub.docker.com上下载ddns-go镜像 配置compose文件 --- versio…
阅读更多...
【刷题日志】深度理解除（/）与取模（%）附水仙花数以及变种水仙花数题解

【刷题日志】深度理解除（/）与取模（%）附水仙花数以及变种水仙花数题解

文章目录 &#x1f680;前言&#x1f680;除与取模&#x1f680;水仙花数&#x1f680;变种水仙花数 &#x1f680;前言 本专栏文章都直奔刷题主题&#xff0c;阿辉都不会在废话了&#xff0c;加油&#xff0c;少年&#xff01;&#xff01;&#xff01; &#x1f680;除与取…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023