Question

我正在尝试使用以下代码在python上测试一个简单的mpi代码：

from scipy.sparse import csr_matrix
from mpi4py import MPI

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

if rank == 0:
    data = [1, 2, 3, 4, 5]
    indices = [1, 3, 2, 1, 0]
    indptr = [0, 2, 3, 4, 5]
    #A=csr_matrix((data,indices,indptr),shape=(4,4))                                                                                                              

data=comm.bcast(data, root=0)
indices=comm.bcast(indices, root=0)
indptr=comm.bcast(indptr, root=0)

print rank,data,indices,indptr

返回以下错误：

Traceback (most recent call last):
  File "test.py", line 14, in <module>
    data=comm.bcast(data, root=0)
NameError: name 'data' is not defined
Traceback (most recent call last):
  File "test.py", line 14, in <module>
    data=comm.bcast(data, root=0)
NameError: name 'data' is not defined
Traceback (most recent call last):
  File "test.py", line 14, in <module>
    data=comm.bcast(data, root=0)
NameError: name 'data' is not defined
0 [1, 2, 3, 4, 5] [1, 3, 2, 1, 0] [0, 2, 3, 4, 5]
-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[10263,1],1]
  Exit code:    1

似乎错误是由于我没有正确使用comm.bcast，但这正是它在文档中使用的方式。

Answer 1

您正在data块中定义if。 if阻止为false时会发生什么？变量data未定义。

from scipy.sparse import csr_matrix
from mpi4py import MPI

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

data = []
indices = []
indptr = []

if rank == 0:
    data = [1, 2, 3, 4, 5]
    indices = [1, 3, 2, 1, 0]
    indptr = [0, 2, 3, 4, 5]
    #A=csr_matrix((data,indices,indptr),shape=(4,4))                                                                                                              

data=comm.bcast(data, root=0)
indices=comm.bcast(indices, root=0)
indptr=comm.bcast(indptr, root=0)

print rank,data,indices,indptr

现在应该可以了。

comm.bcast不能正常工作

1 个答案: