如何在numpy中重现scipy.linalg.block_diagonal？

Question

我有一个使用的代码

import numpy, scipy
r = numpy.random.random((3,2,2))
d = scipy.linalg.block_diagonal(*r)

我希望消除对scipy的依赖。在numpy中创建块对角矩阵的最典型方法是什么，它涵盖所有块都具有相同大小的情况（在上面的示例中为2 * 2）？

（我确实调查过numpy.diag，但它没有用于对角线的功能...）