关于openmp并行域数组传递的问题(fortran代码)

您好,我是一名力学研究生,目前在使用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