检查 Fortran 调用 C 函数中的 ISO_C_BINDING 用法
Check for ISO_C_BINDING usage in Fortran calling C function
我在ISO_C_BINDING
上问了几个问题,我做了一个例程。由于我在处理 Fortran POINTER 时总是犯了很多错误,所以我想确保没有错误或一些不推荐的奇怪点。
(Fortran 部分)
program test
use iso_c_binding
implicit none
interface
subroutine get_value_array(in, num) bind(C, name='get_value_Array')
use iso_c_binding
implicit none
type(C_PTR), intent(inout) :: in
integer(C_INT), value, intent(in) :: num
end subroutine
end interface
real(C_DOUBLE), allocatable, target :: array(:)
real(C_DOUBLE), pointer :: array_fptr(:)
type(C_PTR) :: array_cptr
integer :: array_len
allocate(array(12))
array_len = size(array,1)
array_cptr = C_LOC(array)
call get_value_array (array_cptr, array_len)
call C_F_POINTER(array_cptr, array_fptr, [array_len])
print *, 'array_fptr'
print *, array_fptr
print *, 'array'
print *, array
end program
(C部分)
void get_value_Array(double **in, int num) {
int i;
for (i = 0; i < num; i++) {
(*in)[i] = i+1;
}
}
(输出)
array_fptr
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000 6.00000000000000
7.00000000000000 8.00000000000000 9.00000000000000
10.0000000000000 11.0000000000000 12.0000000000000
array
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000 6.00000000000000
7.00000000000000 8.00000000000000 9.00000000000000
10.0000000000000 11.0000000000000 12.0000000000000
为了确定,流程是
array_cptr
与具有 C_LOC
功能的 array
关联。 (也许 'associate' 在这里不是一个合适的术语,因为它是 Fortran 指针相关的术语。我应该怎么称呼它?)
array
在我调用 get_value_array
时被修改,因为子例程(或 C 函数)修改了 array_cptr
指向的位置。
当我调用 C_F_POINTER
时,array_fptr
通过 array_cptr
关联到 array
最后 array_cptr
和 array_fptr
是彼此不同的对象 (?),同时指向同一个目标 array
,对吗?
问题中的陈述含糊不清。准确地说。
用Fortran语句array_cptr = C_LOC(array)
,array_cptr
被赋值为array
.
的C地址
调用get_value_array
时,array
的值最终被C函数修改。 Fortran 代码还允许子例程(或 C 函数)修改存储在 array_cptr
中的 C 地址 - 伪参数声明为 INTENT(INOUT) 并且没有 VALUE 属性,并且与C 函数通过引用获取指向数组的指针(即 C 形式参数被声明为指向双精度指针的指针)。但是C函数不会对array_cptr
中存储的C地址做这样的修改。对于显示的代码,存在不必要的间接级别。
Fortran 代码将 array_fptr
与 array_cptr
中的 C 地址引用的事物相关联。因为array_cptr
中的C地址没有被C函数改变,所以array_fptr继续与array
关联。在这种情况下执行 C_F_POINTER 不会改变程序的状态。
array_cptr
与 array_fptr
不同 - 但它们具有相同的目标(array_cptr 是一个对象,array_fptr 与对象相关联 - 中的指针与 C 不同,Fortran 本身不是对象。
我在ISO_C_BINDING
上问了几个问题,我做了一个例程。由于我在处理 Fortran POINTER 时总是犯了很多错误,所以我想确保没有错误或一些不推荐的奇怪点。
(Fortran 部分)
program test
use iso_c_binding
implicit none
interface
subroutine get_value_array(in, num) bind(C, name='get_value_Array')
use iso_c_binding
implicit none
type(C_PTR), intent(inout) :: in
integer(C_INT), value, intent(in) :: num
end subroutine
end interface
real(C_DOUBLE), allocatable, target :: array(:)
real(C_DOUBLE), pointer :: array_fptr(:)
type(C_PTR) :: array_cptr
integer :: array_len
allocate(array(12))
array_len = size(array,1)
array_cptr = C_LOC(array)
call get_value_array (array_cptr, array_len)
call C_F_POINTER(array_cptr, array_fptr, [array_len])
print *, 'array_fptr'
print *, array_fptr
print *, 'array'
print *, array
end program
(C部分)
void get_value_Array(double **in, int num) {
int i;
for (i = 0; i < num; i++) {
(*in)[i] = i+1;
}
}
(输出)
array_fptr
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000 6.00000000000000
7.00000000000000 8.00000000000000 9.00000000000000
10.0000000000000 11.0000000000000 12.0000000000000
array
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000 6.00000000000000
7.00000000000000 8.00000000000000 9.00000000000000
10.0000000000000 11.0000000000000 12.0000000000000
为了确定,流程是
array_cptr
与具有C_LOC
功能的array
关联。 (也许 'associate' 在这里不是一个合适的术语,因为它是 Fortran 指针相关的术语。我应该怎么称呼它?)array
在我调用get_value_array
时被修改,因为子例程(或 C 函数)修改了array_cptr
指向的位置。
当我调用 array_fptr
通过array_cptr
关联到array
C_F_POINTER
时,最后 array_cptr
和 array_fptr
是彼此不同的对象 (?),同时指向同一个目标 array
,对吗?
问题中的陈述含糊不清。准确地说。
用Fortran语句
的C地址array_cptr = C_LOC(array)
,array_cptr
被赋值为array
.调用
get_value_array
时,array
的值最终被C函数修改。 Fortran 代码还允许子例程(或 C 函数)修改存储在array_cptr
中的 C 地址 - 伪参数声明为 INTENT(INOUT) 并且没有 VALUE 属性,并且与C 函数通过引用获取指向数组的指针(即 C 形式参数被声明为指向双精度指针的指针)。但是C函数不会对array_cptr
中存储的C地址做这样的修改。对于显示的代码,存在不必要的间接级别。Fortran 代码将
array_fptr
与array_cptr
中的 C 地址引用的事物相关联。因为array_cptr
中的C地址没有被C函数改变,所以array_fptr继续与array
关联。在这种情况下执行 C_F_POINTER 不会改变程序的状态。
array_cptr
与 array_fptr
不同 - 但它们具有相同的目标(array_cptr 是一个对象,array_fptr 与对象相关联 - 中的指针与 C 不同,Fortran 本身不是对象。