使用unsafe_store!向Fortran写入结构体数组的内存对齐问题

问题根源 #

出现值不匹配和大小错误的核心原因有两个：

1. 类型大小不兼容：Fortran默认real为单精度（4字节，对应Julia的Float32），而Julia结构体用了Float64（8字节）；Fortran默认logical为4字节，Julia的Bool为1字节，直接导致内存布局错位，读取到垃圾值。
2. Fortran可分配数组的特殊性：Fortran的allocatable数组并非单纯的内存指针，它包含额外的元数据（如数组大小、分配状态等），直接写入裸指针会让Fortran无法正确解析数组信息。

解决方案 #

1. 统一结构体类型与内存布局 #

在Fortran中使用iso_c_binding模块定义与C/Julia兼容的类型，确保类型大小和布局完全匹配；Julia结构体保持对应类型即可，默认布局与C兼容，无需额外调整。

2. 替换可分配数组为指针类型 #

将Fortran中的allocatable数组改为pointer类型，配合单独的变量存储数组大小，这样可以直接接收Julia传入的裸指针并正确解析。

修改后的代码 #

Fortran模块（testmodule.f90） #

module testmodule

  use, intrinsic :: iso_c_binding, only: c_double, c_bool

  ! 用C绑定类型定义结构体，确保与Julia类型匹配
  type :: JulParam
    real(c_double) :: &
      f1, &
      f2
    logical(c_bool) :: &
      testbool
  end type JulParam

  ! 用pointer类型替代allocatable，直接接收外部指针
  type(JulParam), pointer, dimension(:) :: test_param
  ! 单独存储数组大小
  integer :: test_param_size

  contains

  module subroutine test_func()
    print *, "test size", test_param_size
    print *, "test val1", test_param(1)%f1
    print *, "test bool1", test_param(1)%testbool
  end subroutine test_func

end module testmodule

Julia代码（test.jl） #

using Base.Libc.Libdl
using Base.Libc

! 结构体类型与Fortran的c_double、c_bool完全对应
mutable struct jul_param
    f1::Float64
    f2::Float64
    testbool::Bool
end

try
    arr_size = 2
    testparam = jul_param(1.0, 2.0, true)
    ! 分配内存存储结构体数组
    ptr_jul_struct_array = convert(Ptr{jul_param}, calloc(arr_size, sizeof(jul_param)))
    unsafe_store!(ptr_jul_struct_array, testparam, 1)
    unsafe_store!(ptr_jul_struct_array, testparam, 2)

    testmodule_bin = Libdl.dlopen("testmodule.so")

    ! 设置Fortran的pointer数组指针
    test_param_sym = Libdl.dlsym(testmodule_bin, Symbol("__testmodule_MOD_test_param"))
    ptr_fortran_ptr_array = convert(Ptr{Ptr{jul_param}}, test_param_sym)
    unsafe_store!(ptr_fortran_ptr_array, ptr_jul_struct_array)

    ! 设置数组大小变量
    test_size_sym = Libdl.dlsym(testmodule_bin, Symbol("__testmodule_MOD_test_param_size"))
    ptr_fortran_size = convert(Ptr{Cint}, test_size_sym)
    unsafe_store!(ptr_fortran_size, arr_size)

    ! 调用Fortran子程序
    function_sym = Libdl.dlsym(testmodule_bin, Symbol("__testmodule_MOD_test_func"))
    ccall(function_sym, Cvoid, ())

    ! 清理资源
    free(ptr_jul_struct_array)
    Libdl.dlclose(testmodule_bin)
catch e
    println("Error: ", e)
    rethrow(e)
end

编译运行命令 #

gfortran -fpic -c testmodule.f90 && gfortran -shared -o testmodule.so testmodule.o && julia test.jl

运行后将正确输出：

test size           2
 test val1   1.0000000000000000
 test bool1 T

内容的提问来源于stack exchange，提问作者Yes