BF: scipy 0.12.0 rightfully uses dia_matrix but doesn't provide functional setdiag for it -- use 'csr' matrix instead

https://github.com/scipy/scipy/issues/2931