2

我正在尝试在 Fortran 中创建一个数组数组。

类似于以下内容

TYPE :: array_of_arrays
   REAL, DIMENSION(:), POINTER :: p => NULL()
END TYPE

TYPE(array_of_arrays), DIMENSION(2) :: some_array

这样我就可以做到：

REAL, DIMENSION(3), TARGET :: some_vector1 = (/1.0, 2.1, 4.3/)
REAL, DIMENSION(3), TARGET :: some_vector2 = (/3.0, 1.2, 9.6/)

some_array(1)%p => some_vector1
some_array(2)%p => some_vector2

WRITE(*,*) some_array(1)%p ! I see some_vector1
WRITE(*,*) some_array(2)%p ! I see some_vector2

现在，实际上将这些some_vector数组中的每一个声明为对应于我的数组数组中的每个元素对我来说是很麻烦的。

我想做的是在一个单独的子例程中设置一个临时向量作为目标，并且该子例程设置我的数组数组以指向该临时向量。

这样我就可以拥有匿名数组。

但是，这似乎不起作用，我想知道是否首先我正在做一些 Fortran 不支持的事情。

那么 Fortran 是否支持匿名数组，也就是说（如果我的术语有误），一个只能通过引用访问的数组？

4

1 回答 1

4

当然; 正如 IanH 建议的那样，您可以让指针直接引用分配的内存，而不是引用变量；这是分配的内存一旦超出范围就不会自动释放的少数情况之一。

例如，

module arrays


TYPE :: array_of_arrays
   REAL, DIMENSION(:), POINTER :: p => NULL()
END TYPE

contains

subroutine alloc(aa)
    type(array_of_arrays), intent(inout) :: aa(:)
    integer :: i

    allocate( aa(1)%p(1) )
    aa(1) % p = [1.]

    allocate( aa(2)%p(5) )
    aa(2) % p = [ (i, i=1,5) ]

end subroutine alloc

end module arrays



program usearrays
    use arrays

    TYPE(array_of_arrays), DIMENSION(2) :: some_array

    call alloc(some_array)

    WRITE(*,*) some_array(1)%p ! I see some_vector1
    WRITE(*,*) some_array(2)%p ! I see some_vector2

    deallocate( some_array(1) )
    deallocate( some_array(2) )

end program usearrays

并运行它给出

$ gfortran -o arrays arrays.f90 
$ ./arrays 
   1.0000000    
   1.0000000       2.0000000       3.0000000       4.0000000       5.0000000
于 2012-08-17T03:36:47.970 回答