c - 如何使用 Fortran 接口调用包含用户定义类型的 C 函数

Question

事实上，我想从 fortran 中调用岩浆。所以我添加了magma.lib并创建了一个接口来使用magma的C函数：

Interface
      Integer function magma_dpotrf(uplo, n, a, lda, info) BIND (C, NAME="magma_dpotrf")
        use iso_c_binding
        Implicit none
        !character (c_char), value :: uplo????
        integer (c_int), value ::n
        real (c_double) ::a(*)
        integer (c_int), value ::lda
        integer (c_int)::info
      end function
   end Interface

但是参数 uplo 是 C 代码中的用户定义类型（magma_uplo_t uplo）：

typedef enum {
    MagmaUpper         = 121,
    MagmaLower         = 122,
    MagmaUpperLower    = 123,
    MagmaFull          = 123,  /* lascl, laset */
    MagmaHessenberg    = 124   /* lascl */
} magma_uplo_t;

magma_int_t
magma_dpotrf(
    magma_uplo_t uplo, magma_int_t n,
    double *A, magma_int_t lda,
    magma_int_t *info);

magma_int_t = int，有谁知道如何为它创建接口？提前致谢

Answer 1

magma_uplo_t是一个枚举。Fortran 2003 中对它们有一些支持，但您可以完全放心地假设它是一个integer(c_int)可以将值从 121 到 124 的值。在您的情况下，它是按值传递的。

integer(c_int), value :: uplo

您实际上可以使用 Fortran 2003 枚举创建常量

 enum, bind( C )
    enumerator :: MagmaUpper         = 121, &
                  MagmaLower         = 122, &
                  MagmaUpperLower    = 123, &
                  MagmaFull          = 123, &
                  MagmaHessenberg    = 124
  end enum

但是变量然后你也可以尝试integer(kind=kind(MagmaUpper))完全安全。这将在 GCC 的--short-enums选项之类的情况下继续存在。