Question

我正在学习MPI-IO并遵循教程(PDF download here)。

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

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
*

Answer 1

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

如果您阅读MPI_File_write_at的{{3}}，则必须将偏移量指定为MPI_Offset类型。

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

integer(kind=MPI_OFFSET_KIND) :: offset

工作正常。

您是否注意到生成的test.dat文件的大小？

为什么这两个MPI-IO代码的工作方式不同？

1 个答案: