【OpenMP学习笔记】1:基本使用和critical指令

2023-12-09 14:48

本文主要是介绍【OpenMP学习笔记】1:基本使用和critical指令,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

MPI是可以针对分布式内存,在进程级别实现并行的API;OpenMP则是针对共享内存,在线程级别实现并行的API。

基本使用

不同于MPI的init和finalize,OpenMP用携带了parallel指令的预处理指令指示接下来的一个代码块被多个线程执行。

OpenMP预处理指令-(携带)->OpenMP指令-(携带)->OpenMP指令的子句。

隐式路障是OpenMP的一个特点,所有执行这个代码块的线程在同一个线程组中以隐式路障隐式同步。

#include<stdio.h>
#include<stdlib.h>//检查是否定义了_OPENMP预处理宏,即编译器是否支持OpenMP
#ifdef _OPENMP 
#include<omp.h> //支持时包含OpenMP头文件
#endifint main(int argc,char *argv[])
{//stdlib的strtol()函数把字符串转换成整形数//参数1:字符串首地址,跳过前面的空格符,遇到数字或正负符号才开始做转换//参数2:用于引用返回遇到非数字或'\0'结束转换时的字符指针所在地址//参数3:转换至几进制,取值2~36或0,取0时即默认为10进制int thrdCnt=strtol(argv[1],NULL,0);printf("strtol()函数从主函数参数解析得的线程数:%d\n",thrdCnt);
//# pragma omp指示后面接的是OpenMP的预处理指令//parallel指令表明之后的结构化代码块被多个线程执行//num_threads(线程数)子句指定执行后代码块的线程数
#   pragma omp parallel num_threads(thrdCnt){//如果不支持OpenMP,前面没有包含omp.h,则里面的函数也不可用//所以要在所有使用了omp.h的代码处做包含检查
#       ifdef _OPENMP//omp_get_thread_num()函数获得该线程的编号int myThrdNm=omp_get_thread_num();//omp_get_num_threads()函数获得总的线程数int allNmThrds=omp_get_num_threads();
#       else//当不能使用OpenMP时,只有一个线程,线程号是0int myThrdNm=0;//当不能使用OpenMP时,只有一个线程,线程数目是1int allNmThrds=1;
#       endif//为了验证"隐式路障",此处sleep()不同的时间,不妨按线程编号延时sleep(myThrdNm*3);//C语言里sleep()的单位是秒printf("我是%d/%d\n",myThrdNm,allNmThrds);}//隐式路障:完成代码块的线程要等待其它所有线程完成此代码块printf("路障解除了\n");return 0;
}

输出

[lzh@hostlzh OpenMP]$ !gcc
gcc -fopenmp -o hello.o hello.c
[lzh@hostlzh OpenMP]$ ./hello.o 4
strtol()函数从主函数参数解析得的线程数:4
我是0/4
我是1/4
我是2/4
我是3/4
路障解除了
[lzh@hostlzh OpenMP]$

critical指令

OpenMP里变量的作用域可以按照块内是否共享分为共享作用域私有作用域。对于并行代码块而言,在块前声明的变量缺省作用域是共享的(不过这个”共享”当然不能穿透函数)。

某些地方未使用critical时可能存在的问题

当块内共同操作了共享的资源时,不对其做互斥保护就可能会在运行时出问题。因为只能保证一次汇编级的指令是原子的,甚至不能保证一条C语言语句在并发执行过程中不会发生线程的切换,在并行的情况下就更加危险了:

#include<stdio.h>
#include<stdlib.h>
#include<omp.h>int main(int argc,char *argv[])
{int a=20;int thrdCnt=strtol(argv[1],NULL,10);
#   pragma omp parallel num_threads(thrdCnt)    a=a-1;//并行块中访问了共享变量aprintf("a=%d\n",a);//最终输出看一下return 0;
}

输出

[lzh@hostlzh OpenMP]$ !gcc
gcc -fopenmp -o test1.o test1.c
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=6
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5

才执行了这么几次就发生了出错,可见问题严重。

使用critical时
#include<stdio.h>
#include<stdlib.h>
#include<omp.h>int main(int argc,char *argv[])
{int a=20;int b=20;int thrdCnt=strtol(argv[1],NULL,10);
#   pragma omp parallel num_threads(thrdCnt){//OpenMP指令critical安排线程对接下来的代码块互斥访问
#       pragma omp criticala=a-1;//共享变量a受critical保护互斥访问b=b-1;//共享变量b未受critical保护}printf("a=%d\nb=%d\n",a,b);//最终输出看一下return 0;
}

输出

[lzh@hostlzh OpenMP]$ !gcc
gcc -fopenmp -o test1.o test1.c
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=6
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ ./test1.o 15
a=5
b=5
[lzh@hostlzh OpenMP]$ 

可见受critical指令保护的a能够正常执行了,而不受保护的b仍然会发生前面的问题。

这篇关于【OpenMP学习笔记】1:基本使用和critical指令的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/474087

相关文章

Linux中的计划任务(crontab)使用方式

《Linux中的计划任务(crontab)使用方式》:本文主要介绍Linux中的计划任务(crontab)使用方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、前言1、linux的起源与发展2、什么是计划任务(crontab)二、crontab基础1、cro

kotlin中const 和val的区别及使用场景分析

《kotlin中const和val的区别及使用场景分析》在Kotlin中,const和val都是用来声明常量的,但它们的使用场景和功能有所不同,下面给大家介绍kotlin中const和val的区别,... 目录kotlin中const 和val的区别1. val:2. const:二 代码示例1 Java

C++变换迭代器使用方法小结

《C++变换迭代器使用方法小结》本文主要介绍了C++变换迭代器使用方法小结,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录1、源码2、代码解析代码解析:transform_iterator1. transform_iterat

C++中std::distance使用方法示例

《C++中std::distance使用方法示例》std::distance是C++标准库中的一个函数,用于计算两个迭代器之间的距离,本文主要介绍了C++中std::distance使用方法示例,具... 目录语法使用方式解释示例输出:其他说明:总结std::distance&n编程bsp;是 C++ 标准

vue使用docxtemplater导出word

《vue使用docxtemplater导出word》docxtemplater是一种邮件合并工具,以编程方式使用并处理条件、循环,并且可以扩展以插入任何内容,下面我们来看看如何使用docxtempl... 目录docxtemplatervue使用docxtemplater导出word安装常用语法 封装导出方

Linux换行符的使用方法详解

《Linux换行符的使用方法详解》本文介绍了Linux中常用的换行符LF及其在文件中的表示,展示了如何使用sed命令替换换行符,并列举了与换行符处理相关的Linux命令,通过代码讲解的非常详细,需要的... 目录简介检测文件中的换行符使用 cat -A 查看换行符使用 od -c 检查字符换行符格式转换将

使用Jackson进行JSON生成与解析的新手指南

《使用Jackson进行JSON生成与解析的新手指南》这篇文章主要为大家详细介绍了如何使用Jackson进行JSON生成与解析处理,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录1. 核心依赖2. 基础用法2.1 对象转 jsON(序列化)2.2 JSON 转对象(反序列化)3.

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

Elasticsearch 在 Java 中的使用教程

《Elasticsearch在Java中的使用教程》Elasticsearch是一个分布式搜索和分析引擎,基于ApacheLucene构建,能够实现实时数据的存储、搜索、和分析,它广泛应用于全文... 目录1. Elasticsearch 简介2. 环境准备2.1 安装 Elasticsearch2.2 J

使用C#代码在PDF文档中添加、删除和替换图片

《使用C#代码在PDF文档中添加、删除和替换图片》在当今数字化文档处理场景中,动态操作PDF文档中的图像已成为企业级应用开发的核心需求之一,本文将介绍如何在.NET平台使用C#代码在PDF文档中添加、... 目录引言用C#添加图片到PDF文档用C#删除PDF文档中的图片用C#替换PDF文档中的图片引言在当