3

我使用 mpi4py 在几个 procs 之间拆分一些计算。基本上我只是计算一些凸包的体积，我使用 tvtk 和 mayavi 创建。

只有第一个 proc 导入这些库：

...
if rank==0:
   from tvtk.api import tvtk
   from mayavi Import mlab
   ...
   mlab.figure(size=(1024,768),bgcolor=(1,1,1))
   ...

然后，我尝试在所有 proc 之间共享对象 mlab 和 tvtk：

for i in range(comm.Get_size()):
    comm.send(mlab,dest=i)
    comm.send(tvtk,dest=i)
....

以下步骤看起来像......

Points=local_data
ug=tvtk.UnstructuredGrid(Points=Points)
...
dataname="Data %s " % rank
ds=mlab.pipeline.add_dataset(ug,name=dataname)
delaunay=mlab.pipeline.delaunay3d(ds,name=dataname)
... calc volume...

但是，显然不可能发送实例/类（或任何 mlab 和 tvtk ），因为我总是收到以下错误：

comm.send(mlab,dest=i)
File "Comm.pyx", line 753, in mpi4py.MPI.Comm.send (src/mpi4py.MPI.c:53848)
File "pickled.pxi", line 122, in mpi4py.MPI.PyMPI_send (src/mpi4py.MPI.c:20409)
 File "pickled.pxi", line 39, in mpi4py.MPI._p_Pickle.dump (src/mpi4py.MPI.c:19503)
cPickle.PicklingError: Can't pickle <type 'module'>: attribute lookup __builtin__.module failed

有没有办法在所有 procs 之间“共享”实例 mlab 和 tvtk？

编辑：简短的例子；你能让这个工作吗？

from mpi4py import MPI
comm=MPI.COMM_WORLD

size=comm.Get_size()
rank=comm.Get_rank()

if rank==0:
        from tvtk.api import tvtk
        from mayavi import mlab

if __name__=='__main__':
        if rank==0:
                for i in range(size):
                        comm.send(mlab,dest=i)
                        comm.send(tvtk,dest=i)
        else:
                local_mlab=comm.recv(mlab,source=0)
                local_tvtk=comm.recv(tvtk,source=0)
4

2 回答 2

3

我想以下可能适用于您的情况

from mpi4py import MPI
MPI._p_pickle.dumps = dill.dumps
MPI._p_pickle.loads = dill.loads
于 2014-02-13T19:54:25.227 回答
0

我认为问题是你不能以这种方式发送包，你可以在所有正在运行的实例上导入它们并在需要的地方使用它们，这似乎可以让你做你想做的事

但是，这不会共享任何全局模块数据，这可能是您正在尝试执行的操作，因此这需要您单独发送它，我建议使用： 

global_data = None

if rank == 0:
    global_data = {'value_of_interest': mlab.some_parameter}

global_data = MPI.COMM_WORLD.bcast(global_data, root=0)
于 2014-09-19T10:43:57.957 回答