Compiler construction 传递包含不同编译器（PGI和Intel）之间可分配数组的fortran派生类型
我们有一个项目，它发展了Nvidia GPU和英特尔至强Phi。主机代码和GPU代码用Fortran编写，并由pgfortran编译。为了将部分工作转移到Phi，我们必须创建一个由ifort编译的共享库（静态链接无法工作），并从代码的pgfortran部分调用共享子例程。通过这样做，我们可以将阵列从pgfortran部分的代码卸载到可以与Xeon Phi通信的intel fortran共享库 现在，我尝试将一个派生类型（包含可分配数组）从代码的pgfortran部分传递到ifort共享库。看来有些问题 下面是一个简单的示例（此处没有Xeon Phi卸载指令）： caller.f90：Compiler construction 传递包含不同编译器（PGI和Intel）之间可分配数组的fortran派生类型,compiler-construction,fortran,intel,pgi,Compiler Construction,Fortran,Intel,Pgi,我们有一个项目，它发展了Nvidia GPU和英特尔至强Phi。主机代码和GPU代码用Fortran编写，并由pgfortran编译。为了将部分工作转移到Phi，我们必须创建一个由ifort编译的共享库（静态链接无法工作），并从代码的pgfortran部分调用共享子例程。通过这样做，我们可以将阵列从pgfortran部分的代码卸载到可以与Xeon Phi通信的intel fortran共享库 现在，我尝试将一个派生类型（包含可分配数组）从代码的pgfortran部分传递到ifort共享库。看来有
program caller
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell) :: cl(2)
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
call offload(cl(1))
print *, cl(1)%c
end program caller
program caller
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell),target :: cl(2)
type(cell_c) :: cl_c
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
cl_c%a = c_loc(cl(1)%a)
cl_c%b = c_loc(cl(1)%b)
cl_c%c = c_loc(cl(1)%c)
cl_c%na = size(cl(1)%a)
cl_c%nb = size(cl(1)%b)
cl_c%nc = size(cl(1)%c)
cl_c%id = cl(1)%id
call offload(cl_c)
print *, cl(1)%c
end program caller
subroutine offload(cl)
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
type(cell) :: cl
integer :: n
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
subroutine offload(cl_c)
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, pointer :: a(:)
real, pointer :: b(:)
real, pointer :: c(:)
end type cell
type(cell) :: cl
type(cell_C) :: cl_C
integer :: n
cl%id = cl_C%id
call c_f_pointer(cl_C%a, cl%a, [cl_c%na])
call c_f_pointer(cl_C%b, cl%b, [cl_c%nb])
call c_f_pointer(cl_C%c, cl%c, [cl_c%nc])
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
run: caller.o libcalled.so
pgfortran -L. caller.o -lcalled -o $@
caller.o: caller.f90
pgfortran -c caller.f90
libcalled.so: called.f90
ifort -shared -fPIC $^ -o $@
cl%a（1:10）（1:10）cl（1）%acl（1）%bcl%a（1:10）print*，cl%b（1:10）cl%b（1:10）谢谢大家! 编译器的ABI可能不同。您不应该直接传递结构，而是在子例程中构建它们并使用指针，这些指针应该作为
类型（c_ptr）sequencebind（C）subroutine offload(cl)
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
type(cell) :: cl
integer :: n
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
subroutine offload(cl_c)
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, pointer :: a(:)
real, pointer :: b(:)
real, pointer :: c(:)
end type cell
type(cell) :: cl
type(cell_C) :: cl_C
integer :: n
cl%id = cl_C%id
call c_f_pointer(cl_C%a, cl%a, [cl_c%na])
call c_f_pointer(cl_C%b, cl%b, [cl_c%nb])
call c_f_pointer(cl_C%c, cl%c, [cl_c%nc])
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
program caller
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell) :: cl(2)
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
call offload(cl(1))
print *, cl(1)%c
end program caller
program caller
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell),target :: cl(2)
type(cell_c) :: cl_c
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
cl_c%a = c_loc(cl(1)%a)
cl_c%b = c_loc(cl(1)%b)
cl_c%c = c_loc(cl(1)%c)
cl_c%na = size(cl(1)%a)
cl_c%nb = size(cl(1)%b)
cl_c%nc = size(cl(1)%c)
cl_c%id = cl(1)%id
call offload(cl_c)
print *, cl(1)%c
end program caller
>gfortran called.f90 -c -o called.o
>ifort caller.f90 -c -o caller.o
>ifort -o a.out called.o caller.o -lgfortran
>./a.out
11.0000000 12.0000000 13.0000000 14.0000000 15.0000000 16.0000000 17.0000000 18.0000000 19.0000000 20.0000000
11.0000000 21.0000000 31.0000000 41.0000000 51.0000000 61.0000000 71.0000000 81.0000000 91.0000000 101.000000
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
c_intc_float您还可以重载
单元格单元格_C类型（c_ptr）sequencebind（C）subroutine offload(cl)
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
type(cell) :: cl
integer :: n
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
subroutine offload(cl_c)
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, pointer :: a(:)
real, pointer :: b(:)
real, pointer :: c(:)
end type cell
type(cell) :: cl
type(cell_C) :: cl_C
integer :: n
cl%id = cl_C%id
call c_f_pointer(cl_C%a, cl%a, [cl_c%na])
call c_f_pointer(cl_C%b, cl%b, [cl_c%nb])
call c_f_pointer(cl_C%c, cl%c, [cl_c%nc])
print *, cl%a(1:10)
print *, cl%b(1:10)
end subroutine offload
program caller
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell) :: cl(2)
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
call offload(cl(1))
print *, cl(1)%c
end program caller
program caller
use iso_c_binding
type, bind(C) :: cell_C
integer :: id
integer :: na, nb, nc
type(c_ptr) :: a,b,c
end type cell_C
type cell
integer :: id
real, allocatable :: a(:)
real, allocatable :: b(:)
real, allocatable :: c(:)
end type cell
integer :: n,i,j
type(cell),target :: cl(2)
type(cell_c) :: cl_c
n=10
do i=1,2
allocate(cl(i)%a(n))
allocate(cl(i)%b(n))
allocate(cl(i)%c(n))
end do
do j=1, 2
do i=1, n
cl(j)%a(i)=10*j+i
cl(j)%b(i)=10*i+j
end do
end do
cl_c%a = c_loc(cl(1)%a)
cl_c%b = c_loc(cl(1)%b)
cl_c%c = c_loc(cl(1)%c)
cl_c%na = size(cl(1)%a)
cl_c%nb = size(cl(1)%b)
cl_c%nc = size(cl(1)%c)
cl_c%id = cl(1)%id
call offload(cl_c)
print *, cl(1)%c
end program caller
>gfortran called.f90 -c -o called.o
>ifort caller.f90 -c -o caller.o
>ifort -o a.out called.o caller.o -lgfortran
>./a.out
11.0000000 12.0000000 13.0000000 14.0000000 15.0000000 16.0000000 17.0000000 18.0000000 19.0000000 20.0000000
11.0000000 21.0000000 31.0000000 41.0000000 51.0000000 61.0000000 71.0000000 81.0000000 91.0000000 101.000000
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
c_intc_float您还可以重载
单元格单元格_C