python将numpy矩阵延对角线方向扩展

在做数学运算的时候经常碰到这样的需求：
给定矩阵A，希望得到一个对角线方向上增广的矩阵
B =
| A 0 0 |
| 0 A 0 |
| 0 0 A …|

查了很久，没有看到numpy原生的函数可以实现，只能自己写。
朴素的方法就是写循环，不停stack加零，然后再赋值。不过这样做就少了一些python代码的优雅。
一种可能的写法：

def diag_stack(A : np.ndarray, n : int=2):
  return np.tile(A, (n,n)) * \
         np.repeat(np.repeat(np.eye(n), A.shape[0], axis=0), A.shape[1], axis=1).astype(int)

本质上其实是np.tile和np.repeat的拼接

测试一下：

a = np.array([[1,2,3],
              [4,5,6]])
print(diag_stack(a, 3))

输出：

[[1 2 3 0 0 0 0 0 0]
 [4 5 6 0 0 0 0 0 0]
 [0 0 0 1 2 3 0 0 0]
 [0 0 0 4 5 6 0 0 0]
 [0 0 0 0 0 0 1 2 3]
 [0 0 0 0 0 0 4 5 6]]

优点是代码简单。缺点是算力消耗不一定是最优。
如有更简便的方法，欢迎分享。