我在将二维数组从Fortran传递到C时遇到问题。下面是我的C函数,它只在屏幕上显示数组元素。
#include <stdio.h>
void print2(double *arr , int *n)
{
int y = *n;
printf("\n y = %d", y);
for(int i =0; i<y; i++)
{
for (int j = 0; j < y; j++)
printf("%.6g", *((arr + i*y) + j));
printf("\n");
}
}
到目前为止,我的Fortran代码如下:
program linkFwithC
use, intrinsic :: iso_c_binding
implicit none
real, dimension(3,3)::a
a(1,1)=1
a(1,2)=2
a(1,3)=3
a(2,1)=4
a(2,2)=5
a(2,3)=6
a(3,1)=7
a(3,2)=8
a(3,3)=9
interface
subroutine print2(a,n) bind( c )
use, intrinsic :: iso_c_binding
type(c_ptr)::a
integer(C_INT)::n
end subroutine print2
end interface
call print2(c_loc(a),3)
end program linkFwithC
我链接两个文件的方式是为C函数创建一个静态库并构建.lib文件。生成.lib文件后,我将其添加到fortran项目并运行fortran项目。代码运行时没有错误,n值显示正确;但是,显示的数组值都是错误的。
请帮忙!
谢谢,
阿纳斯
最佳答案
代码中有一些问题[当前]显示。
print2的Fortran接口中的n参数没有VALUE属性,但是C函数中的相应参数是按VALUE获取的。考虑给Fortran声明增加值。
指向数组的指针也会出现同样的问题。Fortran接口传递一个没有值的指针,C函数需要一个“逐值指针”(与指向指针的指针相反)。注意,这里不需要显式地使用C_PTR—您可以使用数组的实际类型构造一个可互操作的接口。
在大多数平台上,Fortran默认实数与C double不一样-考虑使用ISO_C_BINDING中的kind常量,以确保Fortran端的实数类型与C端的实数类型匹配。
C_LOC要求其参数具有TARGET属性。将该属性添加到主程序中变量的声明中。
关于c - 将二维数组从Fortran传递到C,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/27582715/