c - 如何在 C 中访问(动态分配的)Fortran 数组

Question

我的主要问题是为什么数组会做如此奇怪的事情，以及是否有任何方法可以以“干净”的方式执行以下操作。

我目前有一个 C 程序 foo.c 通过 dlopen/dlsym 与 Fortran 程序 bar.f90 接口(interface)，大致类似于代码下面:

foo.c:

#include <dlfcn.h>
#include <stdio.h>

int main()
{
int i, k = 4;
double arr[k];
char * e;

void * bar = dlopen("Code/Test/bar.so", RTLD_NOW | RTLD_LOCAL);

void (*allocArray)(int*);
*(void **)(&allocArray) = dlsym(bar, "__bar_MOD_allocarray");
void (*fillArray)(double*);
*(void **)(&fillArray) = dlsym(bar, "__bar_MOD_fillarray");
void (*printArray)(void);
*(void **)(&printArray) = dlsym(bar, "__bar_MOD_printarray");
double *a = (double*)dlsym(bar, "__bar_MOD_a");

for(i = 0; i < k; i++)
    arr[i] = i * 3.14;

(*allocArray)(&k);
(*fillArray)(arr);
(*printArray)();
for(i = 0; i < 4; i++)
    printf("%f ", a[i]);
printf("\n");

return 0;
}

bar.f90:

module bar

integer, parameter :: pa = selected_real_kind(15, 307)
real(pa), dimension(:), allocatable :: a
integer :: as

contains

subroutine allocArray(asize)
    integer, intent(in) :: asize

    as = asize
    allocate(a(asize))

    return
end subroutine

subroutine fillArray(values)
    real(pa), dimension(as), intent(in) :: values

    a = values
    return
end subroutine

subroutine printArray()
    write(*,*) a
    return
end subroutine

end module

运行主要 yield

0.0000000000000000        3.1400000000000001        6.2800000000000002        9.4199999999999999     
0.000000 -nan 0.000000 0.000000 

这表明 Fortran 正确地分配了数组，甚至正确地存储了给定的值，但它们不能再通过 dlsym 访问(处理该数据会导致段错误)。我也对固定大小的数组进行了尝试 - 结果保持不变。

有人知道这种行为的原因吗？我个人希望事情要么双向解决，要么根本不解决 - 这个“Fortran 接受 C 数组但反之亦然”让我想知道我在以这种方式从 C 访问数组时是否犯了一些基本错误。

另一个(甚至更重要的)问题是，如何“以正确的方式”进行数组访问。目前我什至不确定坚持使用“Fortran as .so”接口(interface)是否是一个好方法——我认为在这种情况下也可以尝试混合编程。尽管如此，数组问题仍然存在——我读到这可以使用 ISO C 绑定(bind)以某种方式解决，但我还不知道如何解决(我还没有用 Fortran 工作过很多，但是，尤其是没有说绑定(bind)) , 因此，我们将不胜感激对此问题的帮助。

编辑:

好的，所以我进一步阅读了 ISO C 绑定(bind)并找到了一个非常有用的方法 here .使用 C_LOC 我可以获得指向我的 Fortran 结构的 C 指针。不幸的是，指向数组的指针似乎是指向指针的指针，需要在 C 代码中解除引用，然后才能将它们视为 C 数组 - 或类似的东西。

编辑:

按照 Vladimir F 指出的方式使用 C 绑定(bind)使我的程序现在可以运行，至少在大多数情况下是这样。 C 文件和 Fortran 文件现在链接在一起，所以我可以避免 libdl 接口(interface)，至少对于 Fortran 部分 - 我仍然需要加载一个动态 C 库，获取指向其中一个符号的函数指针并传递它作为指向 Fortran 的函数指针，稍后会调用该函数作为其计算的一部分。正如所说的函数需要 double*s [arrays]，我无法使用 C_LOC 传递我的 Fortran 数组，奇怪的是 - C_LOC(array) 和 C_LOC(array(1)) 将正确的指针传递回 C 函数。 array(1) 成功了。可悲的是，这不是执行此操作的“最干净”的方法。如果有人给我提示如何使用 C_LOC 函数执行此操作，那就太好了。尽管如此，我还是接受了 Vladimir F 的回答，因为我认为这是更安全的解决方案。

Answer 1

许多 Fortran 编译器在内部使用称为数组描述符的东西 - 保存数组形状的结构(即每个维度的大小和范围以及指向实际数据的指针)。它允许实现假设形状数组参数、数组指针和可分配数组等工作。您通过 __bar_MOD_a 符号访问的是可分配数组的描述符，而不是其数据。

数组描述符是特定于编译器的，依赖于特定描述符格式的代码不可移植。示例描述符:

• GNU Fortran
• Intel Fortran

请注意，即使是那些也特定于这些编译器的某些版本。例如，Intel 声明他们当前的描述符格式与 Intel Fortran 7.0 中使用的格式不兼容。

如果您查看这两个描述符，您会发现它们非常相似，并且第一个元素是指向数组数据的指针。因此，您可以使用 double ** 而不是 double * 轻松读取数据:

double **a_descr = (double**)dlsym(bar, "__bar_MOD_a");
...
for(i = 0; i < 4; i++)
    printf("%f ", (*a_descr)[i]);

再一次，这不是可移植的，因为这些描述符的格式将来可能会改变(尽管我怀疑数据指针会被移动到描述符开头以外的其他地方)。有一份规范草案试图统一所有描述符格式，但尚不清楚不同编译器供应商将如何以及何时采用它。

编辑: 下面是如何使用访问器函数，该函数使用 ISO_C_BINDING 模块中的 C_LOC() 可移植地获取指向可分配数组:

Fortran 代码:

module bar
  use iso_c_binding
  ...
  ! Note that the array should be a pointer target
  real(pa), dimension(:), allocatable, target :: a
  ...
contains
  ...

  function getArrayPtr() result(cptr)
    type(c_ptr) :: cptr

    cptr = c_loc(a)
  end function

end module

C 代码:

...
void * (*getArrayPtr)(void);
*(void **)(&getArrayPtr) = dlsym(bar, "__bar_MOD_getarrayptr");
...
double *a = (*getArrayPtr)();
for(i = 0; i < 4; i++)
    printf("%f ", a[i]);
...

结果:

$ ./prog.x
   0.0000000000000000        3.1400000000000001        6.2800000000000002
 9.4199999999999999
0.000000 3.140000 6.280000 9.420000