Fortran: KIND=8数组到KIND=4的高效转换

Question

我使用的是一个外部fortran库，它以REAL(kind=8)类型返回大量数据。但是，我需要将数据写成REAL(kind=4)。

将大型数据数组转换为不同KIND类型的最佳方法是什么？

我现在只是天真地这么做：

REAL(KIND=8), ALLOCATABLE :: data64(:,:)
REAL(KIND=4), ALLOCATABLE :: data32(:,:)

# allocate arrays and call external function
CALL some_external_function(data64)

# Convert to 32bit using simple assignment
data32 = data64

# Write out 32bit data....
WRITE(UOUT) data32

然而，这种方法要求我分配两个大数组而不是一个数组，并对数据进行大量的复制--在内存和处理方面似乎都有点效率低下.

有更聪明的方法吗？

顺便说一句:矩阵的大小不一定很大--它们变化很大，但很可能是1000×1000左右(但这将在一个循环中发生100次)。不是现代机器的问题，但正如其中的评论之一，我正在寻找一个优雅的解决方案，也是有效的！只做seems.....lazy和inelagant的指定？

Answer 1

这是给你的优雅

write(uout) real(data64,kind(data32))