2

我正在学习 MPI-IO 并遵循教程（此处下载 PDF）。

对于一项练习，正确的代码是：

Program MPI_IOTEST
Use MPI
Implicit None
Integer :: wsize,wrank
Integer :: ierror
Integer :: fh,offset
Call MPI_Init(ierror)
Call MPI_Comm_rank(MPI_COMM_WORLD,wrank,ierror)
Call MPI_Comm_size(MPI_COMM_WORLD,wsize,ierror)
offset=4*wrank; ! because 4 bytes is one signed int
! --- open the MPI files using a collective call
Call MPI_File_Open(MPI_COMM_WORLD,'test.dat',MPI_MODE_RDWR+MPI_MODE_CREATE,MPI_INFO_NULL,fh,ierror);

Write(*,*)'rank',wrank


Call MPI_FILE_WRITE_AT(fh, offset, wrank,1,MPI_INTEGER,mpi_status_ignore,ierror);

Call MPI_File_close(fh,ierror)
Call MPI_Finalize(ierror)
End Program MPI_IOTEST

然后您只需构建并运行它作为 24 个 MPI 任务。然后进行验证，只需执行 od -i test/dat 您将获得与教程中完全相同的结果，如下所示。

0000000           0           1           2           3
0000020           4           5           6           7
0000040           8           9          10          11
0000060          12          13          14          15
0000100          16          17          18          19
0000120          20          21          22          23
0000140

但是，如果我将1更改为num：

Call MPI_FILE_WRITE_AT(fh, offset, wrank,1,MPI_INTEGER,mpi_status_ignore,ierror);

进入

Call MPI_FILE_WRITE_AT(fh, offset, wrank,num,MPI_INTEGER,mpi_status_ignore,ierror);

在此之前定义

integer :: num
num=1

在 rm test.dat 之后，然后重新构建文件并运行它，你会得到：

0000000           0           0           0           0
*
4

1 回答 1

4

您的错误实际上不在规范或使用中，num而是在规范中offset。

如果您阅读 的手册页，MPI_File_write_at则必须将偏移量指定为 MPI_Offset 类型。

因此，如果您更改程序以使用：

integer(kind=MPI_OFFSET_KIND) :: offset

它工作正常。

您没有注意到生成的 test.dat 文件的大小吗？

于 2015-08-04T21:57:17.540 回答