Question

Fortran内部函数transfer可用于将派生类型转换为实数或整数数组。当在遗留系统中工作时，这可能非常有用，遗留系统依赖于基本类型的数组（整数，实数等）来实现持久性。

以下代码至少在ifort和gfortran上运行，并将简单的派生类型示例转换为整数数组（使用解决方案更新）：

program main
    implicit none

    integer, parameter :: int_mem_size = storage_size(1)

    type subtype
       integer a
       double precision b
    end type subtype

    type :: mytype
         integer :: foo
         double precision :: bar
         type(subtype) :: some_type
    end type

    type(mytype)  :: my_var
    type(subtype) :: my_subtype

    ! Old version: integer :: x(30)
    integer, allocatable :: x(:)
    integer :: mem_size

    !Allocate array with required size
    mem_size = storage_size(my_var)
    allocate(x(mem_size/int_mem_size))

    my_subtype%a = 1
    my_subtype%b = 2.7

    my_var%foo = 42
    my_var%bar = 3.14
    my_var%some_type = my_subtype

    write(*,*) "transfering..."
    x = transfer(my_var, x)
    write(*,*) "Integer transformation:", x

end program main

在我的电脑上，这是输出（此结果至少取决于平台）：

 transfering...
 Integer transformation:                    42           0  1610612736  1074339512         
 999           0    -1610612736  1074108825

我的问题是我已经猜到了＃34;一个30元素的长整数数组足以存储这个数据结构。有没有办法可以确定数组需要多大才能存储整个数据结构？

Answer 1

如果你有一个符合Fortran 2008标准的编译器，或者一个足够兼容的编译器，你会发现内部函数storage_size，它返回用于存储其参数的位数。如果不熟悉我熟悉的大多数编译器都会执行非标准功能来执行此操作;英特尔Fortran编译器有一个名为sizeof的函数，它返回存储其参数所需的字节数。

确定内存中的Fortran派生类型大小

1 个答案: