Apologies if this has been addressed since my SciPy v0.13.3, but inputting empty arrays into linalg.block_diag yields nonintuitive outputs (Python v2.7.6): >>> scipy.linalg.block_diag([[1,0],[0,1]]...