0

当我尝试构建一个简单的数组并使用 MPIIO 写入输出文件时，输出文件非常大且不正确。简单数组“buff”为 2 * 3，每个元素都等于 1。我首先创建了 1 * 1 的子数组文件类型，然后创建了一个新的输出文件“out.dat”。我使用 set_view 函数和 write_all 来编写这个 buff 数组。我编译它并使用 6 个内核运行它，以便一个内核编写一个元素。我希望我能得到一个“out.dat”和一个值为 1 的 2*3 数组。

我成功编译运行这段代码，结果都是0，而且体积很大。

program test
use mpi
implicit none
integer::rank,nproc,ierr,buffsize,status(MPI_STATUS_SIZE),intsize,i,j,filetype,cart_comm
integer::fh
integer(kind=mpi_offset_kind):: offset=0
integer::buff(2,3)
character:: filename*50
integer::sizes(2)
integer::gsize(2)
integer::start(2)
integer::subsize(2)
integer::coords(2)
integer:: nprocs_cart(2)=(/2,3/)
logical::periods(2)
intsize=sizeof(rank)
call MPI_init(ierr)
call MPI_COMM_SIZE(MPI_COMM_WORLD, nproc, ierr)
call MPI_COMM_RANK(MPI_COMM_WORLD, rank, ierr)
CALL MPI_Dims_create(nproc, 2, nprocs_cart, ierr)
CALL MPI_Cart_create(MPI_COMM_WORLD, 2, nprocs_cart, periods, .TRUE., &
cart_comm, ierr)
CALL MPI_Comm_rank(cart_comm, rank, ierr)
CALL MPI_Cart_coords(cart_comm, rank, 2, coords, ierr)


start=coords
gsize=(/2,3/)
subsize=(/1,1/)
buff=1


call MPI_TYPE_CREATE_SUBARRAY(2,gsize,subsize,start,MPI_ORDER_FORTRAN,&
MPI_int,filetype,ierr)
call MPI_TYPE_COMMIT(filetype,ierr)

call MPI_File_open(MPI_COMM_WORLD,'out.dat',&
        MPI_MODE_WRONLY + MPI_MODE_CREATE, MPI_INFO_NULL, fh,ierr)


        call MPI_File_set_view(fh,0,MPI_int,filetype,&
        "native",MPI_INFO_NULL,ierr)
  
        CALL MPI_FILE_WRITE_all(fh, buff,6, filetype, MPI_STATUS_ignore, ierr)


call MPI_File_close(fh,ierr)
call MPI_FINALIZE(ierr)

结束程序测试

4

0 回答 0