Manera eficiente de normalizar una matriz escasa Scipy

Me gustaría escribir una función que normalice las filas de una matriz dispersa grande (de manera que sumen una).

from pylab import *
import scipy.sparse as sp

def normalize(W):
    z = W.sum(0)
    z[z < 1e-6] = 1e-6
    return W / z[None,:]

w = (rand(10,10)<0.1)*rand(10,10)
w = sp.csr_matrix(w)
w = normalize(w)

Sin embargo, esto da la siguiente excepción:

File "/usr/lib/python2.6/dist-packages/scipy/sparse/base.py", line 325, in __div__
     return self.__truediv__(other)
File "/usr/lib/python2.6/dist-packages/scipy/sparse/compressed.py", line 230, in  __truediv__
   raise NotImplementedError

¿Hay soluciones razonablemente simples? He miradoesta, pero todavía no tengo claro cómo hacer realmente la división.