Question

我在scipy.sparse.csr_matrix中发现了一个意想不到的行为，这对我来说似乎是个错误。谁能证实这不正常？我不是稀疏结构的专家，所以我可能误解了正确的用法。

>>> import scipy.sparse
>>> a=scipy.sparse.csr_matrix((1,1))
>>> b=scipy.sparse.csr_matrix((1,1))
>>> b[0,0]=1
/home/marco/anaconda3/envs/py35/lib/python3.5/site-packages/scipy/sparse/compressed.py:730: SparseEfficiencyWarning: Changing the sparsity structure of a csr_matrix is expensive. lil_matrix is more efficient.
  SparseEfficiencyWarning)
>>> a/b
matrix([[ nan]])

另一方面，numpy正确处理了这个问题：

>>> import numpy as np
>>> a=np.zeros((1,1))
>>> b=np.ones((1,1))
>>> a/b
array([[ 0.]])

由于

Answer 1

对于稀疏矩阵/稀疏矩阵，

SciPy的/稀疏/ compressed.py

    if np.issubdtype(r.dtype, np.inexact):
        # Eldiv leaves entries outside the combined sparsity
        # pattern empty, so they must be filled manually. They are
        # always nan, so that the matrix is completely full.
        out = np.empty(self.shape, dtype=self.dtype)
        out.fill(np.nan)
        r = r.tocoo()
        out[r.row, r.col] = r.data
        out = np.matrix(out)

本节将解释该行动。

尝试使用稍大的矩阵

In [69]: a=sparse.csr_matrix([[1.,0],[0,1]])
In [70]: b=sparse.csr_matrix([[1.,1],[0,1]])
In [72]: (a/b)
Out[72]: 
matrix([[  1.,  nan],
        [ nan,   1.]])

因此a有0（没有稀疏值），除法为nan。它返回一个密集矩阵，并填入nan。

如果没有这段代码，按元素划分的稀疏元素会生成一个稀疏矩阵，其中包含那些“空”对角线槽。

In [73]: a._binopt(b,'_eldiv_')
Out[73]: 
<2x2 sparse matrix of type '<class 'numpy.float64'>'
    with 2 stored elements in Compressed Sparse Row format>
In [74]: a._binopt(b,'_eldiv_').A
Out[74]: 
array([[ 1.,  0.],
       [ 0.,  1.]])

反向可能是有益的

In [76]: b/a
Out[76]: 
matrix([[  1.,  inf],
        [ nan,   1.]])
In [77]: b._binopt(a,'_eldiv_').A
Out[77]: 
array([[  1.,  inf],
       [  0.,   1.]])

看起来combined sparsity pattern由分子决定。在eliminate_zeros之后，进一步的测试看起来像这样。

In [138]: a1=sparse.csr_matrix(np.ones((2,2)))
In [139]: a1
Out[139]: 
<2x2 sparse matrix of type '<class 'numpy.float64'>'
    with 4 stored elements in Compressed Sparse Row format>
In [140]: a1[0,1]=0
In [141]: a1
Out[141]: 
<2x2 sparse matrix of type '<class 'numpy.float64'>'
    with 4 stored elements in Compressed Sparse Row format>
In [142]: a1/b
Out[142]: 
matrix([[  1.,  nan],
        [ inf,   1.]])

Scipy稀疏csr矩阵在0.0 / 1.0上返回nan

1 个答案: