c - 将 Fortran 数组传递给 C

Question

我在将 Fortran 数组传递给 C 程序时遇到了很多麻烦。从我从以前的帖子中收集到的信息是包含界面。这摆脱了我的一些问题。但是，我似乎无法弄清楚如何正确传递这些数组或在 C 中正确访问它们的值。

program f_call_c
  implicit none


  interface
     subroutine cfun(x,len) bind( c )
       use,intrinsic :: iso_c_binding
       implicit none
       integer( c_int) :: len
       real(c_double) :: x(0:3,0:len)
     end subroutine   cfun

     subroutine vec(r,len) bind(c)
       use,intrinsic :: iso_c_binding
       implicit none
       integer(c_int) :: len
       real(c_double) :: r(0:len)
     end subroutine vec


  end interface


  double precision, allocatable :: x(:,:),r(:)
  integer :: len,row,i,j

  len = 7
  allocate(x(0:3,0:len))
  allocate(r(0:len))

  do i =0,len
     r(i) = i
  enddo


  do i = 0,len
     do j = 0,3
        x(j,i) = j+i
     enddo
  enddo

  call vec(r,%val(len) )

  row = 3
  call cfun(x,%val(len))
end program f_call_c


#include <stdio.h>
void cfun(double **x,const int len)
{
  printf("%d\n", len);
  printf("This is in C function cfun...\n");


  for(int i=0; i<len; i++)
    {
      printf(" %d\n  %d\n  %d\n", x[0][i]);

    }
}

void vec(  double *r[],const int len )
{
  printf("This is in C function vec...\n");

  printf("%d\n", len);

  for(int i=0; i<len; i++)
    {
       printf(" %d\n", r[i]);

    }

}

目前，输出是

    Fortran calling C, passing
 r at           0 is   0.0000000000000000     
 r at           1 is   1.0000000000000000     
 r at           2 is   2.0000000000000000     
 r at           3 is   3.0000000000000000     
 r at           4 is   4.0000000000000000     
 r at           5 is   5.0000000000000000     
 r at           6 is   6.0000000000000000     
 r at           7 is   7.0000000000000000     
This is in C function vec...
7
 0
 0
 0
 0
 0
 0
 0
7
This is in C function cfun...
Segmentation fault (core dumped)

有什么建议吗？

Answer 1

你不能在 c 中接收一个 fortran 数组作为 double **，它应该是 double *，所以试试这个

#include <stdio.h>

void cfun(double *x, const int len)
{
    printf("%d\n", len);
    printf("This is in C function cfun...\n");

    for (int i = 0 ; i < len ; i++)
    {
        printf(" %d\n  %d\n  %d\n", x[i]);
    }
}

事实上，如果你有一个 c double ** 指针数组，你应该将这些数组连接成一个数组以将其传递给 fortran，例如参见如何使用 Lapack在 c.

原因是在 fortran 中，二维数组是连续存储的，而在 c 中，double ** 是一个指针数组，因此值不是连续存储的。

请注意，当打印值时，您将打印错误的值，因为您没有为 double 使用适当的格式说明符，您应该修复 printf 行以使其成为看起来像这样

printf(" %f\n  %f\n  %f\n", x[i]);