提取R中的块对角矩阵

Question

我想要一个只有大矩阵D的块对角线的稀疏矩阵

目前，我的操作如下：

getting-started

问题在于D很大，并且此元素明智的乘法需要大量时间。有没有更好的方法可以做到这一点。