gpt4 book ai didi

performance - 读取未格式化文件的最有效方法

转载 作者:行者123 更新时间:2023-12-04 03:04:20 24 4
gpt4 key购买 nike

现在我正在使用 Fortran 处理 100,000 个文件。这些数据由 HPC 使用 MPI I/O 生成。现在我只能想出以下读取原始数据的方法,效率不高。有没有可能一次读完ut_yz(:,J,K),而不是一个一个读?谢谢
旧代码如下,效率没那么高。

  OPEN(10,FILE=trim(filename)//".dat",FORM='UNFORMATTED',&
ACCESS='DIRECT', RECL=4, STATUS='OLD')
!,CONVERT='big_endian'
COUNT = 1
DO K=1,nz
DO J=1,ny
DO I=1,nxt
READ(10,REC=COUNT) ut_yz(I,J,K)
COUNT = COUNT + 1
ENDDO
ENDDO
ENDDO
CLOSE(10)

想要的是
 OPEN(10,FILE=trim(filename)//".dat",FORM='UNFORMATTED', RECL=4, STATUS='OLD')
!,CONVERT='big_endian'
COUNT = 1
DO K=1,nz
DO J=1,ny
READ(10,REC=COUNT) TEMP(:)
COUNT = COUNT + 153
ut_yz(:,J,K)=TEMP(:)
ENDDO
ENDDO
CLOSE(10)

然而,它总是失败。任何人都可以对此发表评论吗?谢谢。

最佳答案

如果我没记错的话,直接 IO 读取将读取单个记录。因此,在您的新代码版本中,您需要相应地增加记录长度:

   inquire(iolength=rl) ut_yz(:,1,1)
open(10, file=trim(filename)//'.dat', form='UNFORMATTED', recl=rl, status='OLD', action='READ')
count = 1
do k=1,nz
do j=1,ny
read(10, rec=count) ut_yz(:,j,k)
count = count + 1
end do
end do
close(10)

当然,在这个例子中你也可以一次读取完整的数组,这应该是最快的选择:
   inquire(iolength=rl) ut_yz
open(10, file=trim(filename)//'.dat', form='UNFORMATTED', recl=rl, status='OLD', action='READ')
read(10, rec=1) ut_yz
close(10)

关于performance - 读取未格式化文件的最有效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/35433141/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com