性能优化-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_nested()：测试当前块是否嵌套在其他并行区域内。
  7. omp_set_nested()：设置当前块允许嵌套在其他并行区域内。
  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.mzph.cn/news/595084.shtml

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

相关文章

Mysql count统计去重的数据

Mysql count统计去重的数据

不去重&#xff0c;是4 &#xff1a; SELECT COUNT(NAME) FROM test2 明显里面包含了2个 name 等于 mike的数据&#xff0c; 所以需要做去重 &#xff1a; 通过结合 count 函数和 DISTINCT 关键字 SELECT COUNT(DISTINCT NAME) FROM test2 好了就到这。
阅读更多...
消息中间件 —— ActiveMQ 使用及原理详解

消息中间件 —— ActiveMQ 使用及原理详解

目录 一. 前言 二. JMS 规范 2.1. 基本概念 2.2. JMS 体系结构 三. ActiveMQ 使用 3.1. ActiveMQ Classic 和 ActiveMQ Artemis 3.2. Queue 模式&#xff08;P2P&#xff09; 3.3. Topic 模式&#xff08;Pub/Sub&#xff09; 3.4. 持久订阅 3.5. 消息传递的可靠性 …
阅读更多...
数模学习day06-主成分分析

数模学习day06-主成分分析

主成分分析(Principal Component Analysis,PCA)主成分分析是一种降维算法&#xff0c;它能将多个指标转换为少数几个主成分&#xff0c;这些主成分是原始变量的线性组合&#xff0c;且彼此之间互不相关&#xff0c;其能反映出原始数据的大部分信息。一般来说当研究的问题涉及到…
阅读更多...
当hashCode相同时，equals是否也相同？

当hashCode相同时，equals是否也相同？

在Java中&#xff0c;理解对象的这两个基本方法—hashCode和equals对于编码是至关重要的&#xff0c;尤其是在处理集合类如HashMap和HashSet时。然而&#xff0c;一个常见的误解是&#xff0c;如果两个对象有相同的哈希码&#xff08;hashCode&#xff09;&#xff0c;那么它们…
阅读更多...
iec104和iec61850

iec104和iec61850

iec104和iec61850 IEC104 规约详细解读(一) 协议结构 IEC104 规约详细解读(二)交互流程以及协议解析 61850开发知识总结与分享【1】 Get the necesarry projects next to each other in the same directory; $ git clone https://github.com/robidev/iec61850_open_server.g…
阅读更多...
ES（Elasticsearch）的基本使用

ES（Elasticsearch）的基本使用

一、常见的NoSQL解决方案 1、redis Redis是一个基于内存的 key-value 结构数据库。Redis是一款采用key-value数据存储格式的内存级NoSQL数据库&#xff0c;重点关注数据存储格式&#xff0c;是key-value格式&#xff0c;也就是键值对的存储形式。与MySQL数据库不同&#xff0…
阅读更多...
DNS安全与访问控制

DNS安全与访问控制

一、DNS安全 1、DNSSEC原理 DNSSEC依靠数字签名保证DNS应答报文的真实性和完整性。权威域名服务器用自己的私有密钥对资源记录&#xff08;Resource Record, RR&#xff09;进行签名&#xff0c;解析服务器用权威服务器的公开密钥对收到的应答信息进行验证。如果验证失败&…
阅读更多...
数字信号处理期末复习——计算小题（二）

数字信号处理期末复习——计算小题（二）

个人名片&#xff1a; &#x1f981;作者简介&#xff1a;一名喜欢分享和记录学习的在校大学生 &#x1f42f;个人主页&#xff1a;妄北y &#x1f427;个人QQ&#xff1a;2061314755 &#x1f43b;个人邮箱&#xff1a;2061314755qq.com &#x1f989;个人WeChat&#xff1a;V…
阅读更多...
图像清晰度评估指标

图像清晰度评估指标

图像清晰度评估涉及多个指标&#xff0c;这些指标可用于定量测量图像的清晰度和质量。 以下是一些常见的图像清晰度评估指标&#xff1a; 均方根误差&#xff08;Root Mean Square Error&#xff0c;RMSE&#xff09;&#xff1a; 通过计算原始图像和处理后图像之间的像素差异的…
阅读更多...
【实用工具】Gradio快速部署深度学习应用1：图像分类

【实用工具】Gradio快速部署深度学习应用1：图像分类

前言 在AI快速发展的今天&#xff0c;我们作为算法开发人员&#xff0c;也应该有一些趁手的工具帮助我们快速开发并验证自己的想法&#xff0c;Gradio可以实现快速搭建和共享的功能&#xff0c;能够展示出一个前端界面&#xff0c;把我们的算法包裹起来&#xff0c;快速验证算…
阅读更多...
DNS主从服务器、转发（缓存）服务器

DNS主从服务器、转发（缓存）服务器

一、主从服务器 1、基本含义 DNS辅助服务器是一种容错设计&#xff0c;考虑的是一旦DNS主服务器出现故障或因负载太重无法及时响应客户机请求&#xff0c;辅助服务器将挺身而出为主服务器排忧解难。辅助服务器的区域数据都是从主服务器复制而来&#xff0c;因此辅助服务器的数…
阅读更多...
Vue2多入口，独立打包配置

Vue2多入口，独立打包配置

提示&#xff1a;Vue2多入口&#xff0c;独立打包配置 文章目录 前言一、修改二、打包和效果预览三、具体操作的文件四、代码包总结 前言 需求&#xff1a;富文本编辑器 一、修改 1、复制&#xff1a;index.html文件并改名share.html。 路径&#xff1a;工程文件夹/index.html …
阅读更多...
Android studio BottomNavigationView 应用设计

Android studio BottomNavigationView 应用设计

一、新建Bottom Navigation Activity项目&#xff1a; 二、修改bottom_nav_menu.xml: <itemandroid:id"id/navigation_beijing"android:icon"drawable/ic_beijing_24dp"android:title"string/title_beijing" /><itemandroid:id"i…
阅读更多...
Excel 插件：ASAP Utilities Crack

Excel 插件：ASAP Utilities Crack

ASAP Utilities是一款功能强大的 Excel 插件&#xff0c;填补了 Excel 的空白。在过去的 20 年里&#xff0c;我们的加载项已经发展成为世界上最受欢迎的 Microsoft Excel 加载项之一。 ASAP Utilities 中的功能数量&#xff08;300 多个&#xff09;可能看起来有点令人眼花缭乱…
阅读更多...
PNG图片导入Abaqus建模：Abaqus Image To Part 2D插件

PNG图片导入Abaqus建模：Abaqus Image To Part 2D插件

插件介绍 Abaqus Image To Part 2D - AbyssFish 插件可将图像导入Abaqus内并通过对网格单元集进行材料指定&#xff0c;实现基于图像的模型部件生成。 插件支持JPEG、JPG、PNG、GIF、TIFF、BMP、PCX、ICO等多种图像格式&#xff0c;兼容彩图、灰度图、二值图像等类型&#x…
阅读更多...
一起学docker（六）| Dockerfile自定义镜像 + 微服务模块实战

一起学docker（六）| Dockerfile自定义镜像 + 微服务模块实战

DockerFile 是什么 Dockerfile是用来构建Docker镜像的文本文件&#xff0c;是由一条条构建镜像所需的指令和参数构成的脚本。 构建步骤 编写Dockerfile文件docker build命令构建镜像docker run运行镜像 Dockerfile构建过程 基础知识 每个保留字指令都必须为大写字母且后面…
阅读更多...
5.云原生安全之ingress配置域名TLS证书

5.云原生安全之ingress配置域名TLS证书

文章目录 cloudflare配置使用cloudflare托管域名获取cloudflare API Token在cloudflare中配置SSL/TLS kubesphere使用cert-manager申请cloudflare证书安装证书管理器创建Secret资源创建cluster-issuer.yaml创建cert.yaml申请证书已经查看申请状态 部署harbor并配置ingress使用证…
阅读更多...
Alibaba Cloud Linux镜像操作系统超详细测评！兼容CentOS

Alibaba Cloud Linux镜像操作系统超详细测评！兼容CentOS

Alibaba Cloud Linux是基于龙蜥社区OpenAnolis龙蜥操作系统Anolis OS的阿里云发行版&#xff0c;针对阿里云服务器ECS做了大量深度优化&#xff0c;Alibaba Cloud Linux由阿里云官方免费提供长期支持和维护LTS&#xff0c;Alibaba Cloud Linux完全兼容CentOS/RHEL生态和操作方式…
阅读更多...
PET塑料粘接时，要求强力粘接性能，那么怎么选胶呢？

PET塑料粘接时，要求强力粘接性能，那么怎么选胶呢？

聚对苯二甲酸乙二醇酯 简称PET。 需要强力粘接塑料PET时&#xff0c;可以选择以下几种胶水&#xff1a; 1.环氧树脂胶 具有较高的强度&#xff0c;硬度和耐久性&#xff0c;能够有效地粘合PET材料。但是&#xff0c;对于某些特殊环境和温度条件&#xff0c;可能需要选择耐高…
阅读更多...
QT_01 安装、创建项目

QT_01 安装、创建项目

QT - 安装、创建项目 1. 概述 1.1 什么是QT Qt 是一个跨平台的 C图形用户界面应用程序框架。 它为应用程序开发者提供建立艺术级图形界面所需的所有功能。 它是完全面向对象的&#xff0c;很容易扩展&#xff0c;并且允许真正的组件编程。 1.2 发展史 1991 年 Qt 最早由奇…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023