1

我正在编写一段代码，它将使用mpi4py. 我收到一个给我带来问题的错误。

我放了 2 个默认变量m, n，用户将数字插入到rank 0过程中以填充矩阵。该过程rank 0将矩阵发送到process 1和process 2。在process 1我初始化矩阵和行列式，计算行列式，并将其发送回rank 0. 在process 2中，我初始化矩阵，找到逆矩阵，并将逆矩阵发送回进程 0。

然后我process 0收到行列式和逆矩阵。如果行列式为0，则发送消息并退出程序（即使行列式为0也可以计算逆矩阵，但不正确）。如果不为 0，程序将打印逆矩阵。

请注意，我知道当前的实现不是并行化此问题的最佳方法，但我需要从某个地方开始。

行列式给了我一些错误，所以我将它初始化为一个包含 2 个元素的空numpy数组，其中第一个元素是行列式，第二个元素是 0。我将代码从我的母语编辑为英语，以便可以阅读和更容易理解，因此它可能包含一些错误。

问题是当我想给用户写矩阵大小的权限时。我检查了一些相关的答案，并尝试使用map()、raw_input()和其他选项，但没有一个有效。

编码：

import numpy as np
from numpy.linalg import inv
from scipy.linalg import det
from mpi4py import MPI



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

m = int(input())
n = int(input())

if rank==0:
    matrix = np.zeros((m,n))
    determinant = np.zeros((2))
    for i in range (0,m):
        for j in range(0,n):
            print("Enter the value for the field: m = ", i+1, ", n = ", j+1,"\n")
            matrix[i][j] = float(input())

    comm.Send(matrix, dest=1, tag=0)
    comm.Send(matrix, dest=2, tag=0)
    comm.Recv(determinant, source=1, tag=0)
    comm.Recv(matrix, source=2, tag=0)
    if(determinant[0]==0):
        print("There is no inverse matrix since the determinant is 0!")
    else:
        print("Inverse matrix:\n",matrica)
elif rank==1:
    matrix = np.zeros((m,n))
    determinant = np.zeros((2))
    comm.Recv(matrix, source=0, tag=0)
    determinant = np.array([[det(matrix)],[0]])
    comm.Send(determinant, dest=0)
elif rank==2:
    matrix = np.zeros((m,n))
    comm.Recv(matrix, source=0, tag=0)
    matrix = inv(matrix)
    comm.Send(matrix, dest=0)
else:
    exit()

错误：

Traceback (most recent call last):
  File "lambda.py", line 13, in <module>
    m = int(input())
EOFError: EOF when reading a line
Traceback (most recent call last):
  File "lambda.py", line 13, in <module>
    m = int(input())
EOFError: EOF when reading a line

编辑：我正在使用 PuTTY 从大学连接到 Debian 操作系统。

解释器版本：sys.version_info(major=3, minor=5, micro=1, releaselevel='final', serial=0)

我输入命令行：mpirun -np 3 python3 lambda.py，我得到一个 EoFError 并插入数字直到它应该做某事的地方，而不是这样做，它继续作为无限循环运行。（更好看在图片中） 图片

添加了跟踪输入的代码：

m = int(input("Enter number of rows, m = \n"))
n = int(input("Enter number of columns, n = \n"))

图二

4

1 回答 1

0

你的问题与矩阵无关：这是你的线

m = int(input())

正在所有节点上发生。您需要在根节点上执行此操作，然后以您喜欢的方式共享结果，例如通过广播。例如，如果我们从

from mpi4py import MPI

x = input("enter x: ")

我们得到

(3.5) dsm@winter:~/coding/mpipar$ mpirun -np 3 python m1.py
enter x: enter x: Traceback (most recent call last):
  File "m1.py", line 3, in <module>
    x = input("enter x: ")
EOFError: EOF when reading a line
enter x: Traceback (most recent call last):
  File "m1.py", line 3, in <module>
    x = input("enter x: ")
EOFError: EOF when reading a line

并注意有两个错误。但是如果我们保护它：

from mpi4py import MPI

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

x = None
if rank == 0:
    x = input("enter x: ")

x = comm.bcast(x, root=0)

print("rank", rank, "thinks x is", x)

然后我们得到

(3.5) dsm@winter:~/coding/mpipar$ mpirun -np 3 python m2.py
enter x: abc
rank 0 thinks x is abc
rank 1 thinks x is abc
rank 2 thinks x is abc
于 2016-02-22T00:53:11.200 回答