Added function block_diag, since it is not yet available in SciPy 0.7.x (upgrade difficult on ATLAS cluster)