Python Scipy线性构造器数据类型未指定
我正在尝试使用Scipy的Python Scipy线性构造器数据类型未指定,python,numpy,matrix,scipy,Python,Numpy,Matrix,Scipy,我正在尝试使用Scipy的gmres命令。我的矩阵密集且较大，因此我的计划是使用LinearOperator命令仅返回矩阵向量积。请看我前面的问题。在这个问题的帮助下，我能够构建一个LinearOperator对象，它成功地进行了计算a*x，其中a是矩阵，x是向量 问题是当我调用gmres时。我运行命令： x, info = scipy.sparse.linalg.gmres(A, b) 并返回运算符A没有dtype的错误。这是正确的，因为A.dtype返回错误。我的问题是我不知道为dtype
gmresLinearOperatorLinearOperatora*xaxgmresx, info = scipy.sparse.linalg.gmres(A, b)
AdtypeA.dtypedtypedtype我试图传递
dtype='float64'尝试
dtype=NoneAdtypeA.dtype=NoneA.dtypeAA = sparse.scipy.linalg.linearoperator(shape = (n,n), matvec = mv, dtype = 'float64')
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/usr/lib/python2.7/dist-packages/scipy/sparse/linalg/isolve/iterative.py", line 393, in gmres
A,M,x,b,postprocess = make_system(A,M,x0,b,xtype)
File "/usr/lib/python2.7/dist-packages/scipy/sparse/linalg/isolve/utils.py", line 119, in make_system
M = LinearOperator(A.shape, matvec=psolve, rmatvec=rpsolve, dtype=A.dtype)
AttributeError: LinearOperator instance has no attribute 'dtype'
A=sparse.scipy.linalg.linearoperator（shape=（n，n），matvec=mv，dtype='float64'）
回溯（最近一次呼叫最后一次）：
文件“”，第1行，在
文件“/usr/lib/python2.7/dist-packages/scipy/sparse/linalg/isolve/iterative.py”，第393行，在gmres中
A、 M，x，b，后处理=生成系统（A，M，x0，b，xtype）
文件“/usr/lib/python2.7/dist packages/scipy/sparse/linalg/isolve/utils.py”，第119行，在make_系统中
M=线性构造器（A.shape，matvec=psolve，rmatvec=rpsolve，dtype=A.dtype）
AttributeError:LineRoperator实例没有属性“dtype”
LinearOperatordtypedtype='float64'In [1]: import numpy as np
In [2]: from scipy.sparse.linalg import LinearOperator, gmres
In [3]: def mymatvec(v):
...: a = np.array([[4,2,1],[2,2,1],[1,1,1]])
...: return a.dot(v)
...:
In [4]: A = LinearOperator((3,3), mymatvec, dtype='float64')
In [5]: b = np.array([1,2,3])
In [6]: gmres(A, b)
Out[6]: (array([-0.5, -0.5, 4. ]), 0)
您使用稀疏模块进行密集矩阵运算有什么特别的原因吗？不，我不想，但我只能在稀疏库中找到它。有一个scipy.linalg.gmres，但它似乎要么被弃用，要么没有出现在我的scipy版本中。文档声称它适用于稀疏、密集或线性构造器类型的矩阵。请注意，不需要使用线性构造器，您可以直接将密集矩阵传递给
gmres（a，b）