您好,我是一名力学研究生,目前在使用vs的openmp将我写的串行代码并行化,但是有一个很关键的问题遇到了麻烦,之前的串行代码,函数之间的数组传递我使用的是动态数组(指针),也就是假定大小数组(assumed-size),使用openmp需要将数组私有化,但是私有变量的设置好像不能使用assumed-size类型的数组,想问一下,是否openmp私有变量数组的形式只能是大小确定的数组,即外层函数传递给内层函数的数组的维数的大小必须是确定的,不能使用指针这种动态数组。
1、建议你给出示例代码,简单的能够体现你的问题的代码。
2、假定大小数组,也就是形参大小为星号,程序没法直接获得它的长度(元素个数),因此没办法私有化。
3、如果是自己的代码,建议改为假定形状数组,omp可以私有化。
4、从经验而言,一般没有私有化形参数组的(我都是用shared),你思考一下是否可以改变算法。
module m
use omp_lib
contains
subroutine sub1(a)
implicit none
integer :: a(*) !假定大小
! error: An assumed sized array cannot be specified in an OpenMP* PRIVATE
! $OMP parallel num_threads(2) private(i, a)
! $OMP end parallel
end subroutine
subroutine sub2(a)
implicit none
integer :: a(:) !假定形状
integer i
! An assumed sized array cannot be specified in an OpenMP* PRIVATE
!$OMP parallel num_threads(2) private(i, a)
a = 0
i = omp_get_thread_num() + 1
a(i) = i
print*,a ! 1 0
! 0 2
!$OMP end parallel
end subroutine
end module
program test
use m
integer,allocatable :: a(:)
allocate(a(2))
a=0
call sub2(a)
end program