Segmentation fault 为什么MPI isend irecv不工作？

我有一个用户定义类型的数组。两个处理器中的每一个都对数组的一部分进行唯一的写入，最后，我希望使用组合数组。为什么我的代码会产生分段错误

program type_derived_recv_send
  use nodeinfo
  implicit none
  include 'mpif.h'
  integer(4) :: ierr
  integer(4) :: istatus1(MPI_STATUS_SIZE)
  integer(4) :: istatus2(MPI_STATUS_SIZE)
  integer(4) :: i,j,k,l
  integer(4) :: iblock(8)
  integer(4) :: idisp(8)
  integer(4) :: itype(8)
  integer(4) :: inewtype
  integer(4) :: iextent
  integer(4) :: itag, isend, irecv
  integer(4) :: n,m,p,q

    TYPE :: newXYZ
        real :: x, u
        real :: y, v
        real :: z, w
        real,dimension(3) :: uvw     
    END TYPE

    TYPE (newXYZ) :: xyz, xyz2

    TYPE (newXYZ), allocatable :: xyzArray(:,:,:,:) 
    n = 1
    m = 1 
    p = 1
    q = 2

    allocate( xyzArray(m,n,p,q) )


  call MPI_INIT(ierr)
  icomm = MPI_COMM_WORLD
  call MPI_COMM_SIZE(icomm,nproc,ierr)
  call MPI_COMM_RANK(icomm,iproc,ierr)

do l=1,q
    do k=1,p
        do j=1,n
            do i=1,m 
                    xyzArray(i,j,k,m)%x = 0
                    xyzArray(i,j,k,m)%y = 0
                    xyzArray(i,j,k,m)%z = 0
                    xyzArray(i,j,k,m)%u = 0
                    xyzArray(i,j,k,m)%v = 0
                    xyzArray(i,j,k,m)%w = 0
                    xyzArray(i,j,k,m)%uvw = (/0,0,0/)
end do
end do
end do
end do

    if (iproc == 1) then
        xyzArray(1,1,1,2)%x = 1.1
        xyzArray(1,1,1,2)%y = 2.1
        xyzArray(1,1,1,2)%z = 3.1
        xyzArray(1,1,1,2)%u = 4.1
        xyzArray(1,1,1,2)%v = 5.1
        xyzArray(1,1,1,2)%w = 6.1
        xyzArray(1,1,1,2)%uvw = (/10.1,10.1,10.1/)
    else
        xyzArray(1,1,1,1)%x = 0.1
        xyzArray(1,1,1,1)%y = 0.1       
        xyzArray(1,1,1,1)%z = 0.1
        xyzArray(1,1,1,1)%u = 0.1
        xyzArray(1,1,1,1)%v = 0.1       
        xyzArray(1,1,1,1)%w = 0.1 
        xyzArray(1,1,1,1)%uvw = (/0.1,0.1,0.1/)
    endif

    itag = 1
    xyz = xyzArray(1,1,1,2)
    xyz2 = xyzArray(1,1,1,1)

    if ( iproc == 1) call MPI_ISEND(xyzArray(1,1,1,2),36,MPI_INTEGER,0,1,icomm,isend,ierr)
    if ( iproc == 0) call MPI_IRECV(xyzArray(1,1,1,2),36,MPI_INTEGER,1,1,icomm,irecv,ierr)
    call MPI_WAIT(isend,istatus1,ierr)
    call MPI_WAIT(irecv,istatus2,ierr)

     if ( iproc == 0 )write(6,*)'iproc = ',iproc
     if ( iproc == 0 ) write(6,*)xyzArray

     if ( iproc == 1 )write(6,*)'iproc = ',iproc
     if ( iproc == 1 ) write(6,*)xyzArray

  call MPI_FINALIZE(ierr)
end program type_derived_recv_send

输出：

--------------------------------------------------------------------------
 newtype extent =           36
 xyz =    0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000    
 xyz2 =   0.1000000      0.1000000      0.1000000      0.1000000      0.1000000      0.1000000      0.1000000      0.1000000      0.1000000    
[west0334:09950] *** Process received signal ***
[west0334:09950] Signal: Segmentation fault (11)
[west0334:09950] Signal code: Address not mapped (1)
[west0334:09950] Failing at address: 0x40
[west0334:09950] [ 0] /lib64/libpthread.so.0 [0x3e00e0eb70]
[west0334:09950] [ 1] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0(ompi_request_wait+0x17) [0x2b33f8628607]
[west0334:09950] [ 2] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0(PMPI_Wait+0x88) [0x2b33f86566e8]
[west0334:09950] [ 3] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi_f77.so.0(PMPI_WAIT+0x3a) [0x2b33f83eac8a]
[west0334:09950] [ 4] ./type_derived_recv_send.x(MAIN__+0x730) [0x401cd0]
 newtype extent =           36
 xyz =    1.100000       4.100000       2.100000       5.100000       3.100000       6.100000       10.10000       10.10000       10.10000    
 xyz2 =    0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000    
[west0334:09951] *** Process received signal ***
[west0334:09951] Signal: Segmentation fault (11)
[west0334:09951] Signal code: Address not mapped (1)
[west0334:09951] Failing at address: 0x40
[west0334:09950] [ 5] ./type_derived_recv_send.x(main+0xe) [0x4049ae]
[west0334:09950] [ 6] /lib64/libc.so.6(__libc_start_main+0xf4) [0x3e0021d994]
[west0334:09950] [ 7] ./type_derived_recv_send.x [0x4014e9]
[west0334:09950] *** End of error message ***
[west0334:09951] [ 0] /lib64/libpthread.so.0 [0x3e00e0eb70]
[west0334:09951] [ 1] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0(ompi_request_wait+0x17) [0x2b73c1aca607]
[west0334:09951] [ 2] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0(PMPI_Wait+0x88) [0x2b73c1af86e8]
[west0334:09951] [ 3] /n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi_f77.so.0(PMPI_WAIT+0x3a) [0x2b73c188cc8a]
[west0334:09951] [ 4] ./type_derived_recv_send.x(MAIN__+0x748) [0x401ce8]
[west0334:09951] [ 5] ./type_derived_recv_send.x(main+0xe) [0x4049ae]
[west0334:09951] [ 6] /lib64/libc.so.6(__libc_start_main+0xf4) [0x3e0021d994]
[west0334:09951] [ 7] ./type_derived_recv_send.x [0x4014e9]
[west0334:09951] *** End of error message ***
Nov 16 20:17:14 2012 9939 4 7.06 handleTSRegisterTerm(): TS reports task <0> pid <9950> on host<west0334> killed or core dumped
Job  /lsf/7.0/linux2.6-glibc2.3-x86_64/bin/openmpi_wrapper ./type_derived_recv_send.x

--------------------------------------------------------------------------
新类型范围=36
xyz=0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.000000
xyz2=0.1000000 0.1000000 0.1000000 0.1000000 0.1000000 0.1000000 0.1000000 0.1000000 0.1000000
[west0334:09950]***处理接收到的信号***
[west0334:09950]信号：分段故障（11）
[west0334:09950]信号代码：地址未映射（1）
[west0334:09950]在地址0x40处失败
[west0334:09950][0]/lib64/libpthread.so.0[0x3e00e00eb70]
[west0334:09950][1]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0（ompi请求等待+0x17）[0x2b33f8628607]
[west0334:09950][2]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0（PMPI_Wait+0x88）[0x2b33f86566e8]
[west0334:09950][3]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi_f77.so.0（PMPI_WAIT+0x3a）[0x2B33f838a]
[west0334:09950][4]。/键入导出的接收发送.x（MAIN+0x730）[0x401cd0]
新类型范围=36
xyz=1.100000 4.100000 2.100000 5.100000 3.100000 6.100000 10.10000 10.10000 10.10000
xyz2=0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.000000
[west0334:09951]***处理接收到的信号***
[west0334:09951]信号：分段故障（11）
[west0334:09951]信号代码：地址未映射（1）
[west0334:09951]在地址0x40处失败
[west0334:09950][5]。/键入导出的接收发送.x（main+0xe）[0x4049ae]
[west0334:09950][6]/lib64/libc.so.6（u libc_start_main+0xf4）[0x3e0021d994]
[west0334:09950][7]。/键入导出的接收发送.x[0x4014e9]
[west0334:09950]***错误消息结束***
[west0334:09951][0]/lib64/libpthread.so.0[0x3e00e0eb70]
[west0334:09951][1]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0（ompi请求等待+0x17）[0x2b73c1aca607]
[west0334:09951][2]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi.so.0（PMPI_Wait+0x88）[0x2b73c1af86e8]
[west0334:09951][3]/n/sw/openmpi-1.2.5-gcc-4.1.2/lib64/libmpi_f77.so.0（PMPI_WAIT+0x3a）[0x2b73c188cc8a]
[west0334:09951][4]。/键入派生的接收发送.x（MAIN+0x748）[0x401ce8]
[west0334:09951][5]。/键入导出的接收发送.x（main+0xe）[0x4049ae]
[west0334:09951][6]/lib64/libc.so.6（libc\u start\u main+0xf4）[0x3e0021d994]
[west0334:09951][7]。/键入导出的接收发送.x[0x4014e9]
[west0334:09951]***错误消息结束***
2012年11月16日20:17:14 9939 4 7.06 handleTSRegisterTerm（）：TS报告主机上的任务pid已终止或内核已转储
Job/lsf/7.0/linux2.6-glibc2.3-x86\u 64/bin/openmpi\u包装器。/type\u派生\u recv\u send.x

为什么要发送带有9个实数的类型和整数MPI数据类型？为什么不发送9个MPI_实数或36个MPI_字节

为什么要发送带有9个实数和整数MPI数据类型的类型？为什么不发送9 MPI_实字节或36 MPI_字节？错误如下：

if ( iproc == 1) call MPI_ISEND(xyzArray(1,1,1,2),36,MPI_INTEGER,0,1,icomm,isend,ierr)
if ( iproc == 0) call MPI_IRECV(xyzArray(1,1,1,2),36,MPI_INTEGER,1,1,icomm,irecv,ierr)
call MPI_WAIT(isend,istatus1,ierr)
call MPI_WAIT(irecv,istatus2,ierr)

irecv

请求句柄在秩

0

中保持未初始化状态，而

isend

请求句柄在秩

1

中保持单元化状态。这里有两个选项-要么将

MPI\u WAIT

调用与其相应的非阻塞MPI操作分组，要么在此之前将

isend

和

irecv

初始化为

MPI\u REQUEST\u NULL

解决方案1：

if ( iproc == 1) then
   call MPI_ISEND(xyzArray(1,1,1,2),36,MPI_INTEGER,0,1,icomm,isend,ierr)
   call MPI_WAIT(isend,istatus1,ierr)
end if
if ( iproc == 0) then
   call MPI_IRECV(xyzArray(1,1,1,2),36,MPI_INTEGER,1,1,icomm,irecv,ierr)
   call MPI_WAIT(irecv,istatus2,ierr)
end if

解决方案2（简化）：

您实际上不需要两个单独的请求句柄以及两个单独的状态数组。通过使用空句柄初始化

ireq

，此代码也可以在2个以上MPI进程中正常运行

请注意，为了便于移植，您应该使用

MPI\u TYPE\u CREATE\u STRUCT

构造派生结构数据类型，并使用它发送和接收

TYPE（newXYZ）

的元素。错误如下：

if ( iproc == 1) call MPI_ISEND(xyzArray(1,1,1,2),36,MPI_INTEGER,0,1,icomm,isend,ierr)
if ( iproc == 0) call MPI_IRECV(xyzArray(1,1,1,2),36,MPI_INTEGER,1,1,icomm,irecv,ierr)
call MPI_WAIT(isend,istatus1,ierr)
call MPI_WAIT(irecv,istatus2,ierr)

irecv

请求句柄在秩

0

中保持未初始化状态，而

isend

请求句柄在秩

1

中保持单元化状态。这里有两个选项-要么将

MPI\u WAIT

调用与其相应的非阻塞MPI操作分组，要么在此之前将

isend

和

irecv

初始化为

MPI\u REQUEST\u NULL

解决方案1：

if ( iproc == 1) then
   call MPI_ISEND(xyzArray(1,1,1,2),36,MPI_INTEGER,0,1,icomm,isend,ierr)
   call MPI_WAIT(isend,istatus1,ierr)
end if
if ( iproc == 0) then
   call MPI_IRECV(xyzArray(1,1,1,2),36,MPI_INTEGER,1,1,icomm,irecv,ierr)
   call MPI_WAIT(irecv,istatus2,ierr)
end if

解决方案2（简化）：

您实际上不需要两个单独的请求句柄以及两个单独的状态数组。通过使用空句柄初始化

ireq

，此代码也可以在2个以上MPI进程中正常运行

---

请注意，为了便于移植，您应该使用

MPI\u TYPE\u CREATE\u STRUCT

构造派生结构数据类型，并使用它发送和接收

TYPE（newXYZ）

@Pippi的元素，欢迎使用。我想您知道如何定义派生的MPI数据类型，因为您已经声明了必要的变量（

iblock

，

idisp

，

itype

，等等），我只想指出，不推荐使用采用普通整数置换的结构构造函数

MPI\u TYPE\u STRUCT

，您应该使用

MPI\u TYPE\u CREATE\u STRUCT

。它具有相同的参数，但位移类型为

INTEGER（KIND=MPI\u ADDRESS\u KIND）

。此外，请不要使用

包含“mpif.h”

。这是非常古老的FORTRAN 77接口，它一定已经死了。应改用Fortran 90+接口：

使用mpi

。甚至更好的是，一旦完整的MPI-3.0实现可用，就应该使用MPI-3.0的Fortran 2008接口（

使用MPI_f08

）。感谢您的建议@皮皮，不客气。我想您知道如何定义派生MPI