*** a/fipy/matrices/pysparseMatrix.py	2013-09-30 19:17:46.000000000 -0700
--- b/fipy/matrices/pysparseMatrix.py	2014-01-16 10:55:13.000000000 -0800
*************** class _PysparseMatrix(_SparseMatrix):
*** 165,170 ****
--- 165,175 ----
                  y = numerix.empty((self.matrix.shape[0],))
                  self.matrix.matvec(other, y)
                  return y
+             elif shape[0]*shape[1] == N:
+                 other = other.ravel()
+                 y = numerix.empty((self.matrix.shape[0],))
+                 self.matrix.matvec(other, y)
+                 return y
              else:
                  raise TypeError