Fortran 中动态加载动态库

在 Fortran 中动态加载动态库（.so或.dll）需借助C 互操作性（ISO_C_BINDING）调用操作系统原生 API（如dlopen/dlsym或LoadLibrary/GetProcAddress）。Fortran 标准本身不提供动态加载接口，但可通过 C 绑定桥接实现。以下分步骤说明并提供完整示例。

一、核心原理

⚠️关键点：Fortran 无法直接操作void*，必须通过 C 封装层中转。

二、完整示例（Linux 平台）

1. 动态库源码（mylib.c）

// mylib.c - 编译为 libmylib.so#include<stdio.h>// 全局变量doubleglobal_value=3.14159;// C 函数（无名称修饰）doublecompute(doublex,doubley){returnx*y+global_value;}// 辅助函数：返回变量地址（避免 Fortran 直接处理 void*）double*get_global_value_ptr(){return&global_value;}

编译动态库：

gcc -shared -fPIC -o libmylib.so mylib.c

2. C 封装层（dlwrapper.c）

// dlwrapper.c - 封装 dlopen/dlsym#define_GNU_SOURCE#include<dlfcn.h>#include<stdint.h>// 返回句柄（转换为 intptr_t 供 Fortran 使用）intptr_topen_library(constchar*path){return(intptr_t)dlopen(path,RTLD_LAZY|RTLD_GLOBAL);}// 获取函数符号void*get_function(intptr_thandle,constchar*name){returndlsym((void*)handle,name);}// 获取变量地址void*get_variable(intptr_thandle,constchar*name){returndlsym((void*)handle,name);}// 关闭库voidclose_library(intptr_thandle){dlclose((void*)handle);}

编译为静态库供 Fortran 链接：

gcc -c -fPIC dlwrapper.c -o dlwrapper.o ar rcs libdlwrapper.a dlwrapper.o

3. Fortran 主程序（main.f90）

program dynamic_load_example use, intrinsic :: iso_c_binding implicit none ! 接口声明 interface function open_library(path) bind(c, name='open_library') import :: c_char, c_intptr_t character(kind=c_char), intent(in) :: path(*) integer(c_intptr_t) :: open_library end function function get_function(handle, name) bind(c, name='get_function') import :: c_intptr_t, c_funptr integer(c_intptr_t), value :: handle character(kind=c_char), intent(in) :: name(*) type(c_funptr) :: get_function end function function get_variable(handle, name) bind(c, name='get_variable') import :: c_intptr_t, c_ptr integer(c_intptr_t), value :: handle character(kind=c_char), intent(in) :: name(*) type(c_ptr) :: get_variable end function subroutine close_library(handle) bind(c, name='close_library') import :: c_intptr_t integer(c_intptr_t), value :: handle end subroutine end interface ! 目标函数类型声明 abstract interface function compute_func(x, y) bind(c) import :: c_double real(c_double), value :: x, y real(c_double) :: compute_func end function end interface ! 变量 integer(c_intptr_t) :: lib_handle type(c_funptr) :: func_ptr type(c_ptr) :: var_ptr procedure(compute_func), pointer :: compute_fptr => null() real(c_double), pointer :: global_value_ptr => null() real(c_double) :: result ! 1. 加载库 lib_handle = open_library('libmylib.so' // c_null_char) if (lib_handle == 0) then print *, 'Error: Failed to load library' stop end if print *, 'Library loaded successfully' ! 2. 获取函数 func_ptr = get_function(lib_handle, 'compute' // c_null_char) if (.not. c_associated(func_ptr)) then print *, 'Error: Failed to get function "compute"' call close_library(lib_handle) stop end if call c_f_procpointer(func_ptr, compute_fptr) ! 3. 获取变量 var_ptr = get_variable(lib_handle, 'global_value' // c_null_char) if (.not. c_associated(var_ptr)) then print *, 'Error: Failed to get variable "global_value"' call close_library(lib_handle) stop end if call c_f_pointer(var_ptr, global_value_ptr) ! 4. 调用函数 & 访问变量 result = compute_fptr(2.0_c_double, 3.0_c_double) print *, 'compute(2.0, 3.0) = ', result print *, 'global_value = ', global_value_ptr ! 5. 修改变量（通过指针） global_value_ptr = 2.71828_c_double result = compute_fptr(1.0_c_double, 1.0_c_double) print *, 'After modification: compute(1.0,1.0) = ', result ! 6. 卸载库 call close_library(lib_handle) print *, 'Library unloaded' end program dynamic_load_example

4. 编译与运行

# 编译 Fortran 程序（链接 libdl）gfortran -o main main.f90 libdlwrapper.a -ldl# 运行（确保 libmylib.so 在 LD_LIBRARY_PATH 中）exportLD_LIBRARY_PATH=.:$LD_LIBRARY_PATH./main

输出示例：

Library loaded successfully compute(2.0, 3.0) = 9.1415900000000003 global_value = 3.1415900000000001 After modification: compute(1.0,1.0) = 3.7182800000000001 Library unloaded

三、Windows 平台适配要点

1. 动态库：编译为.dll，导出函数需加__declspec(dllexport)

   __declspec(dllexport)doublecompute(doublex,doubley){...}

2. 封装层：替换为 Windows API

   #include<windows.h>HMODULEopen_library(constchar*path){returnLoadLibraryA(path);}FARPROCget_function(HMODULE h,constchar*name){returnGetProcAddress(h,name);}

3. Fortran 接口：将intptr_t替换为type(c_ptr)，其他逻辑相同

4. 编译：

   gfortran -o main.exe main.f90 dlwrapper.o -lkernel32

四、关键注意事项

五、高级技巧：直接使用 libdl（无需 C 封装）

部分编译器（如 GCC 10+）支持直接绑定libdl函数：

interface function dlopen(filename, mode) bind(c, name='dlopen') import :: c_char, c_int, c_ptr character(kind=c_char), intent(in) :: filename(*) integer(c_int), value :: mode type(c_ptr) :: dlopen end function end interface

但需注意：

• 需链接-ldl
• RTLD_LAZY等常量需手动定义（如integer(c_int), parameter :: RTLD_LAZY = 1）
• 可移植性较差，推荐使用 C 封装层

六、适用场景

• 插件化架构（运行时加载算法模块）
• 避免硬编码依赖（如可选的 CUDA/MPI 支持）
• 高性能计算中动态切换优化内核
• 与 Python/C/C++ 生态集成（如通过 ctypes 加载 Fortran 库）

💡建议：生产环境推荐使用CMake + 预编译封装库管理依赖，避免手动处理符号和路径问题。