fortran调用子程序的并行调用问题求解答

didi_dudu · 发表于 Post on 2019-8-2 10:34:29

请教一下论坛各位大佬
fortran中调用子程序能不能实现并行调用呢？
比如我有一段代码
do i = 1, 100
call program1
call program2
enddo

call program3

program1和program2都是独立的子程序，各自运算不同的东西。然后我想要program1和program2同时执行，然后都执行完之后再执行program3.有没有什么操作方式？
想用这种方式来实现数据的存储和运算的并行执行。program1是用来读取数据，program2是用来计算上一步program1读取的数据。

万里云 · 发表于 Post on 2019-8-2 11:49:28

rogram1是用来读取数据，program2是用来计算上一步program1读取的数据。

这个没法并行。不读完数据怎么处理？

niobium · 发表于 Post on 2019-8-2 12:54:11

你又说两个程序是独立的，又说2需要1的数据。

像你说的想同时运行1和2，有一个简单的判断方法，如果先1后2和先2后1结果相同，那就可能可以并行，否则不能并行。

didi_dudu · 发表于 Post on 2019-8-2 15:22:41

niobium 发表于 2019-8-2 12:54
你又说两个程序是独立的，又说2需要1的数据。

像你说的想同时运行1和2，有一个简单的判断方法，如果先1 ...

我说的是2需要的是上一轮1读取的数据。呃我可以换一个问题，就是两个子程序program1和program2完全独立，那能不能在fortran中让1和2并行呢？
因为我有两个程序，1个是从硬盘读取数据A，对cpu要求不高，2是对内存里已经有的数据B做运算。1和2并行是不是效率就会高一点？

didi_dudu · 发表于 Post on 2019-8-2 15:24:04

万里云发表于 2019-8-2 11:49
这个没法并行。不读完数据怎么处理？

有个do循环嘛，比如运行到第n轮的时候，第n轮的program2会用到第n-1轮的program1从硬盘读取进内存的数据

niobium · 发表于 Post on 2019-8-2 15:46:04

didi_dudu 发表于 2019-8-2 15:22
我说的是2需要的是上一轮1读取的数据。呃我可以换一个问题，就是两个子程序program1和program2完全独立， ...

可以，最简单的方式就是用openmp的sections

sobereva · 发表于 Post on 2019-8-3 00:28:00

#pragma omp sections
{
{ Work1(); }
#pragma omp section
{ Work2();
Work3(); }
#pragma omp section
{ Work4(); }
}

This code indicates that any of the tasks Work1, Work2 + Work3 and Work4 may
run in parallel, but that Work2 and Work3 must be run in sequence. Each work
is done exactly once.

didi_dudu · 发表于 Post on 2019-8-9 16:18:40

本帖最后由 didi_dudu 于 2019-8-10 10:47 编辑

sobereva 发表于 2019-8-3 00:28
This code indicates that any of the tasks Work1, Work2 + Work3 and Work4 may
run in parallel, b ...

发现问题出在我没在程序最前面加上    include 'omp_lib.h'。用ifort编译时，不加是没问题的，用pgf90, 不加include，跑sections就不行了。但是omp parallel do 还是可以执行。。。。
============================
!$omp sections
!$omp section
            do i = 1, 200
               write(6,*)'i', i
            enddo
!$omp section
            do j = 1, 200
               write(6,*)'j', j
            enddo
!$omp end sections
社长，我用fortran尝试这样写了一段，我理解的这样输出来应该是i和j混合排列的，然而实际输出来的结果还是先打完i再打的j，能不能帮我看看这段代码有什么问题，谢谢！

didi_dudu · 发表于 Post on 2019-12-10 16:35:48

sobereva 发表于 2019-8-3 00:28
This code indicates that any of the tasks Work1, Work2 + Work3 and Work4 may
run in parallel, b ...

社长请问一下，我想实现以下功能，让test1和test2用section并行起来，然后同时让test1和test2内部也可以进行多核并行。section里面有什么相关设置可以实现这种功能么？我目前这种代码结构上，发现test1和test2内部的子程序，即使是利用了openmp parallel do，也只能单核运算了。
====================
!$omp parallel sections
!$omp section
call test1
!$omp section
call test2
!$omp end parallel sections

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[并行运算] fortran调用子程序的并行调用问题求解答

评分 Rate

评分 Rate