在 C 中修改 Fortran 数组

Question

在 Fortran 中，我有一个数组指定为：

dimension insv(*)

现在我想从 C 中编辑这个数组。 我设法通过使用 iso_c_binding 从 C 中读取数组：

  SUBROUTINE userroutine(insv)
  IMPLICIT NONE
  INTERFACE
  SUBROUTINE MODIFYARRAY(insv) BIND(C)
  USE, INTRINSIC::ISO_C_BINDING
     REAL(C_DOUBLE),DIMENSION(*), INTENT(INOUT) :: insv
  END SUBROUTINE MODIFYARRAY
  END INTERFACE

  dimension insv(*)

  WRITE(insv(0))
  call modifyarray(insv)
  WRITE(insv(0))
  END

我的 C 函数：

__declspec(dllexport) void modifyarray(double * insv)

但是在 C 中修改这个数组中的值要么崩溃，要么根本没有改变。

__declspec(dllexport) void modifyarray(double * insv)
{
    insv[0] = 1234.00;
}

我不确定我应该对 iso_c_binding 或 C 函数进行什么更改才能使其工作。

更多上下文：fortran 数组类型无法更改，因为维度变量是用户例程的参数，而不是在我的控制范围内，可以这么说。 进入fortran函数时数组已经被填满，但需要在C中修改。

Answer 1

假设：

• insv 是 real(C_DOUBLE) 的数组
• 你想写insv（而不是现在看来由insv定义的文件单元）
• 调用程序将insv定义为一个包含10个组件的数组

一个可能的工作测试用例如下：

Fortran：

SUBROUTINE userroutine(insv)
    USE, INTRINSIC::ISO_C_BINDING
    IMPLICIT NONE
    INTERFACE
        SUBROUTINE MODIFYARRAY(insv) BIND(C)
            USE, INTRINSIC::ISO_C_BINDING
            REAL(C_DOUBLE),DIMENSION(*), INTENT(INOUT) :: insv
        END SUBROUTINE MODIFYARRAY
    END INTERFACE

    real(C_DOUBLE), dimension(*) :: insv

    WRITE(*,*) insv(1)
    call modifyarray(insv)
    WRITE(*,*) insv(1)
END

PROGRAM MAIN
    USE, INTRINSIC::ISO_C_BINDING
    real(C_DOUBLE) :: insv(10)
    insv = 2.d0
    call userroutine(insv)
END

C:

void modifyarray(double * insv)
{   
    insv[0] = 1234.00;
}