我在OpenACC内核中遇到了一个数组问题。以下是演示代码:
module mpoint
type point
real :: x, y, z
real :: tmp
real :: v(3)
end type point
end module mpoint
program main
use mpoint
implicit none
integer, parameter :: n = 3
real, allocatable :: array(:)
!--------------------------------
allocate(array(n))
array(:) = 1.0
call vecadd
contains
subroutine vecadd()
integer :: i
type(point) :: A
real :: w(3)
real :: scalar
A%v(:) = 0.0
A%v(1) = 1.0
w(:) = 0.0
w(1) = 1.0
scalar = 1.0
write(*,*) 'host: v(1) = ', A%v(1)
write(*,*) 'host: w(1) = ', w(1)
write(*,*) 'host: scalar = ', scalar
!$acc parallel loop firstprivate(A,w)
do i = 1, n
write(*,*) 'device: v(1) = ', A%v(1)
write(*,*) 'device: w(1) = ', w(1)
write(*,*) 'device: scalar = ', scalar
enddo
end subroutine vecadd
end program main当我用nvfortran -acc -Minfo=accel test.f90编译它并运行时,它显示在设备上数组中的值是0.0,而不是我在主机端设置的正确值1.0。这只发生在数组中:标量,如示例所示,具有正确的值。
我想知道这是nvfortran的限制,还是当前的OpenACC标准?
发布于 2020-10-14 04:20:56
看起来像是编译器的问题,我们没有初始化主机上的小数组。在浏览我们现有的错误报告时,我发现了一个几乎相同的错误报告,只是用C而不是Fortran,巧合的是今天早上在我们的开发编译器中修复了这个错误。不幸的是,它似乎并没有解决你的问题。我给负责这个问题的编译器工程师发了一封信,询问他是否可以看一看。
最坏的情况是,如果问题相似但不相关,我将打开一个新的问题报告,并使用跟踪号更新这篇文章。
https://stackoverflow.com/questions/64341827
复制相似问题