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

Question

您好，我是一名力学研究生，目前在使用vs的openmp将我写的串行代码并行化，但是有一个很关键的问题遇到了麻烦，之前的串行代码，函数之间的数组传递我使用的是动态数组（指针），也就是假定大小数组(assumed-size)，使用openmp需要将数组私有化，但是私有变量的设置好像不能使用assumed-size类型的数组，想问一下，是否openmp私有变量数组的形式只能是大小确定的数组，即外层函数传递给内层函数的数组的维数的大小必须是确定的，不能使用指针这种动态数组。

Answer 1

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