Fortran Coder

查看: 6167|回复: 3

[并行] OpenMP并行采用reduction,出现stack overflow问题

[复制链接]

3

帖子

2

主题

0

精华

入门

F 币
122 元
贡献
59 点

规矩勋章

发表于 2015-11-9 19:06:50 | 显示全部楼层 |阅读模式
利用OpenMP进行并行运算2个矩阵P与Q的值,并利用reduction进行归约处理;
先尝试了一个简单算例,计算没有问题;但是移至到另外一个大程序中时,出现了stack overflow问题。
以下是简单算例的程序部分:
[Fortran] 纯文本查看 复制代码
    ! 主程序段
    matrix = 0.d0
    vector = 0.d0
    !$omp parallel private(i, tmp_vector) &
    !$omp          reduction(+:matrix, vector)
    !$omp do              
            do i =1, 48
            call gemm(a, b, matrix, 'n', 'n', 1.d0, 1.d0)
            call gemm(a, c, tmp_vector)
            call gemm(b, tmp_vector, vector, 'n', 'n', 1.d0, 1.d0)
            end do
   !$omp end do 
    !$omp end parallel


在上述简单算例中运行正确,但移至到大型程序中出现以下问题:


Stack overflow 1.png
stack overflow.PNG

953

帖子

0

主题

0

精华

大师

F 币
180 元
贡献
73 点

规矩勋章元老勋章新人勋章水王勋章热心勋章

QQ
发表于 2015-11-9 20:09:02 | 显示全部楼层
OpenMP 可以一定程度解决效率问题。但并不能解决数组太大内存(堆栈)不足的问题。

29

帖子

3

主题

0

精华

熟手

F 币
214 元
贡献
112 点
发表于 2016-3-22 09:15:16 | 显示全部楼层
好像并行编程有 OpenMP 和MPI两种,本人很想学习,但苦于资料比较少

29

帖子

3

主题

0

精华

熟手

F 币
214 元
贡献
112 点
发表于 2016-3-22 09:17:54 | 显示全部楼层
找到了

呵呵,http://fcode.cn/resource_ebook-list-1.html
慢慢啃,入手并行编程!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-3-29 18:56

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表