从Fortran传递2d数组到C

我很难将Fortard中的二维数组传递给C函数。 但是,在所有支持之后,以下代码100%正常运行。

以下是我的C函数:

#include  void print2(void *p, int n) { printf("Array from C is \n"); double *dptr; dptr = (double *)p; for (int i = 0; i < n; i++) { for (int j = 0; j<n; j++) printf("%.6g \t",dptr[i*n+j]); printf("\n"); } } 

以下是我的Fortran代码:

  program linkFwithC use iso_c_binding implicit none interface subroutine my_routine(p,r) bind(c,name='print2') import :: c_ptr import :: c_int type(c_ptr), value :: p integer(c_int), value :: r end subroutine end interface integer,parameter ::n=3 real (c_double), allocatable, target :: xyz(:,:) real (c_double), target :: abc(3,3) type(c_ptr) :: cptr allocate(xyz(n,n)) cptr = c_loc(xyz(1,1)) !Inputing array valyes xyz(1,1)= 1 xyz(1,2)= 2 xyz(1,3)= 3 xyz(2,1)= 4 xyz(2,2)= 5 xyz(2,3)= 6 xyz(3,1)= 7 xyz(3,2)= 8 xyz(3,3)= 9 call my_routine(cptr,n) deallocate(xyz) pause end program linkFwithC 

代码运行正常;但是,C中的数组元素需要重新组织。

注意,为了在Visual Studio环境中链接C函数和FORTRAN代码,应该按照以下步骤操作:

  • 在静态库项目中编写C函数
  • 构建.lib文件
  • 创建FORTRAN项目并编写代码
  • 将.lib文件添加到FORTRAN项目(只需将其添加到源文件)
  • 编译并运行。

谢谢,阿纳斯

存储的数组有两种方式存储 – ROW-MAJOR和COLUMN-MAJOR。 C使用row-major,fortran使用column-major。

您必须按照在fortran中创建的顺序访问C中的数组元素。

在编写代码块之前,您可能会更好地理解它,而不是发布代码块。 为清楚起见,请参阅此页: http : //en.wikipedia.org/wiki/Row-major_order