mpi传递fortran对象问题

wangmianzhi

问题：
* fortran要求必须是non-sequencial的type才能扩展成对象
* mpi要求必须是sequencial的type才能通过计算偏移包装成mpi扩展类型

这就造成要想传递对象，必须：
1. 把对象的数据拷给一个不带procedure的sequencial type
2. 计算偏移，制作自定义mpi类型
3. 发送，用同样的sequencial type接收
4. 拷到新的对象里
这个过程太恶心了。有没有体面一些的方法？
多谢

wangmianzhi

问题已解决。
mpi可以直接传递对象的数据
相关解释在这里：
http://www.mpi-forum.org/docs/mpi21-report/node337.htm#Node341

附上一个例子：

1. module mp
   
2.   type::typeA
   
3.     integer::v1(3)
   
4.     double precision::v2(2)
   
5.   contains
   
6.     procedure::show
   
7.   end type
   
8. contains
   
9.   subroutine show(this)
   
10.     class(typeA)::this
    
11.     write(*,*),this%v1,this%v2
    
12.   end subroutine
    
13. end module
    
14. 
    
15. program test
    
16.   use mp
    
17.   use mpi
    
18.   
    
19.   integer mpiErr,iP,mpiStat(MPI_status_size),mpiLen(2),mpiType(2)
    
20.   integer(MPI_address_kind)::mpiDisp(2),mpiBase
    
21.   type(typeA)::aa
    
22.   
    
23.   call MPI_init(mpiErr)
    
24.   call MPI_comm_rank(MPI_comm_world,iP,mpiErr)
    
25.   
    
26.   mpiLen=[3,2]
    
27.   mpiType=[MPI_integer,MPI_double_precision]
    
28.   call MPI_get_address(aa%v1,mpiDisp(1),mpiErr)
    
29.   call MPI_get_address(aa%v2,mpiDisp(2),mpiErr)
    
30.   mpiBase=mpiDisp(1)
    
31.   mpiDisp(:)=mpiDisp(:)-mpiBase
    
32.   call MPI_type_create_struct(2,mpiLen,mpiDisp,mpiType,MPItypeA,mpiErr)
    
33.   call MPI_type_commit(MPItypeA,mpiErr)
    
34.   
    
35.   if(iP==0)then
    
36.     aa%v1=[1,2,3]
    
37.     aa%v2=[1,2]
    
38.     call MPI_send(aa,1,MPItypeA,1,1,MPI_comm_world,mpiErr)
    
39.   else
    
40.     call MPI_recv(aa,1,MPItypeA,0,1,MPI_comm_world,mpiStat,mpiErr)
    
41.     call aa%show
    
42.     write(*,*),iP
    
43.   end if
    
44.   
    
45.   call MPI_finalize(mpiErr)
    
46. end program

复制代码