Supplied block_diag function for sparse matrices to overcome scipy version limitation.