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)