我在将二维数组从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/

10-13 07:10
查看更多