并行计算工具OpenMP 二

如果你想要变量是每个线程拥有自己的单独的变量,那么就在并行块里面定义变量,对于在进入并行块前已经定义了的变量,是在多个线程之间共享的,需要小心处理,否则程序不会按照我们想要的方式执行。

并行块前定义变量

在上面这个程序中,试图用sum来累加变量和。#pragma omp parallel for告诉编译器要把这个for语句拆开来并行执行。注意 sum 定义在并行块外面,所以在并行块中会共享这个变量。也就是说在四个核中用的是同一个变量,对于sum + = i 语句,是先把sum的值取出来和 i 相加后再把结果赋给sum变量。既然四个线程用的是同一个sum,那么如果sum在线程0中取出来了,另一个线程中也同时取出来了同样值的sum,加完后值在两个线程中都会被送就回到sum,那么就会造成有一次操作被重写了(数据碰撞),也就是说最后的结果很可能小于正确的结果。(实际上我在visual studio2019上运行上面的代码并没有出错,但是不排除是因为vs进行了更正操作)结果如下:

并行结果

为了避免上面中出现的情况,可以进行一些其他操作:比如可以在 sum += i; 前面加上#pragma omp critical 语句,这条语句是告诉编译器下面的代码块很关键,任何时刻都只能有一个线程执行该代码块。不过这样的话代码会在多个核中的运行会相互影响,对性能有一定影响,最好的解决方法是用规约reduction,代码如下图:

规约reduction代码

reduction (+:sum) 告诉编译器,这并行块要用规约,+表示求和操作,sum是目标变量名。这样的话,编译器就会给每个线程一个sum的拷贝并正确初始化为0,然后每个线程执行完之后再合并。规约只适用于固定的操作符。

规约后代码的结果

OpenMP支持的规约操作符如下图:

规约操作符

如果我们不想要已经定义了的变量在多个线程中共享,想要让每个线程有自己的拷贝怎么办呢?下面就简单介绍一个private从句,以及firstprivate和lastprivate。通过使用#pragma omp parallel private(variable list),告诉编译器在下面的并行块中每个线程都对variable list列出的变量进行拷贝。不过值得注意的时,每个线程中拷贝的变量的初始值是不确定的,执行完并行块后该变量的值也不确定。所以就有了firstprivate和lastprivate,firstprivate 是说变量的初始值就是按照该变量进入并行块之前的值初始化;lastprivate 是说退出并行块后,按照串行情况最后的一个值给到主线程中的变量。如果同时使用firstprivate和lastprivate,那么该变量就会有和串行是一样的初始值和结束值。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,254评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,875评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,682评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,896评论 1 285
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,015评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,152评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,208评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,962评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,388评论 1 304
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,700评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,867评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,551评论 4 335
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,186评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,901评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,142评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,689评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,757评论 2 351

推荐阅读更多精彩内容

  • Swift1> Swift和OC的区别1.1> Swift没有地址/指针的概念1.2> 泛型1.3> 类型严谨 对...
    cosWriter阅读 11,094评论 1 32
  • OpenMP2.5 有底层API后,就已经可以实现并行编程;然而,很多时候串行算法已经成型,如果继续使用原有的底层...
    一梦换须臾_阅读 1,396评论 0 0
  • 本篇文章只是记录api的用法和回顾,方便记忆 openMP openMP提供“基于指令”的共享内存API。这就意味...
    Gavinjou大笨象阅读 13,292评论 0 7
  • What's OpenMP The OpenMP API supports multi-platform shar...
    bitzoo阅读 2,824评论 0 0
  • 人生第一次技术面试:紧张到拿到面试通知到面试那天都无法呼吸写代码:两个题,两个题的答案等我今晚补上 一次遍历,链表...
    Allen的光影天地阅读 2,206评论 0 0