arrays - 将派生类型作为数组传递

Question

在 Fortran 中，可以对数组进行操作，但是如何将派生类型的索引也视为数组的一部分呢？代码将解释我最想做的事情:

type mytype
    integer :: b(3,3)
    real :: c(4)
endtype

integer :: a(3,3)
real :: d(2,4)
type(mytype) :: mat(2)

!do stuff so that 'mat' gets values
....

!usually one does this
a = matmul(mat(1)%b, transpose(mat(2)%b))

!multiplying two 3x3 matrices

!but how does one do this? Note the "array"
d = matmul(mat(:)%c, mat(:)%c)

我假设最后一行类似于一个 2x4 矩阵与自身相乘。但是，当我尝试编译时，gfortran 提示

Error: Two or more part references with nonzero rank must not be specified

这可以用 Fortran 实现吗？

Answer 1

您希望编译器将 mat(:)%c 视为 2 x 4 矩阵？它不是那样工作的。 mat 和 c 是不同的对象，它们的等级不会合并到一个数组中。 mat 是用户定义的类型，c 是实数矩阵。仅仅因为您只使用 mat 的 c 组件并不意味着编译器会将 c 提升为更高维的实数数组，基于在 mat 的维度上。

您可以通过 X = [ mat(1)%c, mat(2)%c ] 创建一个新数组。您可以使用 reshape 来控制形状。