2

我正在将 Fortran 代码与 C dll 绑定，并且我希望有一个可与 C 互操作的 Fortran 数组。我目前有以下子例程将 Fortran 数组与 C double* 绑定：

SUBROUTINE Pack_Inputs( Input , In_X )
  TYPE( InputType ) ,               INTENT(INOUT)            :: Input
  REAL(KIND=C_DOUBLE) , ALLOCATABLE , TARGET , INTENT(INOUT) :: In_X(:)

  IF ( .NOT. ALLOCATED(In_X)     ) ALLOCATE( In_X    (Input%Xlen)     )

  DO i = 1,Input%C_obj%Xlen       
     In_X(i) = Input%X(i)
  END DO
  Input%C_obj%X = C_LOC(In_X)    
END SUBROUTINE Pack_Inputs

但是，我不喜欢当前代码的是我不断分配内存，并且在输入 C dll 时必须解包数组（部分原因是我不愿意使用SAVEon 属性In_X(:)）。我宁愿 In_X在 Fortran 派生类型中声明一次。这导致了这篇文章的动机。在这个派生类型中：

USE , INSTRINSIC :: ISO_C_BINDING

TYPE , PUBLIC  :: InputType                               
   TYPE(InputType_C)                          :: C_obj   
   REAL(KIND=C_DOUBLE) , ALLOCATABLE , TARGET :: In_X(:)   
   REAL    , DIMENSION(:) , ALLOCATABLE       :: X                    
   REAL    , DIMENSION(:) , ALLOCATABLE       :: Y                    
   REAL    , DIMENSION(:) , ALLOCATABLE       :: Z                    
   INTEGER , DIMENSION(:) , ALLOCATABLE       :: index
   INTEGER                                    :: Xlen                 
   INTEGER                                    :: Ylen                 
   INTEGER                                    :: Zlen                 
   INTEGER                                    :: indexlen
END TYPE InputType

我得到错误：

REAL(KIND=C_DOUBLE) , ALLOCATABLE  , TARGET :: In_X(:)
                                          1
    Error: Attribute at (1) is not allowed in a TYPE definition

有没有办法消除这个错误？

4

1 回答 1

1

我以前遇到过这个问题，对我有用的解决方案是将组件声明为POINTER而不是ALLOCATABLE, TARGET. 我不确定 Fortran 标准是否不支持它，或者编译器没有实现此功能。我正在使用ifort v12.0.2.137.

这对您来说是一个可以接受的解决方案吗？然后，您就可以将其用作指针目标。

TYPE , PUBLIC  :: InputType                               
   TYPE(InputType_C)                          :: C_obj   
   REAL(KIND=C_DOUBLE),DIMENSION(:),POINTER   :: In_X => NULL()
   REAL    , DIMENSION(:) , ALLOCATABLE       :: X                    
   REAL    , DIMENSION(:) , ALLOCATABLE       :: Y                    
   REAL    , DIMENSION(:) , ALLOCATABLE       :: Z                    
   INTEGER , DIMENSION(:) , ALLOCATABLE       :: index
   INTEGER                                    :: Xlen                 
   INTEGER                                    :: Ylen                 
   INTEGER                                    :: Zlen                 
   INTEGER                                    :: indexlen
END TYPE InputType

然后，您可以将In_X指针与目标数据相关联：

In_X(1:Input%C_obj%Xlen) => Input%X(1:Input%C_obj%Xlen)

请注意，这Input%X也需要有一个TARGET属性。

于 2013-06-14T19:11:10.993 回答