[SciPy-dev] Sparse matvec inconsistency

[Bart Vandereycken]

Hi all,

the matvec (A*x) routine of sparse shows an inconsistency: the shape of 
the returned vector is not the same as the vector x when shape==(n,1). A 
simple change in the _matvec routines in sparse.py should do the trick...

BTW I would have posted a ticket for this, but I lost my password for my 
account... Maybe someone can make a lost password link on the trac?

-- bart


import scipy as SY
import scipy.sparse as SP

A = SP.spidentity(10)

x = SY.ones((10,1))
y = A*x
# returns (10,)
print y.shape

x = SY.ones((10,2))
y = A*x
# returns (10,2)
print y.shape