Fortran90+版本实数数组字节交换内置函数及实现方法问询

内置函数说明 #

标准Fortran 90及以上版本均未提供标准化的内置字节交换函数。部分编译器提供非标准扩展实现，例如Intel Fortran的BSWAP、GNU Fortran的bswap系列函数，但这类扩展不具备跨编译器可移植性，需要兼容多编译环境的代码不建议使用。

基于Do循环的可移植实现 #

实现逻辑与数组维度无关，核心思路为逐元素处理：将单个实数通过TRANSFER函数转为等长的单字节整型数组，倒序交换字节后再转回实数类型，遍历全部元素即可完成操作。

支持任意维度的版本（Fortran 2003及以上） #

使用假定秩语法arr(..)，直接支持1/2/3等任意维度的实数数组输入，无需额外转换。

单精度（4字节）实数数组版本

subroutine bswap_real32(arr)
    use iso_fortran_env, only: int8, real32
    implicit none
    real(real32), intent(inout) :: arr(..)
    integer(int8) :: bytes(4), tmp
    integer :: i

    do i = 1, size(arr)
        bytes = transfer(arr(i), bytes)
        ! 交换字节顺序
        tmp = bytes(1)
        bytes(1) = bytes(4)
        bytes(4) = tmp
        tmp = bytes(2)
        bytes(2) = bytes(3)
        bytes(3) = tmp
        arr(i) = transfer(bytes, arr(i))
    end do
end subroutine bswap_real32

双精度（8字节）实数数组版本

subroutine bswap_real64(arr)
    use iso_fortran_env, only: int8, real64
    implicit none
    real(real64), intent(inout) :: arr(..)
    integer(int8) :: bytes(8), tmp
    integer :: i, j

    do i = 1, size(arr)
        bytes = transfer(arr(i), bytes)
        ! 对称交换字节
        do j = 1, 4
            tmp = bytes(j)
            bytes(j) = bytes(9 - j)
            bytes(9 - j) = tmp
        end do
        arr(i) = transfer(bytes, arr(i))
    end do
end subroutine bswap_real64

Fortran 90/95兼容版本 #

如果编译环境不支持Fortran 2003的假定秩特性，可以使用1维数组接口，多维数组调用时通过reshape转换即可：

subroutine bswap_real32_f90(arr, arr_len)
    use iso_fortran_env, only: int8, real32
    implicit none
    integer, intent(in) :: arr_len
    real(real32), intent(inout) :: arr(arr_len)
    integer(int8) :: bytes(4), tmp
    integer :: i

    do i = 1, arr_len
        bytes = transfer(arr(i), bytes)
        tmp = bytes(1)
        bytes(1) = bytes(4)
        bytes(4) = tmp
        tmp = bytes(2)
        bytes(2) = bytes(3)
        bytes(3) = tmp
        arr(i) = transfer(bytes, arr(i))
    end do
end subroutine bswap_real32_f90

! 调用示例：处理2维单精度数组 real(real32) :: arr_2d(10, 20)
! call bswap_real32_f90(reshape(arr_2d, [size(arr_2d)]), size(arr_2d))

注意事项 #

• 代码默认使用iso_fortran_env模块的类型常量，是Fortran标准支持的可移植写法，无需依赖编译器扩展。
• 如需要处理其他长度的实数类型，仅需调整字节数组长度和交换循环的边界即可，核心逻辑不变。
• 逐元素处理不会修改原数组的维度和形状，处理完成后可直接使用原数组。

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