[英]How to use Fortran interface to call a C function which contains user defined type
实际上,我想从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,有人知道如何为其创建接口吗? 提前致谢
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
选项之类的东西中幸存下来。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.