我很难将Fortard的2d数组传递给C函数。但是,在所有支持之后,以下代码100%正常运行。
以下是我的C函数:
#include <stdio.h>
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代码,应该按照以下步骤操作:
谢谢, 阿纳斯
答案 0 :(得分:-1)
有两种方式可以存储像你一样的数组--ROW-MAJOR和COLUMN-MAJOR。 C使用row-major,fortran使用column-major。
您必须按照在fortran中创建的顺序访问C中的数组元素。
在编写代码块之前,您可能最好先了解一下代码块。为清楚起见,请参见此 http://en.wikipedia.org/wiki/Row-major_order