mpi4py 尝试发送类时出现酸洗错误



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

只有第一个进程导入这些库:

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

然后,我尝试在所有进程之间共享对象 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

有没有办法在所有进程之间"共享"实例 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)
我想

以下内容可能适用于您的情况

from mpi4py import MPI
MPI._p_pickle.dumps = dill.dumps
MPI._p_pickle.loads = dill.loads

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

但是,这

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

global_data = None
if rank == 0:
    global_data = {'value_of_interest': mlab.some_parameter}
global_data = MPI.COMM_WORLD.bcast(global_data, root=0)

相关内容

  • 没有找到相关文章

最新更新