2

coo_matrix我最近在使用 将 a 转换为密集矩阵时遇到了困难scipy。我有一个 dtypefloat16稀疏矩阵并尝试将其转换为密集矩阵。该错误抱怨被赋予了一个 type 的数组char。但是,我很确定我正在传递一个 type 数组float16

错误是:

    self.Xd_train = X_train.todense()
File "C:\Python27\lib\site-packages\scipy\sparse\base.py", line 501, in todense
  return np.asmatrix(self.toarray(order=order, out=out))
File "C:\Python27\lib\site-packages\scipy\sparse\coo.py", line 241, in toarray
  B.ravel('A'), fortran)
File "C:\Python27\lib\site-packages\scipy\sparse\sparsetools\coo.py", line 175, in coo_todense
  return _coo.coo_todense(*args)
TypeError: Array of type 'float' required.  Array of type 'char' given

错误出现在类构造函数中:

self.Xd_train = X_train.todense()

该矩阵X_train似乎是格式良好的,绝对不是类型char

>> X_train.dtype
float16

>> X_train.shape
(6206, 4712)

>> type(X_train)
<class 'scipy.sparse.coo.coo_matrix'>

>> str(X_train)

(0, 63)       2.0
(0, 72)       1.0
(0, 76)       2.0
(0, 100)      1.0
(0, 104)      1.0
(0, 5)        1.0
(0, 10)       2.0
(0, 134)      2.0
(0, 20)       3.0
(0, 263)      1.0
(0, 264)      1.0
(0, 265)      1.0
(0, 27)       1.0
(0, 148)      2.0
(0, 32)       1.0
(0, 275)      1.0
(0, 35)       1.0
(0, 36)       1.0
(0, 279)      1.0
(0, 39)       1.0
(0, 41)       1.0
(0, 42)       1.0
(0, 52)       1.0
(0, 59)       4.0
(1, 72)       1.0
:     :
(6205, 133)   1.0
(6205, 134)   4.0
(6205, 135)   4.0
(6205, 136)   2.0
(6205, 137)   6.0
(6205, 138)   1.0
(6205, 139)   4.0
(6205, 20)    4.0
(6205, 142)   4.0
(6205, 23)    2.0
(6205, 24)    2.0
(6205, 26)    2.0
(6205, 27)    2.0
(6205, 32)    1.0
(6205, 33)    1.0
(6205, 35)    1.0
(6205, 36)    1.0
(6205, 37)    1.0
(6205, 39)    1.0
(6205, 40)    1.0
(6205, 41)    1.0
(6205, 42)    1.0
(6205, 43)    1.0
(6205, 56)    3.0
(6205, 60)    1.0

有什么想法可能是什么问题?另外,如果需要其他详细信息/信息,请告诉我。

我在 Windows 7 上使用 Python 2.7.2,Numpy 1.7 和 Scipy 0.11。谢谢。

4

1 回答 1

2

此错误也发生在最新的 scipy master 分支中。例如

>>> coo_matrix([[0]], dtype=np.float16).todense()

引发相同的异常。数据类型np.float16相对较新,scipy(可能还有其他地方)中有很多代码尚未经过测试。

如果您将稀疏矩阵更改为np.float32,它应该可以工作。

我在 scipy github 网站上为此创建了一个问题:https ://github.com/scipy/scipy/issues/2481

于 2013-05-11T00:43:15.527 回答