4

我想创建一个派生数据类型,其中包含一个表和该表的长度。理想情况下,我想声明如下

type mydata
    integer                :: ndata
    real, dimension(ndata) :: x
end type mydata

但是一开始 ndata 是未知的,因此 x 是一个动态表。

allocatable当 ndata 已知时,我是否必须使用x 并分配它?

type mydata
    integer                         :: ndata
    real, dimension(:), allocatable :: x
end type mydata

....

type(mydata) :: var
var % ndata = 10
allocate(var % x(10))

处理此类案件的最佳方法是什么?

4

1 回答 1

6

解决问题的最佳方法是使用最新的 Fortran 2003 兼容编译器并使用参数化派生类型。如果我没记错的话,目前只有 Cray 和 IBM 编译器支持这个特性。它将允许您定义这样的类型

type mydata(sz)
    integer, len :: sz
    real, dimension(sz) :: x
end type mydata

这似乎正是你想要的。如果您没有这些编译器中的任何一个,那么临时解决方法(例如您建议的那个)将不得不暂时满足您。当然,您可以将首选的解决方法捆绑到非默认类型构造函数中。

如果您在没有参数化派生类型的情况下陷入困境,那么我将从这样的类型定义开始:

type mydata
    real, dimension(:), allocatable :: x
  contains
    procedure :: setSize
    procedure :: getSize
end type mydata

并以明显的(?)方式定义那些类型绑定的过程。这将防止ndata和大小x变得不同。

于 2013-11-06T11:16:27.223 回答